Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 165,039 G→C 5.8% A104P (GCG→CCG) ‡ mrcB → fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913165,0390GC 5.8% 318.0 / 5.0 96A104P (GCG→CCG) ‡mrcBfused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (45/45);  new base C (4/2);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

CAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAG  >  NC_000913/164891‑165188
                                                                                                                                                    |                                                                                                                                                     
caaaggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAActgcctgc                                                                                                                                                      >  2:829032/1‑150 (MQ=255)
 aaaggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCg                                                                                                                                                     <  1:868745/150‑1 (MQ=255)
   aggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTTCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGc                                                                                                                                                   <  2:1696331/150‑1 (MQ=255)
   aggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGc                                                                                                                                                   >  1:1395426/1‑150 (MQ=255)
    gggCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCa                                                                                                                                                  >  2:574122/1‑150 (MQ=255)
     ggCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAg                                                                                                                                                 <  1:953826/150‑1 (MQ=255)
      gCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGt                                                                                                                                                >  2:78525/1‑150 (MQ=255)
      gCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGt                                                                                                                                                >  1:245836/1‑150 (MQ=255)
         tAAGCCTCGTGGCAAACGCGGCTGGCTATTGCTACTGCTAAAACTTGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTa                                                                                                                                             <  1:873269/150‑1 (MQ=255)
         tAAGCCGCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGGGCTGATCGCCATTTACGGCGGTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGCAGGCAGttgc                                                                                                                                             <  1:311959/150‑3 (MQ=255)
          aaGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGCAGGCAGttgcc                                                                                                                                            <  2:1116394/150‑4 (MQ=255)
            gCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGTCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGAAAGGTCTGGCAACTGCCTGCGGCAGTTTATgg                                                                                                                                          >  2:675519/1‑150 (MQ=255)
            gCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATgg                                                                                                                                          <  2:1205378/150‑1 (MQ=255)
            gCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATgg                                                                                                                                          >  1:17859/1‑150 (MQ=255)
              cTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGcc                                                                                                                                        <  1:517734/150‑1 (MQ=255)
               gcGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCg                                                                                                                                       <  1:1248282/149‑1 (MQ=255)
                   ggCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATg                                                                                                                                   <  2:38179/150‑1 (MQ=255)
                    gCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATgg                                                                                                                                  <  2:1639099/150‑1 (MQ=255)
                    gCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGACATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATgg                                                                                                                                  >  2:1399823/1‑150 (MQ=255)
                      aaaCGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGACCGAATGGTc                                                                                                                                <  2:245836/150‑1 (MQ=255)
                         cgctGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAAt                                                                                                                             <  2:36694/150‑1 (MQ=255)
                         cgcgGCTGGCTATGTCTACTGATAAACCTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGGCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAAt                                                                                                                             <  2:97722/150‑1 (MQ=255)
                         cgcgGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAAt                                                                                                                             <  2:42130/150‑1 (MQ=255)
                         cgcgGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAAt                                                                                                                             >  1:195133/1‑150 (MQ=255)
                         cgcgGCTGGCTATGGCTAATGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAAt                                                                                                                             <  2:1211302/150‑1 (MQ=255)
                             gctggctATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTg                                                                                                                         <  1:1624654/150‑1 (MQ=255)
                                ggctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCACCTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGc                                                                                                                      >  2:1701839/1‑150 (MQ=255)
                                gactatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGc                                                                                                                      >  2:1311726/3‑150 (MQ=255)
                                 gctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAgac                                                                                                                     >  1:299958/1‑148 (MQ=255)
                                 gctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGcc                                                                                                                     >  1:745725/1‑150 (MQ=255)
                                   tatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAg                                                                                                                   >  2:46029/1‑150 (MQ=255)
                                   tatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAg                                                                                                                   >  2:1124689/1‑150 (MQ=255)
                                     tggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGAc                                                                                                                 >  2:968483/1‑150 (MQ=255)
                                        ctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGATGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATg                                                                                                              <  2:1279723/150‑1 (MQ=255)
                                           cTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGAcc                                                                                                           <  1:1568031/150‑1 (MQ=255)
                                             gCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGACATTTACGGCGTTTATCTCGATCAAAAAATTCGAAGCAGTATTTATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCaa                                                                                                         >  2:1310416/1‑149 (MQ=255)
                                              cTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATc                                                                                                        <  1:119513/150‑1 (MQ=255)
                                               tAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCa                                                                                                       <  2:1380771/150‑1 (MQ=255)
                                                aaaaCTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAg                                                                                                      <  1:1281792/150‑1 (MQ=255)
                                                 aaaCTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGc                                                                                                     >  2:844256/1‑150 (MQ=255)
                                                    cTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGc                                                                                                     >  2:142080/1‑147 (MQ=255)
                                                    cTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGc                                                                                                     <  1:142080/147‑1 (MQ=255)
                                                      ggCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGaa                                                                                                >  2:698446/1‑150 (MQ=255)
                                                       gCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAAc                                                                                               <  1:1201141/150‑1 (MQ=255)
                                                        cTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTTAGCCAGACATGACCATCAGCAAGAACg                                                                                              >  2:1428205/1‑150 (MQ=255)
                                                        cTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACg                                                                                              >  2:1162622/1‑150 (MQ=255)
                                                         tttCGTTTTTGCCGTGCTGATCGCCATTTTCGGAGTTTATCTCGTTCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCATCTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACga                                                                                             <  2:871190/148‑1 (MQ=255)
                                                            cGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGAt                                                                                          >  2:817915/1‑150 (MQ=255)
                                                            cGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGAt                                                                                          <  2:1480824/150‑1 (MQ=255)
                                                            cGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGAt                                                                                          >  2:701864/1‑150 (MQ=255)
                                                             gTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATg                                                                                         >  1:202393/1‑150 (MQ=255)
                                                              tttttGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATgg                                                                                        <  1:78525/150‑1 (MQ=255)
                                                                   gCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAg                                                                                   >  2:799998/1‑150 (MQ=255)
                                                                   gCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAg                                                                                   >  2:55265/1‑150 (MQ=255)
                                                                    ccGTGCTGATCGGCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATTGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgc                                                                                  <  1:392357/150‑1 (MQ=255)
                                                                    ccGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgc                                                                                  >  1:1028445/1‑150 (MQ=255)
                                                                    ccGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgc                                                                                  >  1:481442/1‑150 (MQ=255)
                                                                     cGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgct                                                                                 >  2:1177415/1‑150 (MQ=255)
                                                                     cGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgct                                                                                 >  1:1529434/1‑150 (MQ=255)
                                                                          ttaTCGCCATTTACGGCGTTTATCTCGATCAAAATATTCGTAGCCTTATTGATGGCAAGGTCGGGCAACTGCCGGCGGCAGTTTATGGCCGAATGTTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTgg                                                                            <  1:1125898/148‑1 (MQ=255)
                                                                          tGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTgg                                                                            <  2:1166351/150‑1 (MQ=255)
                                                                           gATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGa                                                                           >  2:817153/1‑150 (MQ=255)
                                                                           gATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGa                                                                           <  1:829032/150‑1 (MQ=255)
                                                                             tCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAgg                                                                         >  1:1710847/1‑150 (MQ=255)
                                                                             tCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAgg                                                                         <  1:1224294/150‑1 (MQ=255)
                                                                              cGCCATTAACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGc                                                                        <  2:171672/150‑1 (MQ=255)
                                                                               gCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCg                                                                       >  1:1614011/1‑150 (MQ=255)
                                                                                 cATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGAc                                                                     >  1:927740/1‑150 (MQ=255)
                                                                                   tttACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGAccc                                                                   <  1:1176901/150‑1 (MQ=255)
                                                                                     tACGGCGTTTATCTCGATCAAAAAAATCGTAGCCGTATTGACGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAg                                                                 <  1:316030/150‑1 (MQ=255)
                                                                                             ttATCTCGATCAAAAAAATCGTACCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCATTTTATGGCCGAATGGTCAATATTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCa                                                         <  1:675519/150‑1 (MQ=255)
                                                                                               aTCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGa                                                                           <  2:1706256/130‑1 (MQ=255)
                                                                                               aTCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGa                                                                           >  1:1706256/1‑130 (MQ=255)
                                                                                                          aaaaTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATg                                            >  1:1125639/1‑150 (MQ=255)
                                                                                                                        tATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGa                              >  2:32179/1‑150 (MQ=255)
                                                                                                                         aTTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGaa                             >  2:667609/1‑150 (MQ=255)
                                                                                                                          ttGATGGCAAGGTCTGGCAACTGCCTGAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAt                            <  1:1279761/150‑1 (MQ=255)
                                                                                                                           tGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGAATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAtt                           >  1:62272/1‑150 (MQ=255)
                                                                                                                            gATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAttt                          <  1:122571/150‑1 (MQ=255)
                                                                                                                                   aGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGc                   <  1:1311726/150‑1 (MQ=255)
                                                                                                                                    ggTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTAGCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCa                  <  1:46029/150‑1 (MQ=255)
                                                                                                                                      tCTGGCAACTGCCTGCGGCATTTTATGGCCTAATGGTCAATCTTGAGCCTGAGATGACCATCAGCAAGAACGAGATGGTGATGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGATAATGACCCGTCCTGGCGAATTTACCGTGCAgg                <  1:1310416/150‑1 (MQ=255)
                                                                                                                                      tCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAgg                <  1:1124689/150‑1 (MQ=255)
                                                                                                                                      tCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTCGCGAATTTACCGTGCAgg                >  2:829850/1‑150 (MQ=255)
                                                                                                                                      gccataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAgg                >  2:3042/7‑150 (MQ=255)
                                                                                                                                       cTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGCCCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGc               >  1:1214302/1‑150 (MQ=255)
                                                                                                                                       cTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGc               <  1:983459/150‑1 (MQ=255)
                                                                                                                                        gggCAACTGCCTGCGGCAGTTTATGGCCGAATGGCCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGcc              <  2:271018/149‑1 (MQ=255)
                                                                                                                                        cataaaCTGCCGCAGGCAGTTTATGGCCGATTGGTCAATCTTGCGCCAGACATGACCATCCGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGcc              >  1:1253088/5‑150 (MQ=255)
                                                                                                                                        cataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGCCCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGcc              >  2:698299/5‑150 (MQ=255)
                                                                                                                                         ggCCACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCa             <  1:817915/150‑1 (MQ=255)
                                                                                                                                         ggCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCa             >  2:583524/1‑150 (MQ=255)
                                                                                                                                         ggCAACTGCCAGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGTCCa             >  2:36768/1‑150 (MQ=255)
                                                                                                                                         ataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCa             >  1:144679/4‑150 (MQ=255)
                                                                                                                                          gCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCaa            >  2:842110/1‑150 (MQ=255)
                                                                                                                                             aCTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAg         <  1:842110/150‑1 (MQ=255)
                                                                                                                                                gcctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCAt      >  1:1172017/1‑150 (MQ=255)
                                                                                                                                                 cctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCAtt     >  1:914736/1‑150 (MQ=255)
                                                                                                                                                    gcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTgag  >  1:1293688/1‑150 (MQ=255)
                                                                                                                                                    gcGGCAGTTTATGGCCGAATGGTAAATCTTGAGCCCGACATGACCATCAGCAAGAACGCGATGGTGAAGCTGCTGGAGGCGACCCAGTATCTTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTgag  <  2:822964/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAG  >  NC_000913/164891‑165188

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