Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,349,162 C→T 6.9% V216V (GTG→GTA) ‡ mtgA ← biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,349,1620CT 6.9% 407.6 / 11.0 130V216V (GTG→GTA) ‡mtgAbiosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (62/58);  new base T (6/3);  total (68/62)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

Only 100 of 132 total aligned reads displayed.

CTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCAAA  >  NC_000913/3349014‑3349311
                                                                                                                                                    |                                                                                                                                                     
cTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGcc                                                                                                                                                      >  1:994449/1‑149 (MQ=255)
 tATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCCGCTGGTGTTGTTGCATAAACGGCGCACCGCCTCACTGATACCTCTGCCGTAAAATCCCCGCCGGCCGGCTACGCAc                                                                                                                                                     >  1:1033452/1‑150 (MQ=255)
 tATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCAc                                                                                                                                                     >  1:547098/1‑150 (MQ=255)
  aTTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCAAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACg                                                                                                                                                    >  1:1646484/1‑150 (MQ=255)
    tGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTa                                                                                                                                                  >  2:1315460/1‑150 (MQ=255)
       ccTGCAGAAAGCTCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTCCGCACGTAGcc                                                                                                                                               >  1:457507/1‑150 (MQ=255)
       ccTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGcc                                                                                                                                               >  2:1676752/1‑150 (MQ=255)
       ccTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGcc                                                                                                                                               <  1:1142262/150‑1 (MQ=255)
       ccTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGcc                                                                                                                                               <  2:1298294/150‑1 (MQ=255)
       ccTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCGCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGcc                                                                                                                                               >  2:1113172/1‑150 (MQ=255)
        cTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGccc                                                                                                                                              <  2:1373340/150‑1 (MQ=255)
          gCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGTGCGTAGccgtc                                                                                                                                            <  2:1559059/150‑4 (MQ=255)
          gCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGa                                                                                                                                            >  2:453810/1‑150 (MQ=255)
          gCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGa                                                                                                                                            >  1:1027990/1‑150 (MQ=255)
           cAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACCTAGCCCGAt                                                                                                                                           >  2:727131/1‑150 (MQ=255)
              aaaGATCTGCAACCTCGTGAGAGCAAAACGCGGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGc                                                                                                                                        <  1:1655429/150‑1 (MQ=255)
               aaGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGTGCGTAGccgtcag                                                                                                                                          <  1:280754/147‑6 (MQ=255)
               aaGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGTGCGTAGccgtcag                                                                                                                                          >  2:791701/1‑142 (MQ=255)
                   tCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagga                                                                                                                                   <  2:367706/150‑1 (MQ=255)
                    cTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgaggag                                                                                                                                  <  2:178823/150‑1 (MQ=255)
                       cAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACt                                                                                                                               >  1:497398/1‑150 (MQ=255)
                           cTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGa                                                                                                                           <  1:285522/150‑1 (MQ=255)
                           cTCGTGAGAGCAAAACGCTGGCACTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGa                                                                                                                           <  1:1223430/150‑1 (MQ=255)
                             cGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGaaa                                                                                                                         <  1:1233461/150‑1 (MQ=255)
                             cGTGAGAGCAAAACGCTGGCCCTCACATCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGaaa                                                                                                                         <  1:345130/150‑1 (MQ=255)
                              gTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAAc                                                                                                                        <  2:181731/150‑1 (MQ=255)
                                    gCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTg                                                                                                                  >  2:1412257/1‑150 (MQ=255)
                                     cAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTgg                                                                                                                 >  2:434792/1‑150 (MQ=255)
                                     cAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACAACTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGACCGATGGCGAGGAGAATTTGAAACGAAGtct                                                                                                                 >  2:1565182/1‑148 (MQ=255)
                                           gCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAgg                                                                                                           >  1:533878/1‑150 (MQ=255)
                                           gCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAgg                                                                                                           <  2:1494512/150‑1 (MQ=255)
                                            cgggCCCTCACTTCGCGCGAAGCTTAATCCAGCTCGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGAAAAATCCACTCCTGACGTCTAATAACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGt                                                                                                          <  2:253069/148‑1 (MQ=255)
                                               gCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAAt                                                                                                       >  2:1174433/1‑150 (MQ=255)
                                                  cTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACa                                                                                                    >  1:191271/1‑150 (MQ=255)
                                                   tCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCTCGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAg                                                                                                   <  1:1618146/150‑1 (MQ=255)
                                                    cACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGc                                                                                                  <  2:1363281/150‑1 (MQ=255)
                                                    cACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGc                                                                                                  <  2:1544027/150‑1 (MQ=255)
                                                     aCTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCt                                                                                                 >  2:1326037/1‑150 (MQ=255)
                                                      cTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTg                                                                                                <  1:1579122/150‑1 (MQ=255)
                                                      cTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTg                                                                                                >  1:729823/1‑150 (MQ=255)
                                                       ttCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGc                                                                                               <  1:1163680/150‑1 (MQ=255)
                                                        tCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGcc                                                                                              <  2:509435/150‑1 (MQ=255)
                                                         cgcgcgAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTa                                                                                                         >  1:1702710/1‑138 (MQ=255)
                                                         cgcgcgAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTa                                                                                                         <  2:1702710/138‑1 (MQ=255)
                                                           cgcgAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGt                                                                                           >  2:384134/1‑150 (MQ=255)
                                                                  cTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAg                                                                                    >  2:1208114/1‑150 (MQ=255)
                                                                    tAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCt                                                                                  >  2:1503374/1‑150 (MQ=255)
                                                                    tAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCt                                                                                  <  2:1027990/150‑1 (MQ=255)
                                                                       tCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTcc                                                                               <  2:1575451/150‑1 (MQ=255)
                                                                       tCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTcc                                                                               <  2:1256589/150‑1 (MQ=255)
                                                                       tCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTcc                                                                               <  2:467283/150‑1 (MQ=255)
                                                                           gCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGAcc                                                                           >  1:1698300/1‑150 (MQ=255)
                                                                            cTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCa                                                                                             >  2:877954/1‑131 (MQ=255)
                                                                            cTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCa                                                                                             <  1:877954/131‑1 (MQ=255)
                                                                            cTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACAACTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCg                                                                          >  1:459409/1‑150 (MQ=255)
                                                                               gtgtTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGAcggtgt                                                                       >  2:932936/1‑145 (MQ=255)
                                                                                  ttgttgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCGGCTTCCGACCGGGTAAg                                                                    >  1:1667976/1‑150 (MQ=255)
                                                                                   tgttgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGt                                                                   >  2:354228/1‑150 (MQ=255)
                                                                                   tgttgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGt                                                                   >  2:969050/1‑150 (MQ=255)
                                                                                      tgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTg                                                                <  1:1208114/150‑1 (MQ=255)
                                                                                      tgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTg                                                                >  1:650822/1‑150 (MQ=255)
                                                                                        cATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGCCTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCt                                                              >  1:860755/1‑150 (MQ=255)
                                                                                        cATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCt                                                              <  1:1194537/150‑1 (MQ=255)
                                                                                        cATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCt                                                              <  1:384134/150‑1 (MQ=255)
                                                                                                 cTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGttt                                                     <  1:1326037/150‑1 (MQ=255)
                                                                                                 cTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGAAGGCTACGCACGTAGCCCGATGGCGAGGAGACATTGAAACGAAGTGGATTAGGTAATACAGATGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGttt                                                     >  1:472128/1‑150 (MQ=255)
                                                                                                 cTCACCGACTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATCCAGATGCCAGTAATGCAGCTTACTACCGGGTAAGTTAGATCGATGGttt                                                     >  2:914742/1‑150 (MQ=255)
                                                                                                  tCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtg                                                    <  2:1227461/150‑1 (MQ=255)
                                                                                                   cACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgt                                                   <  1:1676752/150‑1 (MQ=255)
                                                                                                   cACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgt                                                   >  1:502278/1‑150 (MQ=255)
                                                                                                    aCCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgtg                                                  >  1:1323956/1‑150 (MQ=255)
                                                                                                     ccGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTgg                                                 <  2:533878/150‑1 (MQ=255)
                                                                                                     ccGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTgg                                                 <  1:1174433/150‑1 (MQ=255)
                                                                                                       gCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGaa                                               >  2:513310/1‑150 (MQ=255)
                                                                                                        ccTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGaaa                                              >  1:685487/1‑150 (MQ=255)
                                                                                                            aCTGATACATCTGCCGTACAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTTATGCAGCTTCCGACCGTGTAAGTTTGCTCGCGGGTTTGTGGAAATAAc                                          <  2:272413/150‑1 (MQ=255)
                                                                                                             cTGCTACATCTGCCGTAAAATCCTCGTCTGACGGCTACGCACGTAGCCCGATGGCGATGAGACTTTGATACGTAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTTGAAATAACg                                         <  1:105735/150‑1 (MQ=255)
                                                                                                               gATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGtt                                       <  1:703036/150‑1 (MQ=255)
                                                                                                                 tACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTtgt                                     >  1:963147/1‑150 (MQ=255)
                                                                                                                  aCATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTtgtg                                    <  2:1646484/150‑1 (MQ=255)
                                                                                                                   cATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGc                                   >  2:1043524/1‑150 (MQ=255)
                                                                                                                   cATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGc                                   >  1:182207/1‑150 (MQ=255)
                                                                                                                   cATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGc                                   <  2:497398/150‑1 (MQ=255)
                                                                                                                       tGCCGTAAAATCCACGCCTGACGGCAACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCa                               <  1:456275/150‑1 (MQ=255)
                                                                                                                         ccGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGc                             >  1:85813/1‑150 (MQ=255)
                                                                                                                          cGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGATACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCt                            <  2:1648871/150‑1 (MQ=255)
                                                                                                                          cGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCt                            >  2:752271/1‑150 (MQ=255)
                                                                                                                           gTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCtt                           <  1:200576/150‑1 (MQ=255)
                                                                                                                               aaTCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGAc                       <  2:191271/150‑1 (MQ=255)
                                                                                                                                  ccACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGATATAACGTTGTGCCGCAGCTTCGACGcc                    >  1:1597880/1‑150 (MQ=255)
                                                                                                                                   cACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGTAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCa                   <  1:32446/150‑1 (MQ=255)
                                                                                                                                   cACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCa                   >  2:480171/1‑150 (MQ=255)
                                                                                                                                   cACCCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGACGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCa                   <  1:1094265/150‑1 (MQ=255)
                                                                                                                                     cGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCaaa                 >  2:445342/1‑150 (MQ=255)
                                                                                                                                        catcgggCTACGTGCGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAacac              >  2:86918/6‑150 (MQ=255)
                                                                                                                                        cTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAacac              >  2:196293/1‑150 (MQ=255)
                                                                                                                                          tcgggCTACGTGCGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAacacac            >  1:1313466/4‑150 (MQ=255)
                                                                                                                                               cTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCg       <  1:969050/150‑1 (MQ=255)
                                                                                                                                               cTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCg       >  2:762071/1‑150 (MQ=255)
                                                                                                                                               cTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCg       >  1:1255926/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCAAA  >  NC_000913/3349014‑3349311

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: