Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,011,926 G→A 6.2% intergenic (+55/‑75) uup → / → pqiA ATP‑binding protein with possible role in replication/intermembrane transport protein PqiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,011,9260GA 6.2% 361.3 / 7.7 113intergenic (+55/‑75)uup/pqiAATP‑binding protein with possible role in replication/intermembrane transport protein PqiA
Reads supporting (aligned to +/- strand):  ref base G (60/46);  new base A (3/4);  total (63/50)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 116 total aligned reads displayed.

AAAAAGTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCT  >  NC_000913/1011777‑1012074
                                                                                                                                                     |                                                                                                                                                    
aaaaaGTGCTTGCTGATATGGCTTCTGCGGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGTAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGGAAAATTGCCtgatg                                                                                                                                                      <  2:1122355/150‑1 (MQ=255)
    aGTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACt                                                                                                                                                  <  2:749922/150‑1 (MQ=255)
     gTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTg                                                                                                                                                 >  1:884104/1‑150 (MQ=255)
       gcttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTggg                                                                                                                                               >  1:687366/1‑150 (MQ=255)
         ttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGtt                                                                                                                                             >  1:458609/1‑150 (MQ=255)
         ttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGtt                                                                                                                                             >  1:914104/1‑150 (MQ=255)
         ttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGtt                                                                                                                                             >  2:1158646/1‑150 (MQ=255)
         ttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGtt                                                                                                                                             <  2:1149634/150‑1 (MQ=255)
          tgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGttt                                                                                                                                            >  1:109898/1‑150 (MQ=255)
           gctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTa                                                                                                                                           >  1:1462095/1‑150 (MQ=255)
           gctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTa                                                                                                                                           >  1:151831/1‑150 (MQ=255)
             tGATATGGCTGCTGAAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATc                                                                                                                                         >  1:1395970/1‑150 (MQ=255)
              gATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCa                                                                                                                                        >  2:313866/1‑150 (MQ=255)
                 atGGCTTCTGCAGATCATGAGCTGGATCAAGCCTTTGTACGCTGGGTGTATCTTGAAGCGTTAAAAAATGGTGGCTGATGGCAAAATATTCGATTTTAACTGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTTTCAGGc                                                                                                                                     <  2:285723/150‑1 (MQ=255)
                 atGGCGGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGc                                                                                                                                     <  2:1439549/150‑1 (MQ=255)
                  tGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGCa                                                                                                                                    <  2:1047472/150‑1 (MQ=255)
                   ggCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCAt                                                                                                                                   >  2:1170224/1‑150 (MQ=255)
                    gctgctGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGCaat                                                                                                                                  <  2:1023661/150‑3 (MQ=255)
                    gctgctGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGCaat                                                                                                                                  <  1:575964/150‑3 (MQ=255)
                       gctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATg                                                                                                                                  <  2:944793/147‑1 (MQ=255)
                       gctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATg                                                                                                                                  >  1:944793/1‑147 (MQ=255)
                       gctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGaaa                                                                                                                               >  2:1095724/1‑150 (MQ=255)
                        ctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGaaaa                                                                                                                              >  1:561465/1‑150 (MQ=255)
                         tgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAc                                                                                                                             >  1:799951/1‑150 (MQ=255)
                         tgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAc                                                                                                                             <  2:1523384/150‑1 (MQ=255)
                          gcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGTATGAAAAcc                                                                                                                            <  1:1144584/150‑1 (MQ=255)
                          gcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAcc                                                                                                                            >  2:1228709/1‑150 (MQ=255)
                           cAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCa                                                                                                                           >  2:1418096/1‑150 (MQ=255)
                             gagCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGGCAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAac                                                                                                                         >  1:1132501/1‑150 (MQ=255)
                               gCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAacac                                                                                                                       <  1:896115/150‑1 (MQ=255)
                               gCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAacac                                                                                                                       <  1:1228709/150‑1 (MQ=255)
                                 aGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACtt                                                                                                                     <  1:1161757/150‑1 (MQ=255)
                                   gAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTg                                                                                                                   <  1:87338/150‑1 (MQ=255)
                                    aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAc                                                                                                                             <  1:147870/139‑1 (MQ=255)
                                    aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAc                                                                                                                             >  2:147870/1‑139 (MQ=255)
                                    aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGt                                                                                                                  >  1:219169/1‑150 (MQ=255)
                                     gCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGtt                                                                                                                 >  2:1195904/1‑150 (MQ=255)
                                     gCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGtt                                                                                                                 >  1:502098/1‑150 (MQ=255)
                                       tGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAg                                                                                                               <  2:561465/150‑1 (MQ=255)
                                        ggAGCAAGCCTTTGAACGCTGGGGGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATg                                                                                                                                  <  2:1016849/130‑1 (MQ=255)
                                        ggAGCAAGCCTTTGAACGCTGGGGGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATg                                                                                                                                  >  1:1016849/1‑130 (MQ=255)
                                        ggAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGt                                                                                                              >  1:164452/1‑150 (MQ=255)
                                            cAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACATTGTTAGTCATc                                                                                                          >  2:1409154/1‑150 (MQ=255)
                                             aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCt                                                                                                         >  2:1098424/1‑150 (MQ=255)
                                              aGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTg                                                                                                        >  1:324748/1‑150 (MQ=255)
                                                ccTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTa                                                                                                      >  1:1377757/1‑150 (MQ=255)
                                                  tttGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTaaa                                                                                                    >  1:953364/1‑150 (MQ=255)
                                                   ttGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTaaaa                                                                                                   >  2:1567824/1‑150 (MQ=255)
                                                         gCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAgcgc                                                                                             <  1:814427/150‑1 (MQ=255)
                                                            gggAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGcatc                                                                                          >  1:1489132/1‑150 (MQ=255)
                                                             ggAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGcatca                                                                                         <  2:506131/150‑1 (MQ=255)
                                                                   tCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAgg                                                                                   <  2:109898/150‑1 (MQ=255)
                                                                       gAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTa                                                                               >  2:1609933/1‑150 (MQ=255)
                                                                           cGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCaa                                                                           >  2:108334/1‑150 (MQ=255)
                                                                            gTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAt                                                                          >  2:1565189/1‑150 (MQ=255)
                                                                            gTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAt                                                                          <  1:1095724/150‑1 (MQ=255)
                                                                              tAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTGAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtgt                                                                        >  1:837718/1‑150 (MQ=255)
                                                                              tAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtgt                                                                        >  1:1089466/1‑150 (MQ=255)
                                                                               aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtgtg                                                                       >  1:386509/1‑150 (MQ=255)
                                                                               aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAAGAGTACCAAtgtg                                                                       >  2:1175106/1‑150 (MQ=255)
                                                                                aaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGc                                                                      >  2:817716/1‑150 (MQ=255)
                                                                                  aaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGa                                                                    <  2:1437893/150‑1 (MQ=255)
                                                                                     tggtggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAca                                                                 <  2:549989/150‑1 (MQ=255)
                                                                                        tggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGATTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGGGAAcatca                                                              <  2:604884/150‑1 (MQ=255)
                                                                                          gCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcatcatc                                                            <  2:687366/150‑1 (MQ=255)
                                                                                            tGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcatcatcat                                                          >  1:445279/1‑150 (MQ=255)
                                                                                              aTCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGc                                                        <  1:173719/150‑1 (MQ=255)
                                                                                              aTCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGc                                                        >  1:120710/1‑150 (MQ=255)
                                                                                               tCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGcc                                                       >  1:736160/1‑150 (MQ=255)
                                                                                                 gCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCACCACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCcgc                                                     >  1:1424739/1‑150 (MQ=255)
                                                                                                   aaaaTAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGa                                                   <  2:914104/150‑1 (MQ=255)
                                                                                                     aaTAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAg                                                 <  1:1372157/150‑1 (MQ=255)
                                                                                                      aTAGTCGATTTAAACAGAGACAGCTGACAACTGTGAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGc                                                <  2:837718/150‑1 (MQ=255)
                                                                                                        aGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGcac                                              <  2:858301/150‑1 (MQ=255)
                                                                                                         gTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATATCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGcaca                                             <  2:324748/150‑1 (MQ=255)
                                                                                                         gTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGcaca                                             >  1:1621973/1‑150 (MQ=255)
                                                                                                              tttAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCtg                                        <  2:1462095/150‑1 (MQ=255)
                                                                                                                tAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGTAGCACATCCtgtg                                      >  2:1439646/1‑150 (MQ=255)
                                                                                                                 aaaCAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATAAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGc                                     <  1:46541/150‑1 (MQ=255)
                                                                                                                   aCAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTc                                   >  2:1368318/1‑150 (MQ=255)
                                                                                                                    cAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCg                                  <  2:151831/150‑1 (MQ=255)
                                                                                                                        gaCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATcc                                          >  1:1640324/1‑138 (MQ=255)
                                                                                                                         aCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAgtg                             >  1:1053316/1‑150 (MQ=255)
                                                                                                                           aGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGtgtg                           <  2:324879/150‑1 (MQ=255)
                                                                                                                             cTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGAc                         <  2:1052971/150‑1 (MQ=255)
                                                                                                                                    cGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtgg                  <  2:216442/150‑1 (MQ=255)
                                                                                                                                      gTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtggtg                >  1:952615/1‑150 (MQ=255)
                                                                                                                                        aaaaTTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGc              <  2:458609/150‑1 (MQ=255)
                                                                                                                                        aaaaTTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGc              <  2:1104103/150‑1 (MQ=255)
                                                                                                                                         ttcatGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCg             >  1:1449436/5‑150 (MQ=255)
                                                                                                                                         aaaTTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCg             <  1:108334/150‑1 (MQ=255)
                                                                                                                                         aaaTTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCg             >  2:720724/1‑150 (MQ=255)
                                                                                                                                         aaaTTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCATGTGCTCGCAGTCTGACATGCTGGTGGCg             >  1:1007680/1‑150 (MQ=255)
                                                                                                                                             tGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTAc         >  1:690817/1‑150 (MQ=255)
                                                                                                                                             tGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTAc         >  1:1496126/1‑150 (MQ=255)
                                                                                                                                              gCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTAcc        >  1:503548/1‑150 (MQ=255)
                                                                                                                                                cTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgc      >  1:1640442/1‑150 (MQ=255)
                                                                                                                                                cTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgc      <  1:1170224/150‑1 (MQ=255)
                                                                                                                                                  gatgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCgcgc    >  1:1194527/1‑150 (MQ=255)
                                                                                                                                                    tgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAAGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCt  <  1:1175106/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAAAGTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCT  >  NC_000913/1011777‑1012074

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