Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,011,926 G→A 8.1% 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 8.1% 298.3 / 11.6 98intergenic (+55/‑75)uup/pqiAATP‑binding protein with possible role in replication/intermembrane transport protein PqiA
Reads supporting (aligned to +/- strand):  ref base G (46/44);  new base A (6/2);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 103 total aligned reads displayed.

AAAAAGTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCT  >  NC_000913/1011777‑1012074
                                                                                                                                                     |                                                                                                                                                    
aaaaaGTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCtgatg                                                                                                                                                      <  1:792264/150‑1 (MQ=255)
aaaaaGTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCtgatg                                                                                                                                                      <  1:553767/150‑1 (MQ=255)
       gcttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTggg                                                                                                                                               <  1:12335/150‑1 (MQ=255)
        cttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGt                                                                                                                                              >  1:969137/1‑150 (MQ=255)
         ttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGtt                                                                                                                                             <  2:1318309/150‑1 (MQ=255)
            ctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTAt                                                                                                                                          >  1:8298/1‑150 (MQ=255)
            ctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTAt                                                                                                                                          >  1:13479/1‑150 (MQ=255)
                tatGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAgg                                                                                                                                      <  1:459282/150‑1 (MQ=255)
                  tGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCa                                                                                                                                    <  1:111313/150‑1 (MQ=255)
                  tGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCa                                                                                                                                    <  2:332900/150‑1 (MQ=255)
                   ggCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGTCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCAt                                                                                                                                   >  2:1282241/1‑150 (MQ=255)
                   ggCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGc                                                                                                                                     <  1:998363/148‑1 (MQ=255)
                     ctgctgCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGa                                                                                                                                 >  2:92425/1‑150 (MQ=255)
                        ctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGaaaa                                                                                                                              >  1:550221/1‑150 (MQ=255)
                         tgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGa                                                                                                                                 <  1:932978/146‑1 (MQ=255)
                         tgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGa                                                                                                                                 >  2:932978/1‑146 (MQ=255)
                         tgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAc                                                                                                                             >  2:749210/1‑150 (MQ=255)
                         tgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAc                                                                                                                             >  2:744288/1‑150 (MQ=255)
                          gcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAcc                                                                                                                            >  1:822767/1‑150 (MQ=255)
                           cAGAGCAGGAGCTTGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCa                                                                                                                           >  1:544523/1‑150 (MQ=255)
                             gagCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAac                                                                                                                         >  2:1168811/1‑150 (MQ=255)
                             gagCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAac                                                                                                                         >  2:271192/1‑150 (MQ=255)
                                cAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACt                                                                                                                      >  1:863826/1‑150 (MQ=255)
                                   gAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTg                                                                                                                   >  2:844385/1‑150 (MQ=255)
                                    aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGt                                                                                                                  >  1:1370259/1‑150 (MQ=255)
                                    aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGt                                                                                                                  <  1:471300/150‑1 (MQ=255)
                                    aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGt                                                                                                                  <  1:711783/150‑1 (MQ=255)
                                    aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGt                                                                                                                  >  1:1398445/1‑150 (MQ=255)
                                     gCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGGCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGtt                                                                                                                 <  2:414016/150‑1 (MQ=255)
                                       tGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAg                                                                                                               <  1:844385/150‑1 (MQ=255)
                                       tGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAg                                                                                                               <  2:552384/150‑1 (MQ=255)
                                        ggAGCAAGCCTTTGAACCCTGGGAGTATCTTGAAGCCTTAAAAAATGGTGGCTGATCCCAAAATAGACCATTTAAACAGAGACAGCTGACAAAGGTAAAATTGCCTCAAGAACTGGGTTTATCAGGCATGAAAAACAACACTTTGTTAga                                                                                                              >  2:735302/1‑149 (MQ=255)
                                           gCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCAt                                                                                                           >  2:589700/1‑150 (MQ=255)
                                             aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCt                                                                                                         <  2:550221/150‑1 (MQ=255)
                                             aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCt                                                                                                         >  2:31282/1‑150 (MQ=255)
                                             aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCt                                                                                                         >  1:221895/1‑150 (MQ=255)
                                             aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGc                                                                                                                                     >  1:346377/1‑122 (MQ=255)
                                             aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGc                                                                                                                                     <  2:346377/122‑1 (MQ=255)
                                                 cTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTaa                                                                                                     >  1:1219619/1‑150 (MQ=255)
                                                 cTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTaa                                                                                                     <  1:272196/150‑1 (MQ=255)
                                                   ttGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTaaaa                                                                                                   <  1:1257255/150‑1 (MQ=255)
                                                    tGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAAt                                                                                                  >  1:600353/1‑150 (MQ=255)
                                                     gAACGCTGGGAGTATCTTGACGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATa                                                                                                 >  1:255507/1‑150 (MQ=255)
                                                       aCGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAgc                                                                                               <  1:745354/150‑1 (MQ=255)
                                                        cGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAgcg                                                                                              <  2:391258/150‑1 (MQ=255)
                                                          cTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGca                                                                                            >  2:132118/1‑150 (MQ=255)
                                                             ggAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGcatca                                                                                         >  1:903275/1‑150 (MQ=255)
                                                             ggAGTATCTTGAAGCGTTAAAAAATGGGGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGcatca                                                                                         <  2:608428/150‑1 (MQ=255)
                                                               aGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCAtt                                                                                       <  2:638757/150‑1 (MQ=255)
                                                                  aTCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAg                                                                                    <  2:600353/150‑1 (MQ=255)
                                                                  aTCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAg                                                                                    <  2:832159/150‑1 (MQ=255)
                                                                   tCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAgg                                                                                   <  2:83163/150‑1 (MQ=255)
                                                                    cTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGa                                                                                  >  1:1077386/1‑150 (MQ=255)
                                                                     ttGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAg                                                                                 <  1:31282/150‑1 (MQ=255)
                                                                      tGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGt                                                                                <  2:903275/150‑1 (MQ=255)
                                                                        aaGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTAc                                                                              <  1:1215516/150‑1 (MQ=255)
                                                                               aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtgtg                                                                       <  2:544523/150‑1 (MQ=255)
                                                                               aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtgtg                                                                       >  1:1136915/1‑150 (MQ=255)
                                                                               aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtgtg                                                                       >  1:1033725/1‑150 (MQ=255)
                                                                                  aaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGa                                                                    >  1:1021181/1‑150 (MQ=255)
                                                                                  aaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGa                                                                    <  2:165413/150‑1 (MQ=255)
                                                                                     tggtggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAca                                                                 <  2:255507/150‑1 (MQ=255)
                                                                                      ggtggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcat                                                                >  1:928736/1‑150 (MQ=255)
                                                                                        tggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcatca                                                              <  1:735302/150‑1 (MQ=255)
                                                                                                cGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAgg                                                                                   >  2:605757/1‑121 (MQ=255)
                                                                                                cGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAgg                                                                                   <  1:605757/121‑1 (MQ=255)
                                                                                                cGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCcg                                                      >  1:399557/1‑150 (MQ=255)
                                                                                                  cAAAATAGTCGATTTAAACAGAGTCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCcgcg                                                    <  1:1282241/150‑1 (MQ=255)
                                                                                                     aaTAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAg                                                 <  1:749210/150‑1 (MQ=255)
                                                                                                     aaTAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAg                                                 <  1:744288/150‑1 (MQ=255)
                                                                                                       tAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGca                                               <  1:1380360/150‑1 (MQ=255)
                                                                                                           cGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATATCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACAta                                           >  2:120512/1‑149 (MQ=255)
                                                                                                            gATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATcc                                          <  1:732779/150‑1 (MQ=255)
                                                                                                            gATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATcc                                          >  1:294122/1‑150 (MQ=255)
                                                                                                            gATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATcc                                          <  1:727814/150‑1 (MQ=255)
                                                                                                              tttAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCtg                                        >  1:1277194/1‑150 (MQ=255)
                                                                                                                 aaaCAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGc                                     <  2:284315/150‑1 (MQ=255)
                                                                                                                    cAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAGAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCg                                  >  2:902637/1‑150 (MQ=255)
                                                                                                                       agaCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAg                               >  2:772124/1‑150 (MQ=255)
                                                                                                                        gaCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAgt                              >  2:1217453/1‑150 (MQ=255)
                                                                                                                            gCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGa                          <  1:723892/150‑1 (MQ=255)
                                                                                                                             cTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGAc                         >  2:1026646/1‑150 (MQ=255)
                                                                                                                              tGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACa                        >  2:37222/1‑150 (MQ=255)
                                                                                                                              tGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACa                        <  1:1217453/150‑1 (MQ=255)
                                                                                                                                  aaCGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCt                    <  1:589700/150‑1 (MQ=255)
                                                                                                                                    cGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtgg                  <  2:1033725/150‑1 (MQ=255)
                                                                                                                                     ggTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtggt                 >  1:1251826/1‑150 (MQ=255)
                                                                                                                                      gTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtggtg                <  2:294122/150‑1 (MQ=255)
                                                                                                                                       tAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtggtgg               <  2:822767/150‑1 (MQ=255)
                                                                                                                                            ttGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTa          <  2:1077386/150‑1 (MQ=255)
                                                                                                                                            atGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAATGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTa          >  1:729475/2‑150 (MQ=255)
                                                                                                                                            atGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTa          >  1:574366/2‑150 (MQ=255)
                                                                                                                                            atGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACAACATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTa          >  1:1027537/2‑150 (MQ=255)
                                                                                                                                            atGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGAATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTa          >  2:976241/2‑150 (MQ=255)
                                                                                                                                              gCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTAcc        <  1:687758/150‑1 (MQ=255)
                                                                                                                                              gCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTAcc        >  1:466939/1‑150 (MQ=255)
                                                                                                                                               ccTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcg       >  2:1194174/1‑150 (MQ=255)
                                                                                                                                                cTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgc      >  1:518244/1‑150 (MQ=255)
                                                                                                                                                 tgatgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgcg     <  2:1277194/150‑1 (MQ=255)
                                                                                                                                                 tgatgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgcc     >  2:517581/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAAAGTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCT  >  NC_000913/1011777‑1012074

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