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%	259.0 / NA	80	*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 (38/42);  total (38/42)

ATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGA  >  NC_000913/3766180‑3766434                                                                                                                                         |                                                                                                                       aTGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAAcc                                                                                                                      >  2:1349313/1‑139 (MQ=255) aTGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAAcc                                                                                                                      >  1:1339882/1‑139 (MQ=255)       aaaaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtgtg                                                                                                                >  1:1416959/1‑139 (MQ=255)        aaaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTgg                                                                                                               <  1:1349313/139‑1 (MQ=255)            cccGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGt                                                                                                           <  1:114468/139‑1 (MQ=255)                     tgtgACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAAt                                                                                                  >  2:733670/1‑139 (MQ=255)                        gACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                               <  2:1416959/139‑1 (MQ=255)                           gTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCaaa                                                                                            <  1:657275/139‑1 (MQ=255)                            tACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAat                                                                                           >  2:806843/1‑139 (MQ=255)                              cTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAatac                                                                                         >  2:1174857/1‑139 (MQ=255)                              cTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGCCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAatac                                                                                         >  1:1559524/1‑139 (MQ=255)                               ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                               >  2:605942/1‑132 (MQ=255)                               ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                               <  1:605942/132‑1 (MQ=255)                                  aGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAatacatac                                                                                     <  1:273120/139‑1 (MQ=255)                                   ggAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATAcc                                                                                    >  2:570669/1‑139 (MQ=255)                                         aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaa                                                                              >  2:1478671/1‑139 (MQ=255)                                               ttgTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAgg                                                                        <  1:570669/139‑1 (MQ=255)                                                tgtgGTGATATCAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                       <  1:1097470/139‑1 (MQ=255)                                                    gtgATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                                                   >  1:1610808/1‑139 (MQ=255)                                                       atatTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                                                <  1:1478671/139‑1 (MQ=255)                                                           tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            <  2:1488194/139‑1 (MQ=255)                                                           tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            <  2:1559524/139‑1 (MQ=255)                                                           tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            <  1:733670/139‑1 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATGAAAAAAGATAACCCTTGTGGAGGt                                                                                                           <  2:1146673/91‑1 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGt                                                                                                           >  1:1146673/1‑91 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTa                                                           >  1:33333/1‑139 (MQ=255)                                                               gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATgg                                                        <  1:806843/139‑1 (MQ=255)                                                                cTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGt                                                       <  2:1610808/139‑1 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                    <  1:98793/91‑1 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                    >  2:98793/1‑91 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                               <  2:296840/96‑1 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                               >  1:296840/1‑96 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTa                                                                  <  2:1511143/125‑1 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTa                                                                  >  1:1511143/1‑125 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCtt                                                                                                                   <  2:1651166/71‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCtt                                                                                                                   >  1:1651166/1‑71 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                                              <  1:107382/76‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                                              >  2:107382/1‑76 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGt                                                                                                           >  2:1517780/1‑79 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGt                                                                                                           <  1:1517780/79‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTg                                                                                                      <  2:1248757/84‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTg                                                                                                      >  1:1248757/1‑84 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            >  1:341160/1‑126 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            <  2:341160/126‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATaa                                               >  2:724681/1‑139 (MQ=255)                                                                         cAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtg                                                                                                                  <  2:1426417/71‑1 (MQ=255)                                                                         cAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtg                                                                                                                  >  1:1426417/1‑71 (MQ=255)                                                                         cAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                                              <  1:724681/139‑1 (MQ=255)                                                                           aaaaGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCa                                            <  2:33333/139‑1 (MQ=255)                                                                                 acacaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGt                                                                                                           >  2:422969/1‑70 (MQ=255)                                                                                 acacaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGt                                                                                                           <  1:422969/70‑1 (MQ=255)                                                                                     aGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                                              >  1:903912/1‑127 (MQ=255)                                                                                     aGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                                              <  2:903912/127‑1 (MQ=255)                                                                                          gCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                                          >  1:869957/1‑94 (MQ=255)                                                                                          gCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                                          <  2:869957/94‑1 (MQ=255)                                                                                            cTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                    <  2:1635508/66‑1 (MQ=255)                                                                                            cTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                    >  1:1635508/1‑66 (MQ=255)                                                                                               gAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAatgatg                        <  2:1339882/139‑1 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAgg                                                                        <  1:430446/85‑1 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAgg                                                                        >  2:430446/1‑85 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAatgatg                        >  1:1425216/1‑133 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAatgatg                        <  2:1425216/133‑1 (MQ=255)                                                                                                          aCTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            <  2:1194767/92‑1 (MQ=255)                                                                                                          aCTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            >  1:1194767/1‑92 (MQ=255)                                                                                                            tCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTCCg           <  1:1174857/139‑1 (MQ=255)                                                                                                               aGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                                                   <  1:594353/80‑1 (MQ=255)                                                                                                               aGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                                                   >  2:594353/1‑80 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGt                                                                                                >  1:1537054/1‑45 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGt                                                                                                <  2:1537054/45‑1 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                                                   >  2:48469/1‑74 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                                                   <  1:48469/74‑1 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            >  2:417723/1‑81 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            >  1:791021/1‑81 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            <  2:791021/81‑1 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            <  1:417723/81‑1 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTtga  >  2:839992/1‑139 (MQ=255)                                                                                                                           gATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            >  1:1251580/1‑75 (MQ=255)                                                                                                                           gATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                            <  2:1251580/75‑1 (MQ=255)                                                                                                                             tAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGt                                                                                                >  2:2038/1‑37 (MQ=255)                                                                                                                             tAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGt                                                                                                <  1:2038/37‑1 (MQ=255)                                                                                                                                         |                                                                                                                       ATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGA  >  NC_000913/3766180‑3766434

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