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% 320.8 / 7.4 100A104P (GCG→CCG) ‡mrcBfused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (42/52);  new base C (3/3);  total (45/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGA  >  NC_000913/164890‑165187
                                                                                                                                                     |                                                                                                                                                    
gcaaaggGCGTAAGCCTCGTGGCAACCGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAActgcctg                                                                                                                                                      >  2:808642/1‑150 (MQ=255)
 caaaggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAATATTCGTAGCCGTATTGATGGCAAGGTCTGGCAActgcctgc                                                                                                                                                     >  1:306841/1‑150 (MQ=255)
 caaaggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAActgcctgc                                                                                                                                                     >  2:66545/1‑150 (MQ=255)
 caaaggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAActgcctgc                                                                                                                                                     >  2:148726/1‑150 (MQ=255)
   aaggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTATGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCTTAGCCGTATTGATGGCAAGCTCTGGCAACTGCCTGCgg                                                                                                                                                   >  2:1796926/1‑150 (MQ=255)
    aggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGc                                                                                                                                                  <  1:995956/150‑1 (MQ=255)
      ggCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAg                                                                                                                                                <  2:803538/150‑1 (MQ=255)
         gTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGttt                                                                                                                                             >  2:363327/1‑150 (MQ=255)
          tAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGCAGGCAGttgc                                                                                                                                            <  2:461752/150‑3 (MQ=255)
           aaGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGCAGGCAGttgcc                                                                                                                                           <  1:1037924/150‑4 (MQ=255)
           aaGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGCAGGCAGttgcc                                                                                                                                           <  1:434284/150‑4 (MQ=255)
            aGCCTCGTGGCAAACGCGGCTTGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATg                                                                                                                                          <  1:1382488/150‑1 (MQ=255)
            aGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATg                                                                                                                                          <  2:1701723/150‑1 (MQ=255)
            aGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATg                                                                                                                                          <  2:1521056/150‑1 (MQ=255)
             gCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATgg                                                                                                                                         <  2:234515/150‑1 (MQ=255)
             gCCTCGGGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTGTTGCCGTTCTGATCGCCATTTACTGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATgg                                                                                                                                         <  2:234520/150‑1 (MQ=255)
              ccTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGc                                                                                                                                        <  1:1468076/150‑1 (MQ=255)
              ccTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGc                                                                                                                                        <  2:1506460/150‑1 (MQ=255)
              ccTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGc                                                                                                                                        >  2:1446456/1‑150 (MQ=255)
               cTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGGCTGGCAACTGCCTGCGGCAGTTTATGGcc                                                                                                                                       <  1:1669053/150‑1 (MQ=255)
                  gTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGaa                                                                                                                                    <  1:1540669/150‑1 (MQ=255)
                    ggCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATg                                                                                                                                  >  1:692886/1‑150 (MQ=255)
                     gCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATgg                                                                                                                                 >  1:884907/1‑150 (MQ=255)
                       ataCGCGGCTGGCTATTGCTACTTCTAATACTGGCTATCGTTTTTGCCTTGCTGATCGCCATTTTCGGCGTTTATCTCGATCAAAAAATTCGTTTCCGTATTGATGGCAAGGTCTTGCAACTGCCTGCTGCAGTTTATGGCCGAATGGTc                                                                                                                               <  1:927510/148‑1 (MQ=255)
                       aaaCGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTc                                                                                                                               <  1:666389/150‑1 (MQ=255)
                            cgGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCt                                                                                                                          <  2:918118/150‑1 (MQ=255)
                              gctggctATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTg                                                                                                                        <  2:1789822/150‑1 (MQ=255)
                              gctggctATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATATACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTg                                                                                                                        <  2:910636/150‑1 (MQ=255)
                                 ggctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGc                                                                                                                     >  1:1556704/1‑150 (MQ=255)
                                 ggctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGc                                                                                                                     >  1:1713909/1‑150 (MQ=255)
                                  gctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGcc                                                                                                                    <  2:1561780/150‑1 (MQ=255)
                                  gctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGcc                                                                                                                    >  2:741868/1‑150 (MQ=255)
                                   ctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCa                                                                                                                   >  2:721049/1‑150 (MQ=255)
                                        gctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACAt                                                                                                              <  2:1537930/150‑1 (MQ=255)
                                           aCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGAc                                                                                                           <  1:66545/150‑1 (MQ=255)
                                             tGCTATAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCTTAGCCGTTTTGATGGCAAGGTCTGGCAACTGCATGCTTCAGTTTATGGACGAATGGTCTATATTGAGAAAGACATGAcaa                                                                                                         >  1:249920/1‑148 (MQ=255)
                                              gCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCAt                                                                                                        >  1:629096/1‑150 (MQ=255)
                                              gCTAAAACTGGCTATCGTTTTTGCCGTGATGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCAt                                                                                                        >  1:1214503/1‑150 (MQ=255)
                                                tAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCa                                                                                                      <  1:343455/150‑1 (MQ=255)
                                                 aaaaCTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAg                                                                                                     <  1:989948/150‑1 (MQ=255)
                                                   aaCTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCa                                                                                                   >  1:107363/1‑150 (MQ=255)
                                                       ggCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGaa                                                                                               >  1:195317/1‑150 (MQ=255)
                                                         cTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACg                                                                                             >  1:965851/1‑150 (MQ=255)
                                                               tttttGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATgg                                                                                       >  1:608120/1‑150 (MQ=255)
                                                               tgtttGCCGTGCTGATCGCCATTTAGGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATgg                                                                                       <  1:335159/148‑1 (MQ=255)
                                                                  ttGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAATATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCATTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGa                                                                                    <  2:306841/150‑1 (MQ=255)
                                                                    gCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCAGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAg                                                                                  >  1:1382352/1‑150 (MQ=255)
                                                                    gCCGTGCGGATCGCCATTTACGGCGTTTCTCTCGATCAAAAAATTCGTAGCCGTATTGATGGCTAGGGCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAg                                                                                  <  2:1005175/150‑1 (MQ=255)
                                                                      cGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgct                                                                                <  2:1243266/150‑1 (MQ=255)
                                                                        tGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgctgc                                                                              >  1:1334384/1‑150 (MQ=255)
                                                                         gCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgctgct                                                                             <  1:1484448/150‑1 (MQ=255)
                                                                          cTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGctgctg                                                                            >  2:662500/1‑150 (MQ=255)
                                                                            gATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGa                                                                          >  2:1299306/1‑150 (MQ=255)
                                                                            gATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGa                                                                          >  1:1641177/1‑150 (MQ=255)
                                                                              tCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAAATGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGACAGACATGACCATCAGCAAGAACTAGATGGTGAAGCTGCTGGAgg                                                                        >  2:1167834/1‑150 (MQ=255)
                                                                               cGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGc                                                                       <  1:1418634/150‑1 (MQ=255)
                                                                                gCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCg                                                                      >  2:150982/1‑150 (MQ=255)
                                                                                  cATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGAAATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGAc                                                                    >  1:128820/1‑150 (MQ=255)
                                                                                         ggCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTAt                                                             <  2:1163365/150‑1 (MQ=255)
                                                                                         ggCGGTTATCTCTATCAAAAAATTCGTAGCCGTAAAGATCGCAAGGTCTGGCAACGGCCTGCGGCAGTTTATGGCCGAATGGGCCAGCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCGGGAGGCGACCCAGTAt                                                             <  2:296516/150‑1 (MQ=255)
                                                                                                 tctcGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGgt                                                     >  1:1396074/1‑150 (MQ=255)
                                                                                                   tcGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGgtgt                                                   <  1:1336103/150‑1 (MQ=255)
                                                                                                     gATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCg                                                 <  2:1382352/150‑1 (MQ=255)
                                                                                                      aTCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGa                                                <  2:1282804/150‑1 (MQ=255)
                                                                                                       tCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGACGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGaa                                               >  2:844822/1‑150 (MQ=255)
                                                                                                        cAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGaaa                                              <  1:184907/150‑1 (MQ=255)
                                                                                                        cAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGaaa                                              >  1:251965/1‑150 (MQ=255)
                                                                                                        cAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAAAGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTAAGGTGTCGaaa                                              >  1:1768056/1‑150 (MQ=255)
                                                                                                        aaaaaaaTTCGTAGCCGTATTGATGGCAAGGTCTGTCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCATGTACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGaaa                                              <  1:808642/149‑1 (MQ=255)
                                                                                                         aaaaaaTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGaaaa                                             <  1:148726/150‑1 (MQ=255)
                                                                                                              aTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGAcc                                        <  1:729111/150‑1 (MQ=255)
                                                                                                                   tAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTcc                                   >  2:528153/1‑150 (MQ=255)
                                                                                                                    aGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCt                                  <  2:1087159/150‑1 (MQ=255)
                                                                                                                    aGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCt                                  >  1:185032/1‑150 (MQ=255)
                                                                                                                     gCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTg                                 >  1:862272/1‑150 (MQ=255)
                                                                                                                     gCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTg                                 >  1:852943/1‑150 (MQ=255)
                                                                                                                      ccGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTgg                                <  1:781082/150‑1 (MQ=255)
                                                                                                                      ccGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTgg                                <  1:736761/150‑1 (MQ=255)
                                                                                                                       cGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGc                               <  1:1757705/150‑1 (MQ=255)
                                                                                                                       cGTATTGATGGCAAGGGCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGc                               <  1:1515164/150‑1 (MQ=255)
                                                                                                                        gTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCg                              >  1:793497/1‑150 (MQ=255)
                                                                                                                            tGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAtt                          <  2:8100/150‑1 (MQ=255)
                                                                                                                              aTGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTa                        >  2:1462048/1‑150 (MQ=255)
                                                                                                                               tGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTAc                       <  1:363327/150‑1 (MQ=255)
                                                                                                                                ggCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTAcc                      <  1:1031434/150‑1 (MQ=255)
                                                                                                                                 gCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCg                     <  2:1821025/150‑1 (MQ=255)
                                                                                                                                   aaGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTg                   <  2:416709/150‑1 (MQ=255)
                                                                                                                                     ggTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGCCCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCa                 >  2:1100268/1‑150 (MQ=255)
                                                                                                                                     ggTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCa                 >  1:1191603/1‑150 (MQ=255)
                                                                                                                                        cTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGAGTCGAAAATGACCCGTACTGGCGAATTTACCGTGCAGGc              >  1:85486/1‑150 (MQ=255)
                                                                                                                                         cataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAGTTACCGTGCAGGcc             >  1:308313/5‑150 (MQ=255)
                                                                                                                                         cataaaCTGCCGCAGGCAGGTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAGTGTACCGTGCCGGcc             >  1:141761/5‑150 (MQ=255)
                                                                                                                                          ggCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCa            >  1:1631142/1‑150 (MQ=255)
                                                                                                                                            cAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGACCGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAAc          <  2:134445/150‑1 (MQ=255)
                                                                                                                                            cAACTGCCTGCGGCAGTTTAAGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTACTGGCGAATTTACCGTGCAGGCCaaa          >  1:1552857/1‑149 (MQ=255)
                                                                                                                                             aaCTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACa         <  2:884907/150‑1 (MQ=255)
                                                                                                                                             aaCTGCCGCAGGCAGTTTATGGACGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATTACCCGTCCTGGCGAATTTACCGTGCAGGCCAACa         >  2:1302351/1‑150 (MQ=255)
                                                                                                                                                 gcctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCAt     >  1:343677/1‑150 (MQ=255)
                                                                                                                                                   gtgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGTGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTg   <  1:1699755/149‑1 (MQ=255)
                                                                                                                                                    tgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTgt  >  1:666009/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGA  >  NC_000913/164890‑165187

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: