Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,349,137 G→A 6.7% R225W (CGG→TGG)  mtgA ← biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,349,1370GA 6.7% 310.8 / 9.2 90R225W (CGG→TGG) mtgAbiosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (39/45);  new base A (3/3);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTG  >  NC_000913/3349000‑3349275
                                                                                                                                         |                                                                                                                                          
ccAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTa                                                                                                                                          <  1:1174960/140‑1 (MQ=255)
ccAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTa                                                                                                                                          >  2:1174960/1‑140 (MQ=255)
ccAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCATa                                                                                                                                          >  2:525476/1‑140 (MQ=255)
ccAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCATa                                                                                                                                          >  2:525105/1‑140 (MQ=255)
ccAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCATa                                                                                                                                          >  2:523816/1‑140 (MQ=255)
     cGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGATCACCGCCTAACTGATACATCTGCCGTAAAATc                                                                                                                                     <  2:1240192/140‑1 (MQ=255)
      ggCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATcc                                                                                                                                    <  1:1021440/140‑1 (MQ=255)
       gcgcTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCa                                                                                                                                   <  2:1392074/140‑1 (MQ=255)
            ggCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCt                                                                                                                              <  2:956412/140‑1 (MQ=255)
                    cccTGCAGCAAGATCTGCACCCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTACTCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTa                                                                                                                      <  1:516014/140‑1 (MQ=255)
                    cccTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTa                                                                                                                      <  1:516024/140‑1 (MQ=255)
                    cccTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTa                                                                                                                      <  1:289730/140‑1 (MQ=255)
                       tGCAGAATGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTGGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACTGCTACGc                                                                                                                   >  2:1436993/1‑140 (MQ=255)
                       tGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGc                                                                                                                   >  1:1480957/1‑140 (MQ=255)
                       tGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGc                                                                                                                   >  1:1533190/1‑140 (MQ=255)
                       tGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGc                                                                                                                   >  1:1300900/1‑140 (MQ=255)
                       tGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGc                                                                                                                   >  2:135803/1‑140 (MQ=255)
                       tGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGGAAAATCCACGCCTGACGGCTACGc                                                                                                                   >  2:1234522/1‑140 (MQ=255)
                       tGCAGAAAGATCTGCAACCACGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCAGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTCACGGCTACGc                                                                                                                   >  2:1508006/1‑140 (MQ=255)
                         cAGAAAGATCTGCGACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCAc                                                                                                                 >  2:840731/1‑140 (MQ=255)
                         cAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCAc                                                                                                                 >  1:564984/1‑140 (MQ=255)
                         cAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCAc                                                                                                                 >  2:300437/1‑140 (MQ=255)
                         cAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCAc                                                                                                                 >  2:239170/1‑140 (MQ=255)
                         cAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCAc                                                                                                                 >  2:1526787/1‑140 (MQ=255)
                         cAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGAAACATCTGCCGTAAAATCCACGCCTGACGGCTACGCAc                                                                                                                 >  1:564398/1‑140 (MQ=255)
                          aGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACg                                                                                                                <  1:1334159/140‑1 (MQ=255)
                          aGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACg                                                                                                                <  1:139368/140‑1 (MQ=255)
                          aGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACg                                                                                                                <  2:87467/140‑1 (MQ=255)
                          aGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACg                                                                                                                <  2:1101011/140‑1 (MQ=255)
                          aGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACg                                                                                                                <  1:138948/140‑1 (MQ=255)
                                   tGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGAt                                                                                                       >  1:1107665/1‑140 (MQ=255)
                                   tGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGAt                                                                                                       <  1:1020449/140‑1 (MQ=255)
                                   tGCAACCGCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCGACGCACGTAGCCCGAt                                                                                                       >  1:1143130/1‑140 (MQ=255)
                                   agCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTAACTCCAGCTGGAGTAGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGAt                                                                                                       <  1:784817/139‑1 (MQ=255)
                                     cAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATgg                                                                                                     <  2:635032/140‑1 (MQ=255)
                                           cgcgagaGCAAAACGCTGGCCCACACCACGCGCGAAGCTCACCCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagga                                                                                               <  1:60783/137‑1 (MQ=255)
                                           cGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagga                                                                                               <  1:784673/140‑1 (MQ=255)
                                           cGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagga                                                                                               <  1:783139/140‑1 (MQ=255)
                                                     aaaCGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGaaa                                                                                     >  1:1637742/1‑140 (MQ=255)
                                                      aaCGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAAc                                                                                    >  2:772290/1‑140 (MQ=255)
                                                        cGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAATCGa                                                                                  <  2:1052721/140‑1 (MQ=255)
                                                        cGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGa                                                                                  <  2:1053163/140‑1 (MQ=255)
                                                         gCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGaa                                                                                 <  2:920586/140‑1 (MQ=255)
                                                           tGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGt                                                                               >  1:405867/1‑140 (MQ=255)
                                                             gCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTgg                                                                             >  1:534707/1‑140 (MQ=255)
                                                              cccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  2:675411/140‑1 (MQ=255)
                                                              cccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  1:236235/140‑1 (MQ=255)
                                                              cccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  1:416376/140‑1 (MQ=255)
                                                              cccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  2:896364/140‑1 (MQ=255)
                                                              cccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  1:527615/140‑1 (MQ=255)
                                                              cccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  1:416390/140‑1 (MQ=255)
                                                              cccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  1:500926/140‑1 (MQ=255)
                                                              cccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  1:527929/140‑1 (MQ=255)
                                                              cccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  1:928900/140‑1 (MQ=255)
                                                              cccTCAATTCGCTCGAAGCATTATCCAGCTGGTGTAGTTGCATAAACGGCTCACCGCCTAACTTATACATCTGCCGTAAAATCCACGACTGACGGCTACGCACGATGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                            <  2:1667326/140‑1 (MQ=255)
                                                                cTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGAtt                                                                          <  2:534707/140‑1 (MQ=255)
                                                                cTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGAtt                                                                          <  1:979317/140‑1 (MQ=255)
                                                                         cgcgAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATAc                                                                 >  1:546353/1‑140 (MQ=255)
                                                                         cgcgAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATAc                                                                 >  1:1328262/1‑140 (MQ=255)
                                                                         cgcgAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATAc                                                                 >  1:1391492/1‑140 (MQ=255)
                                                                         cgcgAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATAc                                                                 >  1:810252/1‑140 (MQ=255)
                                                                             aaGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCt                                                             <  1:518231/140‑1 (MQ=255)
                                                                             aaGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCt                                                             <  1:1656916/140‑1 (MQ=255)
                                                                               gCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGc                                                           <  1:1582426/140‑1 (MQ=255)
                                                                                   aaTCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGt                                                       >  2:1416921/1‑140 (MQ=255)
                                                                                          cTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGTTAATACAGCTGCCAGTAATGCAg                                                >  2:1667128/1‑140 (MQ=255)
                                                                                          cTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAg                                                >  1:622848/1‑140 (MQ=255)
                                                                                          cTGGTGTTGTTGCATAAACGGCTCACCGCCTAACAGATACAACTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAg                                                >  1:1409824/1‑140 (MQ=255)
                                                                                                 tgttgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGa                                         <  2:1260243/140‑1 (MQ=255)
                                                                                                 tgttgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGa                                         <  2:266383/140‑1 (MQ=255)
                                                                                                                  aCCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTcg                        <  2:1070893/140‑1 (MQ=255)
                                                                                                                         aaCTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGttt                 >  2:638875/1‑140 (MQ=255)
                                                                                                                         aaCTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGttt                 >  2:1224775/1‑140 (MQ=255)
                                                                                                                         aaCTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGttt                 >  1:1481175/1‑140 (MQ=255)
                                                                                                                         aaCTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGttt                 >  2:1627318/1‑140 (MQ=255)
                                                                                                                         aaCTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGttt                 >  2:181657/1‑140 (MQ=255)
                                                                                                                           cTGATTAATCTGCCGTATAATCCACGCCTGACGGCTACGCTCGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgt               <  1:1369601/140‑1 (MQ=255)
                                                                                                                           cTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgt               <  1:1452652/140‑1 (MQ=255)
                                                                                                                           cTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgt               <  1:1197514/140‑1 (MQ=255)
                                                                                                                              aTACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGa            >  1:736088/1‑140 (MQ=255)
                                                                                                                                aCATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTACGTTTGCTCGCGGGTTTGTGGaaa          >  2:561678/1‑140 (MQ=255)
                                                                                                                                aCATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGaaa          <  2:766680/140‑1 (MQ=255)
                                                                                                                                aCATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGAATAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGaaa          <  2:1264306/140‑1 (MQ=255)
                                                                                                                                aCATCTGCCGTAAAATCCACGCCGGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGATTGCTCGCGGGTTTGTGGaat          >  2:587719/1‑139 (MQ=255)
                                                                                                                                aCATCTGCCGTAAAATACACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGaaa          >  2:561709/1‑140 (MQ=255)
                                                                                                                                     tGCCATAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACg     <  1:523816/140‑1 (MQ=255)
                                                                                                                                     tGCCATAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACg     <  1:525105/140‑1 (MQ=255)
                                                                                                                                     tGCCATAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACg     <  1:525476/140‑1 (MQ=255)
                                                                                                                                      gCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGt    >  2:239029/1‑140 (MQ=255)
                                                                                                                                        cGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTtg  >  2:899629/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTG  >  NC_000913/3349000‑3349275

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: