Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,912,748 T→G 7.6% intergenic (‑172/+20) htpX ← / ← prc putative endopeptidase/carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,912,7480TG 7.6% 339.2 / 14.8 115intergenic (‑172/+20)htpX/prcputative endopeptidase/carboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base T (54/52);  new base G (5/4);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 117 total aligned reads displayed.

ACAAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTT  >  NC_000913/1912600‑1912896
                                                                                                                                                    |                                                                                                                                                    
aCAAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAtt                                                                                                                                                     >  2:875742/1‑150 (MQ=255)
  aaaaCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTc                                                                                                                                                   >  2:1153437/1‑150 (MQ=255)
    aaCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTGTATTCTGGTTAGTGATTAAGTCACCAGAAAGACACATCTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCtt                                                                                                                                                 <  1:49649/150‑1 (MQ=255)
       gCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTtgtg                                                                                                                                              <  2:1365764/150‑1 (MQ=255)
         aaTATGCGATACCCACATCGTATGGGTTACTCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACGAGAAAGTCACTACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGcc                                                                                                                                            <  1:82977/150‑1 (MQ=255)
         aaTATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGcc                                                                                                                                            <  2:557374/150‑1 (MQ=255)
         aaTATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTGGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGcc                                                                                                                                            <  2:498383/150‑1 (MQ=255)
         aaTATGCGATACACACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGcc                                                                                                                                            <  1:708836/150‑1 (MQ=255)
          atatGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCt                                                                                                                                           <  1:1227725/150‑1 (MQ=255)
          atatGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCt                                                                                                                                           >  2:270966/1‑150 (MQ=255)
          atatGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCt                                                                                                                                           <  1:775002/150‑1 (MQ=255)
           tatGAGATACCCACATCGGATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCACAAAGACACACCTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCAGTTAAAAAATCAGGCACAAGAAATTGTGCCtg                                                                                                                                          <  2:278884/150‑1 (MQ=255)
             tGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgat                                                                                                                                        <  1:550581/150‑1 (MQ=255)
             tGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgat                                                                                                                                        <  2:1175269/150‑1 (MQ=255)
              gCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCGGATTTACGGCTGCTGTTCGATGTTAAAAAATCAGGCACAATTAATTGTGCCtgatt                                                                                                                                       >  1:1719161/1‑150 (MQ=255)
                 aTACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGACAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCATGTCGCGGGTAAAAAATCAGGCACAATTTCTTGTGCCtgattttt                                                                                                                                    <  2:283488/150‑4 (MQ=255)
                   aCCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCtgatttttta                                                                                                                                  <  1:1117462/150‑6 (MQ=255)
                    cccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATAtt                                                                                                                                 >  1:104749/1‑150 (MQ=255)
                    cccACATCGAATGGGTTACGCGACTATATTCAAGTCTGGATACTGCGTAAGTCACCAGAAAGACACAAGTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCAGTTAAAAAATCAGGCACAATTTCTTGTGCCTGAATGATAtt                                                                                                                                 >  2:1207527/1‑150 (MQ=255)
                     ccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGGCACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTCCCTGATTGATATTa                                                                                                                                >  2:337789/1‑150 (MQ=255)
                     ccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTa                                                                                                                                >  1:927640/1‑150 (MQ=255)
                       acaTCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACt                                                                                                                              <  1:1068102/150‑1 (MQ=255)
                            gTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACg                                                                                                                         <  1:613403/150‑1 (MQ=255)
                              aTGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACATTTTCTTGTGCCTGATTGATATTACTTGACggg                                                                                                                       >  2:1499670/1‑150 (MQ=255)
                              aTGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACggg                                                                                                                       >  1:786852/1‑150 (MQ=255)
                                 ggTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGc                                                                                                                    >  1:1125233/1‑150 (MQ=255)
                                          aCTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTc                                                                                                           >  1:1122923/1‑150 (MQ=255)
                                          aCTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTc                                                                                                           >  1:1117190/1‑150 (MQ=255)
                                          aCTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTc                                                                                                           >  1:60658/1‑150 (MQ=255)
                                          aCTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGAATCTGATTTACGGAATCATGTCGCTGTTAAAAAATCAGGCACAATTTATTGTGCCTGATTGATATTACTTGACGGGAGAGGGTTGTTc                                                                                                           >  1:1547907/1‑150 (MQ=255)
                                            tATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCcg                                                                                                         <  2:1082082/150‑1 (MQ=255)
                                                     gTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTgg                                                                                                >  2:745762/1‑150 (MQ=255)
                                                       cTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTTATTGATATTACTTGACGGGAGAGGGTTGTTCCGCGGGTCTGGCt                                                                                              >  2:331166/1‑150 (MQ=255)
                                                            tAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTc                                                                                         >  1:751106/1‑150 (MQ=255)
                                                             aGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGATCGGGTTGTTCCGCGGGTCTGGCTTTTTCa                                                                                        >  1:852341/1‑150 (MQ=255)
                                                              gTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCaa                                                                                       >  2:522461/1‑150 (MQ=255)
                                                               tGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAg                                                                                      >  1:32392/1‑150 (MQ=255)
                                                                gCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTATCATCTGATTTACTGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGTTTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGc                                                                                     <  1:957628/150‑1 (MQ=255)
                                                                   tAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTc                                                                                  <  2:1267541/150‑1 (MQ=255)
                                                                       tCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCa                                                                              <  1:517435/150‑1 (MQ=255)
                                                                        cACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAg                                                                             >  1:692/1‑150 (MQ=255)
                                                                          ccAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGAt                                                                           <  2:186084/150‑1 (MQ=255)
                                                                          ccAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGAt                                                                           <  2:35951/150‑1 (MQ=255)
                                                                          ccAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGAt                                                                           <  1:1334196/150‑1 (MQ=255)
                                                                             gAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGa                                                                        >  2:705184/1‑150 (MQ=255)
                                                                             gAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGAGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGa                                                                        >  1:413858/1‑150 (MQ=255)
                                                                                aGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGCCGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTg                                                                     >  1:1532155/1‑150 (MQ=255)
                                                                                aGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTg                                                                     <  2:1036943/150‑1 (MQ=255)
                                                                                aGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTg                                                                     <  1:875742/150‑1 (MQ=255)
                                                                                   cacaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                  >  1:970691/1‑150 (MQ=255)
                                                                                   cacaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                  >  1:1368535/1‑150 (MQ=255)
                                                                                    acaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGat                                                                 <  1:268840/150‑1 (MQ=255)
                                                                                     caACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTATTGTGCCTGATAGATATTACTAGACGGGAGAGGGTTGATAATCGGTACTGTCTTTTTAAAGCATCGCCAGATCGAGTGCGata                                                                >  1:33885/1‑150 (MQ=255)
                                                                                           tACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAcc                                                          >  2:1122416/1‑150 (MQ=255)
                                                                                             cATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGt                                                        >  2:821980/1‑150 (MQ=255)
                                                                                               ttttGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtct                                                      <  2:1117190/150‑1 (MQ=255)
                                                                                               ttttGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtct                                                      <  2:1122923/150‑1 (MQ=255)
                                                                                               ttttGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtct                                                      <  2:60658/150‑1 (MQ=255)
                                                                                               ttttGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtct                                                      <  2:1258484/150‑1 (MQ=255)
                                                                                                  tGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCAt                                                   >  1:1442112/1‑150 (MQ=255)
                                                                                                   gTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATc                                                  >  1:468584/1‑150 (MQ=255)
                                                                                                   gTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATc                                                  >  2:30436/1‑150 (MQ=255)
                                                                                                     aGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCa                                                >  2:696871/1‑150 (MQ=255)
                                                                                                      gCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAg                                               >  2:589111/1‑150 (MQ=255)
                                                                                                        aTCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGAt                                             >  1:527658/1‑150 (MQ=255)
                                                                                                          ttGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATaa                                           <  1:181585/149‑1 (MQ=255)
                                                                                                          cTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATaa                                           <  1:1718448/150‑1 (MQ=255)
                                                                                                           tGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAg                                          <  2:413858/150‑1 (MQ=255)
                                                                                                           tGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAg                                          >  1:1635405/1‑150 (MQ=255)
                                                                                                           tGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAg                                          >  1:1205432/1‑150 (MQ=255)
                                                                                                           tGATTTACGGCATCTTGTCGCTGTTAAAAAATCACGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAg                                          <  2:734452/150‑1 (MQ=255)
                                                                                                             aTTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGa                                        <  1:1257669/150‑1 (MQ=255)
                                                                                                                tACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATcc                                     >  2:428480/1‑150 (MQ=255)
                                                                                                                tACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATcc                                     >  2:1189986/1‑150 (MQ=255)
                                                                                                                tACGGCATATTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGCGCGCCGGTTGGGCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCCCATCCAGATAAGGATcc                                     <  2:1138539/150‑1 (MQ=255)
                                                                                                                 tcGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCg                                    <  1:1328687/149‑1 (MQ=255)
                                                                                                                  cGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGTGTTGTTACGCGGGTCTGGCTTTTTCCAGCTTCGCCAGATCGAGTGCGATATTAACCGTCTCATCCAGATAAGGATCCgg                                   >  2:235271/1‑150 (MQ=255)
                                                                                                                  cGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCgg                                   >  2:782236/1‑150 (MQ=255)
                                                                                                                              gCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATc                       <  2:1542503/150‑1 (MQ=255)
                                                                                                                              gCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATc                       <  1:522461/150‑1 (MQ=255)
                                                                                                                               cTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCt                      <  1:782236/150‑1 (MQ=255)
                                                                                                                                 gTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCttt                    >  1:940238/1‑150 (MQ=255)
                                                                                                                                   tAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCg                  <  2:927640/150‑1 (MQ=255)
                                                                                                                                    aaaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCgg                 >  1:1118099/1‑150 (MQ=255)
                                                                                                                                     aaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGt                >  2:110497/1‑150 (MQ=255)
                                                                                                                                     aaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGt                <  1:493641/150‑1 (MQ=255)
                                                                                                                                        aaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGa             <  1:696871/150‑1 (MQ=255)
                                                                                                                                         aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat            >  1:1647975/1‑150 (MQ=255)
                                                                                                                                         aTCAGGCACAAGAAATTGTGCCTGATTGATATTACTTGAAGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCGCCTGGTAATCTTTCGGTTGat            >  1:854551/1‑150 (MQ=255)
                                                                                                                                           cAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatca          <  2:991911/150‑1 (MQ=255)
                                                                                                                                           cAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatca          <  2:786852/150‑1 (MQ=255)
                                                                                                                                           cAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatca          >  1:1635479/1‑150 (MQ=255)
                                                                                                                                           cAGGCACAAGAAATTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatca          >  1:1577157/1‑150 (MQ=255)
                                                                                                                                           cAGGCACAAGAAATTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTATCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATACGGCTCCTGGTAATCTTTCGGTAGatca          >  1:4285/1‑150 (MQ=255)
                                                                                                                                             ggCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatcatc        >  2:864796/1‑150 (MQ=255)
                                                                                                                                             ggCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatcatc        >  2:430452/1‑150 (MQ=255)
                                                                                                                                              gCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATcc       <  2:1293364/150‑1 (MQ=255)
                                                                                                                                              gCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATcc       >  2:493561/1‑150 (MQ=255)
                                                                                                                                                  aaTTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAATCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGtt   <  1:1441818/150‑1 (MQ=255)
                                                                                                                                                   aTTTCTTGTGCCTGATTGAAATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGTTCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGttt  <  2:1547907/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACAAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTT  >  NC_000913/1912600‑1912896

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