Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,855 T→G 5.6% intergenic (+72/+80) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8550TG 5.6% 405.6 / 2.1 169intergenic (+72/+80)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (63/92);  new base G (4/9);  total (68/101)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

CAAACCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAAT  >  NC_000913/3504708‑3505003
                                                                                                                                                   |                                                                                                                                                    
 aaaCCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTcc                                                                                                                                                   >  2:931707/1‑150 (MQ=255)
 aaaCCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTGcc                                                                                                                                                   >  2:2137279/1‑150 (MQ=255)
  aaCCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTccc                                                                                                                                                  <  2:1408052/150‑1 (MQ=255)
  aaCCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTccc                                                                                                                                                  >  2:876925/1‑150 (MQ=255)
  aaCCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGATccc                                                                                                                                                  >  2:3530002/1‑150 (MQ=255)
      ggTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCa                                                                                                                                              <  1:1500390/150‑1 (MQ=255)
      ggTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCa                                                                                                                                              <  2:883709/150‑1 (MQ=255)
            cGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTaa                                                                                                                                        <  1:1161658/150‑1 (MQ=255)
              tctcCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAAcc                                                                                                                                      >  2:278571/1‑150 (MQ=255)
              gctcCGAACTCTTCTGCCGCGCCACTCGCATCACCTTAACTATTGATAATACAAGACATTCACCGTAGGCCGGATCAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAAGTTTACATCCCGGACGTGCCCCGCccccgaccc                                                                                                                                      >  1:2231317/2‑145 (MQ=38)
                tcCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCct                                                                                                                                    <  1:1640605/150‑1 (MQ=255)
                  cGATCTCTTCTGCCGCGCCAATCGCATCACCTTCACTATTGATAATAACAGACAGTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCTTCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCctct                                                                                                                                  <  2:2345964/150‑1 (MQ=39)
                   gAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCtctc                                                                                                                                 <  1:3898724/150‑1 (MQ=255)
                         ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                           <  1:2721042/150‑1 (MQ=255)
                         ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                           <  1:2509828/150‑1 (MQ=255)
                         ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAACCCACTCCCCaaa                                                                                                                           >  1:3552457/1‑150 (MQ=255)
                             gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAGCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGACCAAACCCAAAACCTCTTACCAAAgggg                                                                                                                       >  2:1102293/1‑150 (MQ=38)
                             gCCGCGACAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAgggg                                                                                                                       <  2:1429501/150‑1 (MQ=255)
                              ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGc                                                                                                                      <  2:3308938/150‑1 (MQ=255)
                              ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGAAACCAAAACCACTCCCCAAAGGGGc                                                                                                                      >  1:1681878/1‑150 (MQ=39)
                                gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAACCCGCACCCCCAAGGGGCgc                                                                                                                    >  2:653672/1‑149 (MQ=39)
                                 cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                   <  2:3012821/150‑1 (MQ=255)
                                 cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAAACCTCTCCCCAAAGGGGCGAg                                                                                                                   <  2:1478251/150‑1 (MQ=255)
                                 cgcCAATCGCATCAACTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCACTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                   <  2:1696296/150‑1 (MQ=255)
                                  gcCAATCGCATCACCTTAACTATTTATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                  <  1:2816919/150‑1 (MQ=255)
                                   ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                 <  2:3882170/150‑1 (MQ=255)
                                     aagcGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGa                                                                                                               <  1:1854590/147‑1 (MQ=255)
                                     aaTCGCATCACCTTACCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                               <  1:1422416/150‑1 (MQ=255)
                                     aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                               <  1:2157281/150‑1 (MQ=255)
                                     aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                               >  1:3403260/1‑150 (MQ=255)
                                       tCGCATCACCTTCACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGAcc                                                                                                             <  1:3937037/150‑1 (MQ=255)
                                       tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTCACCCTCTCCCCAAAGGGGCGAGGGGAcc                                                                                                             <  1:2320482/150‑1 (MQ=255)
                                        cGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACCCCCAAACCCACTCCCCAAACGGGCGCGGGGAccc                                                                                                            >  1:3882187/1‑149 (MQ=38)
                                          cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACGACCGAAACCGAAACCCCAAGGGGCGGGGGGGcccgg                                                                                                          >  1:3115293/1‑147 (MQ=25)
                                          cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGtt                                                                                                          <  1:2586885/150‑1 (MQ=255)
                                          cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCGAACCCACACCCCCCACGGGCGCGGGGCCCggt                                                                                                          >  1:2240348/1‑148 (MQ=38)
                                          cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCACAACCCAAAGGGGCGAGGGGGcccgt                                                                                                          >  1:99801/1‑147 (MQ=38)
                                          cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTGCCCCTCACCCGAACCCACACCCCCAAGGGGCGAGGGGAcccgt                                                                                                          >  1:3544425/1‑147 (MQ=39)
                                          cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCACACCCAAACCCACACCCCCAACGGGCGAGGGGccccga                                                                                                          >  2:3518087/1‑147 (MQ=38)
                                            tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                        <  2:1789245/150‑1 (MQ=255)
                                            tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                        <  2:1789261/150‑1 (MQ=255)
                                            tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                        <  2:2618536/150‑1 (MQ=255)
                                            tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                        <  1:2098352/150‑1 (MQ=255)
                                              aCCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGt                                                                                                      <  1:1069778/150‑1 (MQ=255)
                                              aCCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGt                                                                                                      <  1:1069770/150‑1 (MQ=255)
                                              aCCTTAACTATTGATAATAAAAGACATTAACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCGCACCCAAACCCGCACCCCCCACGGGCGCGGGGCCCgggccgg                                                                                                      >  2:3857525/1‑144 (MQ=37)
                                                cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGGTGAATTTTAAATCCCGTACGGTGCCCTCACCCTAAACCGCTCCCCCAAGTGGCGGGGGGACCCGTCAGTAc                                                                                                    >  1:466639/1‑150 (MQ=37)
                                                cTTAACTATTGATAATAAAAGACAGTACCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTAc                                                                                                    <  2:1013700/150‑1 (MQ=255)
                                                 ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACg                                                                                                   <  1:212256/150‑1 (MQ=255)
                                                 ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCACTACCCCAAGGGGCGAGGGGCCCCGTCAGTACg                                                                                                   >  2:2895652/1‑150 (MQ=38)
                                                  tAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACCGGGGCGAGGGGACCGTTCAGTACgg                                                                                                  <  2:3071225/150‑1 (MQ=255)
                                                  tAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACgg                                                                                                  <  1:2570355/150‑1 (MQ=255)
                                                  tAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACgg                                                                                                  <  1:640685/150‑1 (MQ=255)
                                                   aaCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTTCGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGa                                                                                                 <  1:2987956/150‑1 (MQ=255)
                                                      tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAAtt                                                                                              <  2:2941724/150‑1 (MQ=255)
                                                       aTTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAAttt                                                                                             <  2:2512858/150‑1 (MQ=255)
                                                        ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCACGGGGCGAGGGGACCGTTCAGTACGGAAtttn                                                                                            <  1:2267518/150‑2 (MQ=255)
                                                         tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                           <  1:2931601/150‑1 (MQ=255)
                                                          gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGc                                                                                          <  2:3176303/150‑1 (MQ=255)
                                                           ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCACCCTAAACCACACCCCCAAGGGGCGAGGGGACCGTTAAGTACGGAATTTAgat                                                                                         >  1:637320/1‑148 (MQ=38)
                                                           ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCACCCAAAAACACACCCCAAAAGGGCGAGGGGGCCCGTAAGTACCGAATTTATCt                                                                                         >  2:747981/1‑150 (MQ=25)
                                                                 aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCGCTCCCCCAAGGGGCGCGGGGACCGGTCAGTACGGAATTTAGCTTGACgg                                                                                   >  2:3683491/1‑150 (MQ=255)
                                                                  aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGt                                                                                  <  1:832771/150‑1 (MQ=255)
                                                                   aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                 <  1:739386/150‑1 (MQ=255)
                                                                   aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                 <  2:2057025/150‑1 (MQ=255)
                                                                           aCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATACCGTACGTTGCCCACAACCAAAACCAAAACCCAAAGGGGCCAGGGGGCCCGTCACTACAGAAATTAACATTTCGGTGcgaaaaac                                                                         >  1:2858727/1‑143 (MQ=16)
                                                                            ccGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGc                                                                        <  1:1937836/150‑1 (MQ=255)
                                                                            ccGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGc                                                                        <  2:2346419/150‑1 (MQ=255)
                                                                             cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                       <  1:189963/150‑1 (MQ=255)
                                                                             cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                       <  2:1509939/150‑1 (MQ=255)
                                                                                aGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTa                                                                    <  1:3315821/150‑1 (MQ=255)
                                                                                     ggATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATAc                                                               <  2:4054102/150‑1 (MQ=255)
                                                                                      gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCACACCCCAAAGGGGCGAGGGGACCCGTCACTAAGGAAGTTAGCATTACCGTGCGATCCACCAGTAAATAca                                                              >  2:55533/1‑149 (MQ=37)
                                                                                       aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCTTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACACTCTCCCCAAAGGGTCGAGGGGACCGTTCATTACGGAATTTAGCTTGACGTTGCTATCCAGCATTAAATACTg                                                             <  2:902773/150‑1 (MQ=255)
                                                                                       aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                             <  1:3561330/150‑1 (MQ=255)
                                                                                       aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                             <  2:1239972/150‑1 (MQ=255)
                                                                                       aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                             <  1:1674260/150‑1 (MQ=255)
                                                                                       aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGATGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                             <  1:4038932/150‑1 (MQ=255)
                                                                                         aaGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAAACCACAACCCAAAGGGGCGAGGGGACCCTTCAGTACGGAAATTATCTTGAAGGTGCTATACAGCAATAAAAACAGGc                                                           >  2:2438078/1‑150 (MQ=255)
                                                                                             tGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATa                                                       <  2:1991587/150‑1 (MQ=255)
                                                                                               cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                     <  2:1232009/150‑1 (MQ=255)
                                                                                                  cAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACtctc                                                  <  2:1480396/150‑1 (MQ=255)
                                                                                                   aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                 <  2:409750/150‑1 (MQ=255)
                                                                                                   aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                 <  1:2004774/150‑1 (MQ=255)
                                                                                                    gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                <  2:62857/150‑1 (MQ=255)
                                                                                                    gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAACCCTCTCCCCAAAGGGGCGCGGGGACCGGTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATAATGGCCATACTCTCAc                                                >  2:3730863/1‑150 (MQ=39)
                                                                                                          catcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTAACCCTAAACCTCTACCCAAATGGGGGAGGGGACCGTTCAGTAAGGAATTTTGCTTTAAGGTGCTATCCAGCATTAAATAATGGCAATAATATAAAGCaaaa                                          >  2:327454/1‑150 (MQ=255)
                                                                                                              cGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCg                                      <  2:1942884/150‑1 (MQ=255)
                                                                                                              cGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCg                                      <  2:2486257/150‑1 (MQ=255)
                                                                                                               ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCAACCTAAACCACTCCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCATGACGGTGCTATCCAGCATTAAATAATGGCAATACTCTCACGCAAAATACgc                                     >  2:2043986/1‑150 (MQ=255)
                                                                                                                gCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgca                                    <  2:3347697/150‑1 (MQ=255)
                                                                                                                  gATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCACTACCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCATGACGGTGCTATCCAGCATTAAAAACTGGCAATACTCTCACGCAAAAACCgcagc                                  >  2:3991836/1‑150 (MQ=39)
                                                                                                                  gATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAAACCTCTCCCCAAAGGGGCGAGGGGACCGGTCAGTAAGGAAATTAGCTTGACGGTGCTATCCAGCATTAAATAATGGCAATAAAATCAAGCAAAAACCgcagc                                  >  2:2245145/1‑150 (MQ=38)
                                                                                                                          cgcgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                          <  1:418033/150‑1 (MQ=255)
                                                                                                                            atgagAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCtt                        <  1:2199460/146‑1 (MQ=255)
                                                                                                                               tGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTACCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTTGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                     >  2:3421815/1‑150 (MQ=255)
                                                                                                                                  aTTTTACATCCCGTACGTTCCCCCCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGc                  <  1:2335405/150‑1 (MQ=255)
                                                                                                                                   ttttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCt                 <  2:127286/150‑1 (MQ=255)
                                                                                                                                   ttttACATCCCGTACGTTCCCCTCACCCGAACCCACACCCCAAAGGGGCGAGGGGACCGTTAAGTACGGAATTTAGCTTGACGGTGCTATACAGCAATAAATACTGGCAATACACTAAAGCAAAATCCCCAGCACCCTCTCGTCAACGcg                 >  2:1715820/1‑149 (MQ=37)
                                                                                                                                               tACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTAAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGtt     >  2:2000381/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAAACCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAAT  >  NC_000913/3504708‑3505003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: