Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,398,875 C→T M1M (GTG→ATG) †
*368* (TAG→TAA) 		 mreD ←
mreC ←		 cell wall structural complex MreBCD transmembrane component MreD
cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,398,8750CT100.0% 259.0 / NA 79M1M (GTG→ATG) †
*368* (TAG→TAA) 		mreD
mreC		cell wall structural complex MreBCD transmembrane component MreD
cell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (46/33);  total (46/33)

AATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCT  >  NC_000913/3398740‑3399008
                                                                                                                                       |                                                                                                                                     
aaTACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTAt                                                                                                                                    >  2:1369660/1‑139 (MQ=255)
  tACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTg                                                                                                                                  <  2:1821280/139‑1 (MQ=255)
       caaGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTc                                                                                                                             >  1:1011362/1‑139 (MQ=255)
       caaGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTc                                                                                                                             >  1:504715/1‑139 (MQ=255)
                 aaCACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCAcg                                                                                                                   >  2:1995983/1‑139 (MQ=255)
                     cccAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCATCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCa                                                                                                               >  2:1946338/1‑139 (MQ=255)
                        aGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGc                                                                                                            >  1:1678990/1‑139 (MQ=255)
                          tttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                          <  2:1225153/139‑1 (MQ=255)
                          tttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                          <  2:504715/139‑1 (MQ=255)
                               ccGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAg                                                                                                     >  1:1869091/1‑139 (MQ=255)
                                    aaaCAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGcacca                                                                                                >  2:1487712/1‑139 (MQ=255)
                                             ggTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggc                                                                                       >  2:1642641/1‑139 (MQ=255)
                                                        cAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCATCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCt                                                                            <  1:1946338/139‑1 (MQ=255)
                                                         aGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCTTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCtt                                                                           >  1:239657/1‑139 (MQ=255)
                                                                gATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGgcg                                                                                                                       >  2:535331/1‑88 (MQ=255)
                                                                gATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGgcg                                                                                                                       <  1:535331/88‑1 (MQ=255)
                                                                gATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAg                                                                                                     <  2:196539/106‑1 (MQ=255)
                                                                gATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAg                                                                                                     >  1:196539/1‑106 (MQ=255)
                                                                gATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGgag                                                                    >  2:1853928/1‑139 (MQ=255)
                                                                gATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCACAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGgag                                                                    >  1:1293282/1‑139 (MQ=255)
                                                                  tttGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCTTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGagag                                                                  <  2:239657/139‑1 (MQ=255)
                                                                          aaCAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagc                                                          <  2:1011362/139‑1 (MQ=255)
                                                                               cgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTAt                                                                                                                                    >  2:1830156/1‑60 (MQ=255)
                                                                               cgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTAt                                                                                                                                    <  1:1830156/60‑1 (MQ=255)
                                                                               cgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCAcgcg                                                                                                                 <  2:1288691/79‑1 (MQ=255)
                                                                               cgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCAcgcg                                                                                                                 >  1:1288691/1‑79 (MQ=255)
                                                                                      aGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTgtgt                                              >  2:836843/1‑139 (MQ=255)
                                                                                       ggAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                          >  2:1828887/1‑78 (MQ=255)
                                                                                       ggAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                          <  1:1828887/78‑1 (MQ=255)
                                                                                         aaaGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCgg                                           >  1:1966244/1‑139 (MQ=255)
                                                                                               agCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                             >  1:959508/1‑67 (MQ=255)
                                                                                               agCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                             <  2:959508/67‑1 (MQ=255)
                                                                                               agCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagc                                     <  2:1122181/139‑1 (MQ=255)
                                                                                               agCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagc                                     <  2:1966244/139‑1 (MQ=255)
                                                                                                gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGtt                                                                                                        >  2:1665565/1‑71 (MQ=255)
                                                                                                gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGtt                                                                                                        <  1:1665565/71‑1 (MQ=255)
                                                                                                gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAgcagca                                                      >  1:754549/1‑121 (MQ=255)
                                                                                                gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAgcagca                                                      <  2:754549/121‑1 (MQ=255)
                                                                                                gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagca                                    >  2:2016503/1‑139 (MQ=255)
                                                                                                 ccAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGtt                                                                                                        >  1:1105574/1‑70 (MQ=255)
                                                                                                 ccAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGtt                                                                                                        <  2:1105574/70‑1 (MQ=255)
                                                                                                     aTTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGgcgc                               <  1:1853928/139‑1 (MQ=255)
                                                                                                        aCCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGTGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGt                            <  1:1369660/139‑1 (MQ=255)
                                                                                                            aGCGTCCCTGGCTACGATAGCTCGCCATTATTGcc                                                                                                                                >  1:1601830/1‑35 (MQ=37)
                                                                                                            aGCGTCCCTGGCTACGATAGCTCGCCATTATTGcc                                                                                                                                <  2:1601830/35‑1 (MQ=37)
                                                                                                              cGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGc                      >  1:838038/1‑139 (MQ=255)
                                                                                                               gTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCa                     <  2:838038/139‑1 (MQ=255)
                                                                                                               gTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCa                     >  1:224872/1‑139 (MQ=255)
                                                                                                                   cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTc                 <  1:1395620/139‑1 (MQ=255)
                                                                                                                   cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTc                 >  1:229969/1‑139 (MQ=255)
                                                                                                                   cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTc                 >  2:1030283/1‑139 (MQ=255)
                                                                                                                   cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTc                 <  1:836843/139‑1 (MQ=255)
                                                                                                                    tGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTcc                >  1:301395/1‑139 (MQ=255)
                                                                                                                    tGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTcc                >  1:717912/1‑139 (MQ=255)
                                                                                                                          cGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTg              <  2:1492345/135‑1 (MQ=255)
                                                                                                                          cGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTg              >  1:1492345/1‑135 (MQ=255)
                                                                                                                           gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                        >  2:1665551/1‑60 (MQ=255)
                                                                                                                           gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                        <  1:1665551/60‑1 (MQ=255)
                                                                                                                           gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGcc         >  1:700328/1‑139 (MQ=255)
                                                                                                                           gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTCGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTa                                                            <  1:203727/88‑1 (MQ=255)
                                                                                                                           gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTCGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTa                                                            >  2:203727/1‑88 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                          <  2:1191606/37‑1 (MQ=38)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                          >  1:1191606/1‑37 (MQ=38)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagcag                                                        >  1:155974/1‑87 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagcag                                                        <  2:155974/87‑1 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagcag                                                        >  1:402964/1‑87 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagcag                                                        <  2:402964/87‑1 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagca                                                         <  2:1292513/86‑1 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagca                                                         >  1:1292513/1‑86 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAgcagca                                                      <  1:484856/89‑1 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAgcagca                                                      >  2:484856/1‑89 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  1:1131285/1‑139 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  1:1059387/1‑139 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  2:610407/1‑139 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  1:1982953/1‑139 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  1:285177/1‑139 (MQ=255)
                                                                                                                                cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg    >  1:305795/1‑139 (MQ=255)
                                                                                                                                  cGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtgct  <  1:1487712/139‑1 (MQ=255)
                                                                                                                                  cGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtgct  <  1:1995983/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
AATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCT  >  NC_000913/3398740‑3399008

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