Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 165,038 T→G 5.1% P103P (CCT→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,0380TG 5.1% 437.4 / 10.9 150P103P (CCT→CCGmrcBfused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (74/68);  new base G (6/2);  total (80/70)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 160 total aligned reads displayed.

GGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGA  >  NC_000913/164889‑165187
                                                                                                                                                     |                                                                                                                                                     
    aaggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCgg                                                                                                                                                   >  1:862064/1‑150 (MQ=255)
     aggGCGTAAGGCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTAGCGTTTTTGCCGTGCTGATCGCGATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGc                                                                                                                                                  <  1:46807/150‑1 (MQ=255)
     aggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTTGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGc                                                                                                                                                  >  2:804341/1‑150 (MQ=255)
     aggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCTTATCCGTATTGATGGCAAGGTATGGCAACTGCCTGCGGc                                                                                                                                                  >  2:166516/1‑150 (MQ=255)
        gCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGt                                                                                                                                               >  1:781163/1‑150 (MQ=255)
         cGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGtt                                                                                                                                              <  1:1056212/150‑1 (MQ=255)
         cGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGtt                                                                                                                                              <  2:553652/150‑1 (MQ=255)
            aaGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGGCTGGCAACTGCCGCAGGCAGttgcc                                                                                                                                           <  1:2014498/150‑4 (MQ=255)
             aGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATg                                                                                                                                          >  2:841622/1‑150 (MQ=255)
              gCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATgg                                                                                                                                         >  1:1395533/1‑150 (MQ=255)
               ccTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGc                                                                                                                                        >  1:766286/1‑150 (MQ=255)
                cTCGTGGCAAACGCGGCTGGCAATGGCTACTGCTAAAACTGGCAATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGcc                                                                                                                                       <  1:1137777/150‑1 (MQ=255)
                 tCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCg                                                                                                                                      <  1:1004082/150‑1 (MQ=255)
                        aaaCGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTc                                                                                                                               <  2:1895901/150‑1 (MQ=255)
                            gcgGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGACGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATc                                                                                                                           >  2:373575/1‑150 (MQ=255)
                               gctggctATGGCTACTGCTAAAACTGGCTATCGTTTTTGCGGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTg                                                                                                                        <  2:1261203/150‑1 (MQ=255)
                                 tggctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAg                                                                                                                      <  1:2014618/150‑1 (MQ=255)
                                   gctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGcc                                                                                                                    <  1:1638112/150‑1 (MQ=255)
                                        ggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACa                                                                                                               >  2:1031446/1‑150 (MQ=255)
                                         gctaCTGCTAAAACTGGCTATCGTTTTAGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACAt                                                                                                              >  2:725907/1‑150 (MQ=255)
                                           taCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGa                                                                                                            >  2:2144222/1‑150 (MQ=255)
                                            aCTGCTAAAACTGGCTATCGTTTTTTCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGAc                                                                                                           <  1:1992403/150‑1 (MQ=255)
                                            aCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGAc                                                                                                           >  2:762078/1‑150 (MQ=255)
                                             cTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATACGTAGCCGTATTGCTGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGGCAATCTTGAGCCAGACATGAcc                                                                                                          <  1:1218314/150‑1 (MQ=255)
                                              tGCTAAAACTTGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCa                                                                                                         <  1:1817868/150‑1 (MQ=255)
                                              tGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCa                                                                                                         <  1:1252956/150‑1 (MQ=255)
                                               gCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGCCATGACCAt                                                                                                        >  2:1854135/1‑150 (MQ=255)
                                               gCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTGTATCTCGATCAAAATATTCGTAGCCGTATTGATGGCACGGTCTGGCAACTGCCTGCGGCAGTTTATGGACGAATGGTCAATCTTGAGACAGACATGACCAt                                                                                                        >  2:258000/1‑150 (MQ=255)
                                                 tAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTTCGGCAGTTTATGGCCGAATGGTCAATCTTGGGCCAGACATGAACATCa                                                                                                      >  1:450963/1‑150 (MQ=255)
                                                 tAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCa                                                                                                      >  2:639075/1‑150 (MQ=255)
                                                 tAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCa                                                                                                      >  1:698958/1‑150 (MQ=255)
                                                  aaaaCTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAg                                                                                                     >  2:1526465/1‑150 (MQ=255)
                                                    aaCTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCa                                                                                                   >  2:2035471/1‑150 (MQ=255)
                                                        ggCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGaa                                                                                               >  2:959186/1‑150 (MQ=255)
                                                        ggCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGaa                                                                                               >  1:423618/1‑150 (MQ=255)
                                                         gCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAAc                                                                                              >  1:1995710/1‑150 (MQ=255)
                                                          cTATCGTGATTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGGCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACg                                                                                             <  2:42418/150‑1 (MQ=255)
                                                            aTCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACgag                                                                                           <  1:1266502/150‑1 (MQ=255)
                                                            aTCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACgag                                                                                           <  2:1012274/150‑1 (MQ=255)
                                                               gTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATg                                                                                        >  1:972217/1‑150 (MQ=255)
                                                               gTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGACTGGCAACTGCCTGCGGCAGTTTATCGCCGAATGTTCAATCTTGAGCCAGACATGACCATCAGCAAGACAGAGATg                                                                                        >  2:159366/1‑150 (MQ=255)
                                                                tttttGCCGTGATGATAGCCATTTACGGCGTTTATCTCGATCAAACTATTCGTAGCAGTATTGATGGCAAGGTCTGGCAACTGCCGGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATgg                                                                                       <  1:407205/150‑1 (MQ=255)
                                                                    tGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGaa                                                                                   <  1:662373/150‑1 (MQ=255)
                                                                    tGCCGTGCTGATCGACATTTACGTCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCTAGAACGAGATGGTGaa                                                                                   <  1:2110559/150‑1 (MQ=255)
                                                                      ccGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgc                                                                                 <  1:592328/150‑1 (MQ=255)
                                                                      ccGTGCTGATCGCCATTTACGGCGTTTATCTCGATAAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgc                                                                                 <  1:1390821/150‑1 (MQ=255)
                                                                       cGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGGCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgct                                                                                >  2:1210669/1‑150 (MQ=255)
                                                                          gCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgctgct                                                                             >  1:714664/1‑150 (MQ=255)
                                                                             gATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGAAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGATATGGTGAAgctgacgga                                                                          >  1:1860069/1‑145 (MQ=255)
                                                                                cGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGc                                                                       <  1:639075/150‑1 (MQ=255)
                                                                                cGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGc                                                                       <  1:1680746/150‑1 (MQ=255)
                                                                                 gCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCg                                                                      >  1:1633205/1‑150 (MQ=255)
                                                                                  ccATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCgc                                                                     >  1:2159242/1‑149 (MQ=255)
                                                                                     tttACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGAccc                                                                  >  1:2046487/1‑150 (MQ=255)
                                                                                           gCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATc                                                            >  2:1816637/1‑150 (MQ=255)
                                                                                                 aTCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTCGCAACAGCCTGCGGCAGTTTATGTCCGAATGGTCAATCTAGAGACAGACATGACCATCAGCAAGATCGAGATGGTGATGCTGCTGGAGGCGACCCAGTCTCTTCAgg                                                      >  2:236062/1‑150 (MQ=255)
                                                                                                 aTCTCGATCAAAAAATTAGTAGCCGTATTGATGGCAAGGTGTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAgg                                                      <  2:1170390/150‑1 (MQ=255)
                                                                                                 aTCTCCATCAAAAAATTCGCCGCCCAATTCATCGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAgg                                                      <  2:2192553/150‑1 (MQ=255)
                                                                                                   ctcGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGgtg                                                    >  1:600142/1‑150 (MQ=255)
                                                                                                   ctcGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGgtg                                                    <  2:781163/150‑1 (MQ=255)
                                                                                                     cGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATTACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGACGCGACCCACTATCGTCAGGTGTc                                                  >  2:684112/1‑150 (MQ=255)
                                                                                                      gATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCg                                                 >  1:1803482/1‑150 (MQ=255)
                                                                                                       aTCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGa                                                <  2:862064/150‑1 (MQ=255)
                                                                                                         cAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGaaa                                              >  1:2210134/1‑150 (MQ=255)
                                                                                                         cAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGAGACCCAGTATCGTCAGGTGTCTaaa                                              >  1:1492852/1‑150 (MQ=255)
                                                                                                           aaaaaTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGCCCATCAGCAAGCACGAGATGGTGCAGCTGCTGGAGGCGCCCCAGTATCGTCCGGTGTCGAAAAt                                            >  1:305438/1‑150 (MQ=255)
                                                                                                            aaacTTCGTAGCCGTATTGAGGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATg                                           <  2:1722133/150‑1 (MQ=255)
                                                                                                               aTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGAcc                                        <  2:714664/150‑1 (MQ=255)
                                                                                                                   gTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTc                                    >  1:294238/1‑150 (MQ=255)
                                                                                                                         gTATTGATGGCAAGGGCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCg                              <  1:1658592/150‑1 (MQ=255)
                                                                                                                           aTTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAACCTGATGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGaa                            >  2:1993111/1‑150 (MQ=255)
                                                                                                                            ttGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAt                           >  2:2141067/1‑150 (MQ=255)
                                                                                                                              gATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGTAttt                         >  1:670688/1‑150 (MQ=255)
                                                                                                                              gATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAttt                         >  1:868545/1‑150 (MQ=255)
                                                                                                                              gATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAttt                         >  1:1221888/1‑150 (MQ=255)
                                                                                                                               aTGGCAAGGGCAGGCAACTGCCTGCGGGAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGACCGAGATGGGGTAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTa                        <  1:1403950/150‑1 (MQ=255)
                                                                                                                                tGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTAc                       >  2:2074008/1‑150 (MQ=255)
                                                                                                                                 ggCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTAcc                      >  1:1649007/1‑150 (MQ=255)
                                                                                                                                  gCAAGGTCTGGCAACTGCCTGCGGCATTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCg                     <  1:1979384/150‑1 (MQ=255)
                                                                                                                                  gCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCg                     >  1:1915420/1‑150 (MQ=255)
                                                                                                                                     aGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGAGCAGCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGc                  <  1:373575/150‑1 (MQ=255)
                                                                                                                                     aGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGc                  <  1:1276451/150‑1 (MQ=255)
                                                                                                                                     aGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGc                  >  1:2192246/1‑150 (MQ=255)
                                                                                                                                      ggTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCa                 >  2:175021/1‑150 (MQ=255)
                                                                                                                                        tcgggCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAgg               <  2:379868/147‑1 (MQ=255)
                                                                                                                                        tCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCATCAAGAACGAGATGGGGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAgg               <  1:69155/150‑1 (MQ=255)
                                                                                                                                         ccataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGCGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGc              >  1:1048041/6‑150 (MQ=255)
                                                                                                                                         ccataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGCATTTACCGTGCAGGc              >  1:984313/6‑150 (MQ=255)
                                                                                                                                         ccataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGc              >  2:210086/6‑150 (MQ=255)
                                                                                                                                          tGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGcc             <  1:253748/150‑1 (MQ=255)
                                                                                                                                          tGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGcc             <  1:1017365/150‑1 (MQ=255)
                                                                                                                                          cataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGcc             >  1:945127/5‑150 (MQ=255)
                                                                                                                                          cataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGcc             >  1:1564149/5‑150 (MQ=255)
                                                                                                                                           ataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCa            >  2:685043/4‑150 (MQ=255)
                                                                                                                                               aCTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAg        <  2:541194/150‑1 (MQ=255)
                                                                                                                                                ctgcctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGc       >  2:2124440/1‑150 (MQ=255)
                                                                                                                                                 tgcctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCa      >  2:1605588/1‑150 (MQ=255)
                                                                                                                                                   cctgcGGCAGTTTATGGCCGAATGGTCTATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGCGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCagt    >  2:490529/1‑148 (MQ=255)
                                                                                                                                                    ctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTg   <  2:459774/150‑1 (MQ=255)
                                                                                                                                                     tgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGCACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTga  >  2:1551421/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGA  >  NC_000913/164889‑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: