Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,104,887 G→A S152N (AGC→AAC)  mltC → membrane‑bound lytic murein transglycosylase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,104,8870GA96.5% 181.0 / NA 58S152N (AGC→AAC) mltCmembrane‑bound lytic murein transglycosylase C
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/23);  minor base T (0/2);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CAATATCAGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGCCGTTCCGATGCGCTGGGATTA  >  NC_000913/3104643‑3105131
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
cAATATCAGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGTTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACtt                                                                                                                                                                                                                                                <  2:312797/251‑1 (MQ=255)
 aagatCAGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGAGTTCGAAAGAACCTTGCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAATCAACTTc                                                                                                                                                                                                                                               <  2:160017/248‑1 (MQ=255)
      cAGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCgg                                                                                                                                                                                                                                           >  2:124170/1‑250 (MQ=255)
       aGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACATAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTGACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGTTCGCGCAAACAACTTCGCGGAt                                                                                                                                                                                                                                         <  1:428027/251‑1 (MQ=255)
                        aaCTTCGATGACGGTACGATTACTATCGACACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGACGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCg                                                                                                                                                                                                                        >  2:300118/1‑251 (MQ=255)
                          cTTCGATGACGGTACGATTACTCTCGAAACCATCGCCGGGACCGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTc                                                                                                                                                                                                                      >  1:117645/1‑251 (MQ=255)
                             cGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTTGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGa                                                                                                                                                                                                                   <  2:269663/251‑1 (MQ=255)
                                                  cGAAACCAGCGCCGGGACAGATCCCTCCGCGCATTTGCGCCGGGCAAGTATCAAAACGTTATTGATGGGGGACGATCCGAGTGCGGTCGATCTCTAGTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAGCTCCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTg                                                                                                                                                                                               <  2:355992/250‑1 (MQ=255)
                                                   gAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCg                                                                                                                                                                                             >  2:329556/1‑251 (MQ=255)
                                                           cGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATGTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTatcatc                                                                                                                                                                                      <  2:376682/250‑1 (MQ=255)
                                                                       aCCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCAcc                                                                                                                                                                          >  2:425550/1‑250 (MQ=255)
                                                                         cTGCCGCGCCTTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCAtt                                                                                                                                                                       >  2:64542/1‑251 (MQ=255)
                                                                                       cgcCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGa                                                                                                                                                          >  1:178405/1‑250 (MQ=255)
                                                                                         ccGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAacc                                                                                                                                                       >  1:402471/1‑251 (MQ=255)
                                                                                                  tATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGAt                                                                                                                                               >  1:15976/1‑250 (MQ=255)
                                                                                                     cAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAAc                                                                                                                                           >  1:74001/1‑251 (MQ=255)
                                                                                                     cAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAAc                                                                                                                                           >  2:317966/1‑251 (MQ=255)
                                                                                                                 gATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAat                                                                                                                               >  1:171321/1‑251 (MQ=255)
                                                                                                                  aTGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAata                                                                                                                              <  1:438180/251‑1 (MQ=255)
                                                                                                                       tGACAATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCg                                                                                                                         >  1:229019/1‑251 (MQ=255)
                                                                                                                        gACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCgg                                                                                                                        <  2:108745/251‑1 (MQ=255)
                                                                                                                          cGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGCCAACAACGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGc                                                                                                                       >  1:260240/1‑250 (MQ=255)
                                                                                                                                gAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTcc                                                                                                                <  2:15976/251‑1 (MQ=255)
                                                                                                                                     cGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGTTGGACAACGCCGGGCAGCCTATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCGACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAg                                                                                                           <  1:72816/251‑1 (MQ=255)
                                                                                                                                         cGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCg                                                                                                        <  1:86280/250‑1 (MQ=255)
                                                                                                                                             ctctATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCCGGCGTCACg                                                                                                   >  2:139243/1‑251 (MQ=255)
                                                                                                                                               ctATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCAc                                                                                                    <  2:171321/248‑1 (MQ=255)
                                                                                                                                               ctATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTAAGGTGGTGGACACCACCGGGCAGCCGATTAGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGAGTATCATCAACAGCGTCACCATTAACATGGAGCCGAACCACCTTGATACACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACgg                                                                                                  >  1:417694/1‑250 (MQ=255)
                                                                                                                                                       aCGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGc                                                                                         >  1:329541/1‑251 (MQ=255)
                                                                                                                                                       aCGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGc                                                                                         <  1:64542/251‑1 (MQ=255)
                                                                                                                                                           tgatgaTATGACGATGTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTTGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAATCAGCTTCGCGGATTATCTGCTGAAAAACCTTCTGAAGAGACGCAGCAACTGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCTCCTTGATAAACTTGCGCACAAATATCTCGGCATGGCCCGCCAGGCGTTACTGAAATATGGCGtt                                                                                      <  2:409644/250‑1 (MQ=255)
                                                                                                                                                                atatTACGATTTAGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCAGATTAGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTtgatg                                                                                  >  2:227723/1‑249 (MQ=255)
                                                                                                                                                                      cGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTg                                                                          >  1:101829/1‑251 (MQ=255)
                                                                                                                                                                       gATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGAGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTg                                                                          >  1:355194/1‑250 (MQ=255)
                                                                                                                                                                              aaaGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGc                                                                  <  1:317966/251‑1 (MQ=255)
                                                                                                                                                                                aGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCa                                                                 >  1:72685/1‑250 (MQ=255)
                                                                                                                                                                                          ctttACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGAc                                                      >  1:249483/1‑251 (MQ=255)
                                                                                                                                                                                             tACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCACCGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCg                                                    >  2:178880/1‑250 (MQ=255)
                                                                                                                                                                                             tACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCTCCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCg                                                    <  1:213233/250‑1 (MQ=255)
                                                                                                                                                                                             tACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCg                                                    >  2:47237/1‑250 (MQ=255)
                                                                                                                                                                                              aCGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATTCAGACCGaa                                                  <  1:308710/251‑1 (MQ=255)
                                                                                                                                                                                              aCGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGa                                                   >  1:155860/1‑250 (MQ=255)
                                                                                                                                                                                                ggTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGCTGAGTCGCTGATTCTGGCAATTATGCAGACCGAATc                                                >  1:171940/1‑251 (MQ=255)
                                                                                                                                                                                                 gTCAGGTGGAGGACAACACCGGGCAGCCGCTTCGCTGGGAAGGTCGCGCAAACAACTTAGCGGATTATCTGCTGAAAAACCTTCTGCAGAGCCGCAGCAACGGGCTGCGTATCCTCTACAGCGTCCCCATTAACATGGTGCCGACCCACCTTGATAAACGTACGCCCAACTATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGCTGAGTCGCTGATTCTGGCAATTATGCAGACCGCATCt                                               >  2:438644/1‑251 (MQ=255)
                                                                                                                                                                                                            gACAACACCGGGCAGCCGATTCGCTGGGAAGGTCTCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAAccc                                    >  2:111009/1‑251 (MQ=255)
                                                                                                                                                                                                               aacaCCGGGCAGCCGATTCGCTGGGAAGGTCGCACAAACAACTTCGCGGATTAGCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATATACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCGCGGAAATATGGCGTTTATGAGTCGCTGATTCTGGCAATTATGCTGACCGAATCTTCCTTTAACCCGTa                                 <  2:310827/251‑1 (MQ=255)
                                                                                                                                                                                                                acacCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTAt                                <  1:327360/251‑1 (MQ=255)
                                                                                                                                                                                                                   ccGGGCAGCCGATTCGATGGGAAGGTCGCGCAAACAACTTCGCCGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAAAGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCg                             >  1:173089/1‑251 (MQ=255)
                                                                                                                                                                                                                         aGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGt                                                                               <  1:216607/195‑1 (MQ=255)
                                                                                                                                                                                                                         aGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGt                                                                               >  2:216607/1‑195 (MQ=255)
                                                                                                                                                                                                                                  gCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGCCGTTCCGAt              <  1:75006/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gggAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATTAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGCCGTTCCGATgc            <  1:281341/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aGGTCGCGCAAACAACTACGCGGATTATCTGCTGAAAAAACGTCTGAAGAGCCGCAGCAACGGGCTGCGCCTCATCTACAGCGTCAACATTAACATGGTGCCGAACCACCTTGCTAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGCTTCTGGCACTTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGACGTTCCGATGCGCTgg       >  2:210193/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ggTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATg                               >  1:383330/1‑226 (MQ=255)
                                                                                                                                                                                                                                          ggTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATg                               <  2:383330/226‑1 (MQ=255)
                                                                                                                                                                                                                                           gTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGCCGTTCCGATGCGCTGGGa     <  1:47237/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cgcgCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAAGTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGCCGTTCCGATGCGCTGGGAtt   <  2:101829/251‑1 (MQ=255)
                                                                                                                                                                                                                                               cgcAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGCCGTTCCGATGCGCTGGGATTa  >  1:200349/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
CAATATCAGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGCCGTTCCGATGCGCTGGGATTA  >  NC_000913/3104643‑3105131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑