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%	261.3 / NA	67	E254K (GAA→AAA)	mrdB	cell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/31);  total (36/31)

AAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAATC  >  NC_000913/665409‑665701                                                                                                                                                  |                                                                                                                                                    aaGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCt                                                                                                                                                >  2:230485/1‑151 (MQ=255)  aGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCtg                                                                                                                                               >  1:160955/1‑151 (MQ=255)       gtTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACt                                                                                                                                          <  2:1072842/151‑1 (MQ=255)       gtTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACt                                                                                                                                          >  2:376537/1‑151 (MQ=255)          tGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTg                                                                                                                                         >  2:393470/1‑149 (MQ=255)           gcgcTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGa                                                                                                                                        <  2:1144272/149‑1 (MQ=255)           gcgcTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGCGACTGAg                                                                                                                                       >  1:318278/1‑150 (MQ=255)             gcTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGc                                                                                                                                    <  1:1052656/151‑1 (MQ=255)                   gcggtTACCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAACCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCg                                                                                                                                  <  1:1132521/147‑1 (MQ=255)                            cacaGACCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCtt                                                                                                                     >  1:916041/1‑151 (MQ=255)                               aGCCCGCGCATGATCAGCATAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGc                                                                                                                  <  2:1040219/151‑1 (MQ=255)                               aGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGc                                                                                                                  >  1:756787/1‑151 (MQ=255)                                gCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCTTGCAGCCAGCCTTTGcc                                                                                                                 >  2:73800/1‑151 (MQ=255)                                    gcgcATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGt                                                                                                             >  2:387494/1‑151 (MQ=255)                                    gcgcATGATCACCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCAATAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCATTTGCCGCGt                                                                                                             <  1:1027939/151‑1 (MQ=255)                                        aTGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATc                                                                                                         <  1:266794/151‑1 (MQ=255)                                          gATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATcc                                                                                                        <  2:917701/150‑1 (MQ=255)                                              agcagAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATccgcc                                                                                                     <  2:340145/149‑1 (MQ=255)                                               gcagAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAGTCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGa                                                                                                  <  1:399079/151‑1 (MQ=255)                                                             gcgagcAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAAtt                                                                                     <  1:24626/150‑1 (MQ=255)                                                                 gcAGAATCAGAATGCCCACTAACCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGa                                                                                <  2:977824/151‑1 (MQ=255)                                                                   agaatcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACt                                                                              >  1:651790/1‑151 (MQ=255)                                                                    gaatcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTg                                                                             <  2:362863/151‑1 (MQ=255)                                                                      atcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaa                                                                           <  2:832407/151‑1 (MQ=255)                                                                        cagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaat                                                                          >  2:800191/1‑150 (MQ=255)                                                                           aatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaataat                                                                       <  1:881332/150‑1 (MQ=255)                                                                           aatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaata                                                                         >  2:1168899/1‑148 (MQ=255)                                                                           aatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAAtg                                                                      <  1:779547/151‑1 (MQ=255)                                                                               cccACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGa                                                                   <  2:72738/150‑1 (MQ=255)                                                                                   cTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGc                                                               <  1:831359/150‑1 (MQ=255)                                                                                     aaTCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGcc                                                              <  2:111374/149‑1 (MQ=255)                                                                                         ccAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGcc                                                        >  2:1048672/1‑151 (MQ=255)                                                                                           aGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCg                                                       >  1:607766/1‑150 (MQ=255)                                                                                              tctTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAGTTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATACTGTGATAGCCCGCGCCg                                                       >  1:1013017/1‑147 (MQ=255)                                                                                               ctTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGt                                                    <  1:655671/149‑1 (MQ=255)                                                                                                ttCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTggg                                                 >  2:905763/1‑151 (MQ=255)                                                                                                 tCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTgg                                                  <  1:720260/149‑1 (MQ=255)                                                                                                  ccGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTc                                               >  1:244462/1‑151 (MQ=255)                                                                                                     ccAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTg                                             >  1:632912/1‑150 (MQ=255)                                                                                                      cAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGAt                                           >  1:61515/1‑151 (MQ=255)                                                                                                       aGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGAtt                                          >  1:323681/1‑151 (MQ=255)                                                                                                       aGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGAtt                                          >  2:512385/1‑151 (MQ=255)                                                                                                       aGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGAtt                                          <  2:724426/151‑1 (MQ=255)                                                                                                        gTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTc                                         >  1:262590/1‑151 (MQ=255)                                                                                                         tACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTcc                                        >  2:537464/1‑151 (MQ=255)                                                                                                          aCCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCg                                       <  2:1070098/151‑1 (MQ=255)                                                                                                             gcgAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGt                                    <  1:704009/151‑1 (MQ=255)                                                                                                                     aaaGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTTCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCAAGTGGGTCTGATTCCGGGTCCAGGAg                             <  1:520197/150‑1 (MQ=255)                                                                                                                     aaaGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGc                            >  2:229078/1‑151 (MQ=255)                                                                                                                     aaaGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGc                            >  1:1004696/1‑151 (MQ=255)                                                                                                                     aaaGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGc                            >  1:950556/1‑151 (MQ=255)                                                                                                                         tCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGcatc                         >  2:660694/1‑150 (MQ=255)                                                                                                                           aGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCAtt                      >  2:941493/1‑151 (MQ=255)                                                                                                                             tATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTAc                    >  1:937964/1‑151 (MQ=255)                                                                                                                              aTGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTAcg                   >  1:1101313/1‑151 (MQ=255)                                                                                                                              aTGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTAc                    <  2:377172/150‑1 (MQ=255)                                                                                                                               tGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTAcgc                  >  1:968046/1‑151 (MQ=255)                                                                                                                               tGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTAcg                   <  1:574652/150‑1 (MQ=255)                                                                                                                                  cGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgct               >  2:863752/1‑151 (MQ=255)                                                                                                                                   gTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctg              <  2:799768/151‑1 (MQ=255)                                                                                                                                   gTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctg              >  2:109010/1‑151 (MQ=255)                                                                                                                                   gTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctg              >  1:343241/1‑151 (MQ=255)                                                                                                                                          gagaAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGcgctggcgctgg       <  1:1165272/151‑1 (MQ=255)                                                                                                                                          gagaAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGcgctggcgctgg       >  1:213771/1‑151 (MQ=255)                                                                                                                                             aaaTTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTaa    <  1:230472/151‑1 (MQ=255)                                                                                                                                             aaaTTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTa     >  1:361884/1‑150 (MQ=255)                                                                                                                                               aTTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAATc  <  1:560415/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAATC  >  NC_000913/665409‑665701

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