Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,766,316	G→A	*1378* (TAG→TAA)	rhsA →	Rhs protein with putative toxin 55 domain, putative polysaccharide synthesis/export protein, putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,766,316	0	G	A	100.0%	272.6 / NA	86	*1378* (TAG→TAA)	rhsA	Rhs protein with putative toxin 55 domain, putative polysaccharide synthesis/export protein, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/47);  total (39/47)

ATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGAT  >  NC_000913/3766180‑3766435                                                                                                                                         |                                                                                                                        aTGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAAcc                                                                                                                       >  2:340804/1‑139 (MQ=255)  tGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAAccc                                                                                                                      >  1:166476/1‑139 (MQ=255)      aaaaaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtgt                                                                                                                  >  1:1008890/1‑139 (MQ=255)       aaaaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtgtg                                                                                                                 <  2:740910/139‑1 (MQ=255)         aaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                                               >  1:1114711/1‑139 (MQ=255)            cccGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGt                                                                                                            >  2:718268/1‑139 (MQ=255)                   gATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     <  2:166476/139‑1 (MQ=255)                   gATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     <  2:700910/139‑1 (MQ=255)                         aCGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCa                                                                                               >  1:949792/1‑139 (MQ=255)                               ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                                               <  1:1006706/117‑1 (MQ=255)                               ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTAGTGGa                                                                                                               >  2:1006706/1‑117 (MQ=255)                                     aaCTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCag                                                                                   >  2:675357/1‑139 (MQ=255)                                       cTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaa                                                                                 >  2:1210274/1‑139 (MQ=255)                                         aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                                               <  2:1314992/107‑1 (MQ=255)                                         aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                                               >  1:1314992/1‑107 (MQ=255)                                             gattgTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                                           <  1:340804/139‑1 (MQ=255)                                                tgtgGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                        <  1:718268/139‑1 (MQ=255)                                                 gtggtgATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCa                                                                       <  2:584235/139‑1 (MQ=255)                                                  tggtgATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCaa                                                                      <  2:1170044/139‑1 (MQ=255)                                                         atTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAataca                                                                                         <  1:223540/113‑1 (MQ=255)                                                         atTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAataca                                                                                         >  2:223540/1‑113 (MQ=255)                                                         atTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                                                 >  1:647212/1‑137 (MQ=255)                                                         atTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                                                 <  2:647212/137‑1 (MQ=255)                                                          ttAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     >  2:384873/1‑100 (MQ=255)                                                          ttAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     <  1:384873/100‑1 (MQ=255)                                                           tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                             <  2:1008890/139‑1 (MQ=255)                                                           tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                             <  2:827127/139‑1 (MQ=255)                                                           tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                             <  2:1114711/139‑1 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGTCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCa                                                                       >  1:215638/1‑128 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGTCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCa                                                                       <  2:215638/128‑1 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAataca                                                                                         <  1:1257097/110‑1 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAataca                                                                                         >  2:1257097/1‑110 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTa                                                            >  2:11364/1‑139 (MQ=255)                                                             gTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAAt                                                                                                   >  1:805971/1‑99 (MQ=255)                                                             gTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAAt                                                                                                   <  2:805971/99‑1 (MQ=255)                                                             gTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGt                                                                                                 >  2:213426/1‑101 (MQ=255)                                                             gTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGt                                                                                                 <  1:213426/101‑1 (MQ=255)                                                               gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                                           >  1:1068848/1‑121 (MQ=255)                                                               gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                                           <  2:1068848/121‑1 (MQ=255)                                                               gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATgg                                                         <  1:1210274/139‑1 (MQ=255)                                                               gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGTCAAGTAAAGAGTATgg                                                         <  2:52586/139‑1 (MQ=255)                                                                cTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAgg                                                                         <  1:198242/122‑1 (MQ=255)                                                                cTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAgg                                                                         >  2:198242/1‑122 (MQ=255)                                                                cTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGt                                                        <  2:444166/139‑1 (MQ=255)                                                                 tAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                             >  2:53830/1‑133 (MQ=255)                                                                 tAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                             <  1:53830/133‑1 (MQ=255)                                                                  aaaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     >  1:879230/1‑92 (MQ=255)                                                                  aaaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     <  2:879230/92‑1 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                        >  1:668831/1‑120 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                        <  2:668831/120‑1 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                                                    >  1:861002/1‑124 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                                                    <  2:861002/124‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAATCGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                        <  2:916225/115‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     <  2:212725/86‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     >  1:212725/1‑86 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaa                                                                                <  1:847749/107‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                        >  1:916225/1‑115 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACAGACCagaaa                                                                                >  2:847749/1‑107 (MQ=255)                                                                         cAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                                               <  1:11364/139‑1 (MQ=255)                                                                           aaaaGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCa                                             <  2:361437/139‑1 (MQ=255)                                                                                   acaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAgg                                                                         <  2:613233/103‑1 (MQ=255)                                                                                   acaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAgg                                                                         >  1:613233/1‑103 (MQ=255)                                                                                   acaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAGTCAATGTAAGa                                     >  1:875225/1‑139 (MQ=255)                                                                                               gAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                             <  2:1298669/103‑1 (MQ=255)                                                                                               gAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                             >  1:1298669/1‑103 (MQ=255)                                                                                               gAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTta                                                     <  2:166389/111‑1 (MQ=255)                                                                                               gAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTta                                                     >  1:166389/1‑111 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                >  1:903362/1‑62 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                <  2:903362/62‑1 (MQ=255)                                                                                                             ccAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGc           <  1:244479/139‑1 (MQ=255)                                                                                                                   aGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     <  2:481873/43‑1 (MQ=255)                                                                                                                   aGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     >  1:481873/1‑43 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                             >  2:3998/1‑81 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                             <  1:3998/81‑1 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                        <  2:618255/69‑1 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                        >  1:618255/1‑69 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                                                 <  2:1160545/76‑1 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                                                 >  1:1160545/1‑76 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATgg                                                         >  1:302589/1‑84 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATgg                                                         <  2:302589/84‑1 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTtgat  <  2:1094966/139‑1 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTtgat  <  2:1106208/139‑1 (MQ=255)                                                                                                                        aaCGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     >  1:378187/1‑38 (MQ=255)                                                                                                                        aaCGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                     <  2:378187/38‑1 (MQ=255)                                                                                                                             tAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                <  2:1225443/38‑1 (MQ=255)                                                                                                                             tAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                >  1:1225443/1‑38 (MQ=21)                                                                                                                                         |                                                                                                                        ATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGAT  >  NC_000913/3766180‑3766435

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