Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 10,621 A→T 5.1% intergenic (‑127/+22) satP ← / ← yaaW succinate‑acetate transporter/UPF0174 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091310,6210AT 5.1% 345.8 / 5.4 116intergenic (‑127/+22)satP/yaaWsuccinate‑acetate transporter/UPF0174 family protein
Reads supporting (aligned to +/- strand):  ref base A (52/57);  new base T (2/4);  total (54/62)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 128 total aligned reads displayed.

TTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCC  >  NC_000913/10472‑10769
                                                                                                                                                     |                                                                                                                                                    
ttAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGgaaaatgaa                                                                                                                                                      <  2:1524941/150‑1 (MQ=255)
 tAGCCAACTTAGTGTTGCCCATAATTCATCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAAATGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTATGAGGGgaaaatgaaa                                                                                                                                                     >  1:1785966/1‑150 (MQ=255)
  aGCCAACTTAGTGCTGCGCATAAAACCTCTCCACGCCTCATCGACTGAATCGTGACAAACTTTCCCTGACTACCTGTAGCGTATTCAGGGCGCGGCCTACTCATCAGCCAGTGGCGCAGTGTCTACGATCTTTTGAGGGgaaaatgaaaa                                                                                                                                                    <  2:1821538/150‑1 (MQ=255)
   gCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTTAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGgaaaatgaaaat                                                                                                                                                   >  1:439538/1‑150 (MQ=255)
   gCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCTGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGCGGGgaaaatgaaaat                                                                                                                                                   >  2:1727304/1‑150 (MQ=255)
   gCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGgaaaatgaaaat                                                                                                                                                   >  2:808144/1‑150 (MQ=255)
     cAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAAttt                                                                                                                                                 >  1:1578572/1‑150 (MQ=255)
     cAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAAttt                                                                                                                                                 >  1:1871892/1‑150 (MQ=255)
        cTTATTGTTGCCCATAATTCCTCAAAAATCATAATCGAATGAATCGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTcc                                                                                                                                              <  1:107241/150‑1 (MQ=255)
        cTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTcc                                                                                                                                              <  1:1021473/150‑1 (MQ=255)
          tAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGGTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATTTTCATTTTcccc                                                                                                                                            <  1:999403/150‑1 (MQ=255)
            gtgtTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCgg                                                                                                                                          >  1:1288665/1‑150 (MQ=255)
            gtgtTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACCGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATTTTCATTTTcccctc                                                                                                                                          <  2:1714121/150‑3 (MQ=255)
              gtTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATTATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGtt                                                                                                                                        >  2:1622208/1‑150 (MQ=255)
                   ccATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCgg                                                                                                                                   <  2:1288665/150‑1 (MQ=255)
                     aTAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTa                                                                                                                                 <  1:20132/150‑1 (MQ=255)
                       aaTTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGAATc                                                                                                                               <  1:1394424/150‑1 (MQ=255)
                           ccTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGAc                                                                                                                           <  2:550227/150‑1 (MQ=255)
                           ccTCAAAAATCATCATCGAATCAATGGAGAAAAAAGTTCCCTGAATAACGGCCCTGTTTTAAGGGCGAGGCAAAATAATCAGCCAGTGGGGCAGTGGCAACGAGCTTTTGAGGCGAAACTGAAAATTGTCCCCGGTTTCCGGTATCAGAc                                                                                                                           <  2:319603/150‑1 (MQ=255)
                            cTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAATATGAAAATTTTCCCCGGTTTCCGGTATCAGAcc                                                                                                                          <  1:15047/150‑1 (MQ=255)
                            cTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGAcc                                                                                                                          >  1:978266/1‑150 (MQ=255)
                            cTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGAcc                                                                                                                          >  2:70495/1‑150 (MQ=255)
                                aaaaTCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAg                                                                                                                      <  2:1369224/150‑1 (MQ=255)
                                    tcatcatcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGgc                                                                                                                  <  1:905937/150‑1 (MQ=255)
                                    tcatcatcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGgc                                                                                                                  <  2:649571/150‑1 (MQ=255)
                                          tcGAATGAATGGTGAAATAATTTCACTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAAc                                                                                                            <  2:1362527/150‑1 (MQ=255)
                                           cGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAAcc                                                                                                           >  2:865479/1‑150 (MQ=255)
                                               tgaatgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATcc                                                                                                       >  1:52161/1‑150 (MQ=255)
                                                gaatgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCg                                                                                                      >  2:1810395/1‑150 (MQ=255)
                                                   tgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGgcg                                                                                                   >  1:916233/1‑150 (MQ=255)
                                                   tgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGgcg                                                                                                   >  1:1273884/1‑150 (MQ=255)
                                                   tgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGgcg                                                                                                   <  1:996267/150‑1 (MQ=255)
                                                      tGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgcag                                                                                                >  2:1032175/1‑150 (MQ=255)
                                                       gAAATAATTTCCCTGAATACCTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgcagg                                                                                               <  1:942229/150‑1 (MQ=255)
                                                       gAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGTGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgcagg                                                                                               >  1:225264/1‑150 (MQ=255)
                                                       gAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGAATCAGACCTGAGTGGCGCTAACCATCCGGCgcagg                                                                                               <  1:1469951/150‑1 (MQ=255)
                                                        aaaTAATTTCCCTGAATAACTGTAGTGTTATCAGGGCGCTGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCACCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCtcaggc                                                                                              >  1:1162774/1‑150 (MQ=255)
                                                           taatTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGcaggcagg                                                                                           <  1:1688410/150‑1 (MQ=255)
                                                             atTTCCCTGAATAACTGTAGTGTTATCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCg                                                                                         <  2:1362319/150‑1 (MQ=255)
                                                                tCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCTAtt                                                                                      <  2:1780747/150‑1 (MQ=255)
                                                                 cccTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGAttt                                                                                     >  2:960663/1‑150 (MQ=255)
                                                                  ccTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTg                                                                                    >  2:1621072/1‑150 (MQ=255)
                                                                  ccTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTg                                                                                    >  2:858782/1‑150 (MQ=255)
                                                                  ccTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTg                                                                                    <  2:1016168/150‑1 (MQ=255)
                                                                  ccTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTg                                                                                    <  2:1610434/150‑1 (MQ=255)
                                                                  ccTGAATAACTGTAGTGAGTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTg                                                                                    <  1:1318979/150‑1 (MQ=255)
                                                                   cTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCGACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGAATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGc                                                                                   <  1:959282/150‑1 (MQ=255)
                                                                    tGACTAACTGTCGTGTTATCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCa                                                                                  <  1:588281/150‑1 (MQ=255)
                                                                    tGAATAACTGTAGTGTTTTCAGGGCGCGGCATACTAATCAGCCAGTGGGGCAGTGTCTACGCTCGTTTGAGGGGAAACTGACAATTTTCCCCGGTGTCCGGTATCAGACCTGAGGGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCa                                                                                  <  1:120357/150‑1 (MQ=255)
                                                                     gAATAACTGTAGTGTTTTCAGGGCGCGGCATAATTATCACGCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAAATTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAg                                                                                 <  2:1163478/150‑1 (MQ=255)
                                                                     gAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTATACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAg                                                                                 <  2:1622568/150‑1 (MQ=255)
                                                                      aaTAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGt                                                                                <  2:18774/150‑1 (MQ=255)
                                                                          aCTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACgg                                                                            <  1:1285345/150‑1 (MQ=255)
                                                                           cTGTAGCGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGc                                                                           <  2:444708/150‑1 (MQ=255)
                                                                             gTAGTGTTTTCAGGGCGCGGCATAATTATCAGCCAGTGGGGGAGTGTCTACGAGCTTTTGAGGGGAAAATGAAAATTTTCGCCGGTTTCCGGTATCTGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTg                                                                         <  1:976534/150‑1 (MQ=255)
                                                                             gTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGTGCAGTCAGGCGATTTGCAGTACGGCTg                                                                         >  2:1019304/1‑150 (MQ=255)
                                                                               aGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGa                                                                       >  2:1792956/1‑150 (MQ=255)
                                                                                gtgtTTTCAGGGCGCCGCAGAATAATCAGCCCGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGaa                                                                      <  1:1213160/150‑1 (MQ=255)
                                                                                 tgtTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAAt                                                                     >  1:123313/1‑150 (MQ=255)
                                                                                 tgtTTTCAGGGAGCGGCATAATAATCAGCCAGTGGGGCATTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAAt                                                                     <  1:1650517/150‑1 (MQ=255)
                                                                                  gtTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATc                                                                    <  2:332782/150‑1 (MQ=255)
                                                                                       cAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGCTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGTCAGGCGATTTGCAGTACGGCTGGAATCGTCAc                                                               >  2:1068728/1‑150 (MQ=255)
                                                                                        aGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCAcg                                                              >  1:1410344/1‑150 (MQ=255)
                                                                                        aGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCAcg                                                              >  2:1250612/1‑150 (MQ=255)
                                                                                             gcgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATa                                                         >  2:574766/1‑150 (MQ=255)
                                                                                                gCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCCGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGgc                                                      >  2:1155956/1‑150 (MQ=255)
                                                                                                    aatcatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTg                                                  <  1:1244987/150‑1 (MQ=255)
                                                                                                     ataatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGCATCGTCACGCGATAGGCGCTGc                                                 >  2:54052/1‑150 (MQ=255)
                                                                                                       aatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCg                                               >  2:1364936/1‑150 (MQ=255)
                                                                                                           aGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGa                                           >  1:230148/1‑150 (MQ=255)
                                                                                                            gCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGAc                                          <  1:1453564/150‑1 (MQ=255)
                                                                                                             ccAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGAcc                                         >  1:598761/1‑150 (MQ=255)
                                                                                                              cAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGAGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCg                                        >  1:38649/1‑150 (MQ=255)
                                                                                                                  ggggCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCttt                                    <  1:472011/150‑1 (MQ=255)
                                                                                                                     gCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCAGCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAc                                 >  2:1531102/1‑150 (MQ=255)
                                                                                                                     gCAGTGTCTACGATCTTTGGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAc                                 <  2:326689/150‑1 (MQ=255)
                                                                                                                          gtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAt                            >  1:1174087/1‑150 (MQ=255)
                                                                                                                            cTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                          >  1:1227880/1‑150 (MQ=255)
                                                                                                                            cTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                          >  2:960124/1‑150 (MQ=255)
                                                                                                                            cTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                          >  1:1088011/1‑150 (MQ=255)
                                                                                                                             tACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTa                         >  1:766481/1‑150 (MQ=255)
                                                                                                                             tACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTa                         >  1:430759/1‑150 (MQ=255)
                                                                                                                             tACGATCATTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCTGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTa                         <  1:1790134/150‑1 (MQ=255)
                                                                                                                                gATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTg                      <  2:439538/150‑1 (MQ=255)
                                                                                                                                gATCTTTTGAGGGGAAAATGAAAATTTTCACCGGTTTCCGGTATCAGACCTGAGTGGCGATAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTg                      >  2:1405961/1‑150 (MQ=255)
                                                                                                                                   cTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCAGCTGACCGCTTTAACCCCATTTAGTGCCg                   <  1:1622208/150‑1 (MQ=255)
                                                                                                                                    ttttGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGc                  <  1:70495/150‑1 (MQ=255)
                                                                                                                                     tttGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCa                 <  2:1398155/150‑1 (MQ=255)
                                                                                                                                       tGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGAGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCAcc               >  2:1722844/1‑150 (MQ=255)
                                                                                                                                         cggggAAAATTTTCATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTa             >  2:83663/2‑150 (MQ=255)
                                                                                                                                         cggggAAAATTTTCATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGAGAGTTGCAGTACGGCTGGAATCGTCACGCGATAGCCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTa             >  1:1228457/2‑150 (MQ=255)
                                                                                                                                          ggggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGCCCGCTTTAACCCCATTTAGTGCCGCACCTAc            >  2:1810389/1‑150 (MQ=255)
                                                                                                                                          ggggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc            <  2:123313/150‑1 (MQ=255)
                                                                                                                                           gggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACa           >  2:1752628/1‑150 (MQ=255)
                                                                                                                                           gggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCGAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCGGGAAGCGGCACGCGATAGGCGCTGCCGCGGACCGCTGTACCCCCATTTAGTGCCGAACCTAcc           >  2:1745437/1‑149 (MQ=255)
                                                                                                                                            ggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAg          >  1:1099686/1‑150 (MQ=255)
                                                                                                                                             gaaaatgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAgg         <  2:1578572/150‑1 (MQ=255)
                                                                                                                                              aaaatgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAggg        <  1:1810395/150‑1 (MQ=255)
                                                                                                                                                   gaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCTCCTACAGGGCCTcc   >  2:1004924/1‑150 (MQ=255)
                                                                                                                                                    aaaatTTTCCCCGGTTTCCGGGATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTccc  <  2:950768/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCC  >  NC_000913/10472‑10769

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑