Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	665,554	C→T	E254K (GAA→AAA)	mrdB ←	cell wall shape‑determining protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	665,554	0	C	T	100.0%	239.3 / NA	62	E254K (GAA→AAA)	mrdB	cell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/31);  total (31/31)

TTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAAT  >  NC_000913/665416‑665700                                                                                                                                           |                                                                                                                                                   tttGCGCTCTGGCGGCTTTCCCCAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACt                                                                                                                                         <  2:869323/150‑1 (MQ=255)    gcgcTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGt                                                                                                                                     >  1:797910/1‑151 (MQ=255)      gcTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGGAGAGCGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTg                                                                                                                                    >  1:132202/1‑150 (MQ=255)           ggcggcTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGc                                                                                                                              >  1:647064/1‑151 (MQ=255)             cggcTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGc                                                                                                                              >  2:297096/1‑149 (MQ=255)               gcTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGcagcc                                                                                                                          <  2:913671/151‑1 (MQ=255)               gcTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGcagcc                                                                                                                          >  1:781560/1‑151 (MQ=255)                cTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGcagcc                                                                                                                          >  1:536112/1‑150 (MQ=255)                        aGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGc                                                                                                                 <  1:37134/151‑1 (MQ=255)                          cccGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCcg                                                                                                               <  1:97463/151‑1 (MQ=255)                           ccGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCcgc                                                                                                              >  1:529876/1‑151 (MQ=255)                                cATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAAt                                                                                                         >  2:998756/1‑151 (MQ=255)                                cATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAAt                                                                                                         >  2:972783/1‑151 (MQ=255)                                 aTGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATc                                                                                                        >  1:268545/1‑151 (MQ=255)                                 aTGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATc                                                                                                        <  2:488205/151‑1 (MQ=255)                                 aTGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATc                                                                                                        >  2:218902/1‑151 (MQ=255)                                   gATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATcc                                                                                                       <  1:727895/150‑1 (MQ=255)                                     tCAGCAGAAGGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATccg                                                                                                      <  2:172952/149‑1 (MQ=255)                                      cagcagAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATccgcc                                                                                                    <  1:50967/150‑1 (MQ=255)                                      cagcagAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATccgcc                                                                                                    <  2:463112/150‑1 (MQ=255)                                        gcagAATGTAGAGAGCGAGCAGAATCAGAATTCCCACTAATCCCAGCCCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATccgccg                                                                                                   <  2:1017868/149‑1 (MQ=255)                                        gcagAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCgg                                                                                                  >  2:107408/1‑150 (MQ=255)                                          agAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGc                                                                                               >  2:1021633/1‑151 (MQ=255)                                            aaTGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGcc                                                                                              <  1:955798/150‑1 (MQ=255)                                             aTGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCa                                                                                             >  2:382196/1‑150 (MQ=255)                                                   agagCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCa                                                                                       <  2:261093/150‑1 (MQ=255)                                                   agagCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCa                                                                                       <  2:408160/150‑1 (MQ=255)                                                               atcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaa                                                                          <  1:364698/151‑1 (MQ=255)                                                               atcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGa                                                                           >  1:7000/1‑150 (MQ=255)                                                                 cagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaata                                                                        <  1:910734/151‑1 (MQ=255)                                                                   gaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaataat                                                                      >  2:620378/1‑151 (MQ=255)                                                                   gaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaataat                                                                      >  1:539044/1‑151 (MQ=255)                                                                    aatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaata                                                                        <  1:403592/148‑1 (MQ=255)                                                                         ccACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGAt                                                                 >  1:1051371/1‑150 (MQ=255)                                                                          cACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAAtgt                                                                    >  1:490163/1‑146 (MQ=255)                                                                              aaaccgaggtctTCCTCCAGTACCGCGGAGCTAAAGGCAGTAGAGCTGTCGGGGAGAAATTTAAGCCGTGACTGAGTACCGGGCAAACAGCCTGTGCCGCGTAATCCGCCGGAGCCAATAGGAAATTTAGACTTAATTATGTGAGAGccc                                                            >  2:116726/10‑150 (MQ=255)                                                                                tCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCcgcg                                                         >  2:599622/1‑151 (MQ=255)                                                                                   cAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCg                                                      <  1:706300/151‑1 (MQ=255)                                                                                         ttCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTggg                                                >  2:204301/1‑151 (MQ=255)                                                                                          tCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGt                                               >  1:142398/1‑151 (MQ=255)                                                                                           ccGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGt                                               <  1:105743/150‑1 (MQ=255)                                                                                                 gTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTc                                        <  2:280228/151‑1 (MQ=255)                                                                                                   aCCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCg                                      <  1:326949/151‑1 (MQ=255)                                                                                                    ccGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCgg                                     <  2:546622/151‑1 (MQ=255)                                                                                                      gcgAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGt                                   <  1:732471/151‑1 (MQ=255)                                                                                                      gcgAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGt                                   <  2:354663/151‑1 (MQ=255)                                                                                                      gcgAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGt                                   >  2:4122/1‑151 (MQ=255)                                                                                                       cgAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTc                                  >  1:997770/1‑151 (MQ=255)                                                                                                       cgAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTc                                  <  2:678935/151‑1 (MQ=255)                                                                                                                   cAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGcatca                       <  2:75502/150‑1 (MQ=255)                                                                                                                       aTGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTAc                   <  1:816290/150‑1 (MQ=255)                                                                                                                       aTGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTAc                   <  1:113009/150‑1 (MQ=255)                                                                                                                            gTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctg             <  1:996264/151‑1 (MQ=255)                                                                                                                            gTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctg             >  2:1079138/1‑151 (MQ=255)                                                                                                                                ggggAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctggcgc         >  2:1105014/1‑151 (MQ=255)                                                                                                                                ggggAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctggcgc         >  2:1051004/1‑151 (MQ=255)                                                                                                                                   gagaAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGcgctggcgctgg      <  1:297108/151‑1 (MQ=255)                                                                                                                                   gagaAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGcgctggcgctgg      >  1:1050586/1‑151 (MQ=255)                                                                                                                                     gaAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGt     >  2:715186/1‑150 (MQ=255)                                                                                                                                     gaAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTa    <  1:138373/151‑1 (MQ=255)                                                                                                                                      aaaTTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTa    >  2:664366/1‑150 (MQ=255)                                                                                                                                       aaTTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAAt  <  1:1098324/151‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                   TTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAAT  >  NC_000913/665416‑665700

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