Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,011,930 T→C 7.5% intergenic (+59/‑71) 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,9300TC 7.5% 316.6 / 12.1 106intergenic (+59/‑71)uup/pqiAATP‑binding protein with possible role in replication/intermembrane transport protein PqiA
Reads supporting (aligned to +/- strand):  ref base T (50/48);  new base C (6/2);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAG  >  NC_000913/1011784‑1012078
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                                        <  2:332900/150‑1 (MQ=255)
           tGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCa                                                                                                                                        <  1:111313/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                                                                                                                                     >  2:932978/1‑146 (MQ=255)
                  tgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAc                                                                                                                                 >  2:744288/1‑150 (MQ=255)
                  tgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAc                                                                                                                                 >  2:749210/1‑150 (MQ=255)
                   gcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAcc                                                                                                                                >  1:822767/1‑150 (MQ=255)
                    cAGAGCAGGAGCTTGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCa                                                                                                                               >  1:544523/1‑150 (MQ=255)
                      gagCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAac                                                                                                                             >  2:271192/1‑150 (MQ=255)
                      gagCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAac                                                                                                                             >  2:1168811/1‑150 (MQ=255)
                         cAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACt                                                                                                                          >  1:863826/1‑150 (MQ=255)
                            gAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTg                                                                                                                       >  2:844385/1‑150 (MQ=255)
                             aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGt                                                                                                                      <  1:711783/150‑1 (MQ=255)
                             aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGt                                                                                                                      >  1:1370259/1‑150 (MQ=255)
                             aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGt                                                                                                                      <  1:471300/150‑1 (MQ=255)
                              gCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGGCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGtt                                                                                                                     <  2:414016/150‑1 (MQ=255)
                                tGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAg                                                                                                                   <  2:552384/150‑1 (MQ=255)
                                tGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAg                                                                                                                   <  1:844385/150‑1 (MQ=255)
                                 ggAGCAAGCCTTTGAACCCTGGGAGTATCTTGAAGCCTTAAAAAATGGTGGCTGATCCCAAAATAGACCATTTAAACAGAGACAGCTGACAAAGGTAAAATTGCCTCAAGAACTGGGTTTATCAGGCATGAAAAACAACACTTTGTTAga                                                                                                                  >  2:735302/1‑149 (MQ=255)
                                    gCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCAt                                                                                                               >  2:589700/1‑150 (MQ=255)
                                      aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCt                                                                                                             >  1:221895/1‑150 (MQ=255)
                                      aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCt                                                                                                             >  2:31282/1‑150 (MQ=255)
                                      aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCt                                                                                                             <  2:550221/150‑1 (MQ=255)
                                      aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGc                                                                                                                                         <  2:346377/122‑1 (MQ=255)
                                      aaGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGc                                                                                                                                         >  1:346377/1‑122 (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:832159/150‑1 (MQ=255)
                                                           aTCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAg                                                                                        <  2:600353/150‑1 (MQ=255)
                                                            tCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAgg                                                                                       >  2:1398167/1‑150 (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                                                                           >  1:1033725/1‑150 (MQ=255)
                                                                        aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtgtg                                                                           <  2:544523/150‑1 (MQ=255)
                                                                        aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtgtg                                                                           >  1:1136915/1‑150 (MQ=255)
                                                                           aaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGa                                                                        <  2:165413/150‑1 (MQ=255)
                                                                           aaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGa                                                                        >  1:1021181/1‑150 (MQ=255)
                                                                              tggtggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAca                                                                     <  2:255507/150‑1 (MQ=255)
                                                                               ggtggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcat                                                                    >  1:928736/1‑150 (MQ=255)
                                                                                 tggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcatca                                                                  <  1:735302/150‑1 (MQ=255)
                                                                                         cGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCcg                                                          >  1:399557/1‑150 (MQ=255)
                                                                                           cAAAATAGTCGATTTAAACAGAGTCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCcgcg                                                        <  1:1282241/150‑1 (MQ=255)
                                                                                              aaTAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAg                                                     <  1:744288/150‑1 (MQ=255)
                                                                                              aaTAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAg                                                     <  1:749210/150‑1 (MQ=255)
                                                                                                tAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGca                                                   <  1:1380360/150‑1 (MQ=255)
                                                                                                    cGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATATCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACAta                                               >  2:120512/1‑149 (MQ=255)
                                                                                                     gATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATcc                                              >  1:294122/1‑150 (MQ=255)
                                                                                                     gATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATcc                                              <  1:727814/150‑1 (MQ=255)
                                                                                                     gATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATcc                                              <  1:732779/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)
                                                                                                                                             tgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCt      >  1:594918/1‑150 (MQ=255)
                                                                                                                                               aaCTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTg    <  2:13479/150‑1 (MQ=255)
                                                                                                                                               aaCTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTg    <  2:8298/150‑1 (MQ=255)
                                                                                                                                               aaCTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTg    >  2:1085297/1‑150 (MQ=255)
                                                                                                                                                 cTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAg  >  2:253386/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAG  >  NC_000913/1011784‑1012078

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: