Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,987,693	G→A	I22I (ATC→ATT)	yecJ ←	DUF2766 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,987,693	0	G	A	100.0%	242.6 / NA	63	I22I (ATC→ATT)	yecJ	DUF2766 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/26);  total (37/26)

AAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATC  >  NC_000913/1987546‑1987837                                                                                                                                                    |                                                                                                                                                 aaGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGAATGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAAt                                                                                                                                                >  1:183778/1‑150 (MQ=255)  aGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTtaataa                                                                                                                                              <  1:410788/151‑1 (MQ=255)    cgcgATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAAcc                                                                                                                                            >  2:307571/1‑151 (MQ=255)              ttttGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGgca                                                                                                                                    <  1:206476/149‑1 (MQ=255)               tttGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGgcaagc                                                                                                                                 >  1:876759/1‑151 (MQ=255)                ttGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACCTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGgcaagc                                                                                                                                 >  2:665541/1‑150 (MQ=255)                 tGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGgcaagca                                                                                                                                <  1:117982/150‑1 (MQ=255)                   gATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAAcc                                                                                                                             <  2:1227378/151‑1 (MQ=255)                   gATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAAcc                                                                                                                             >  1:665693/1‑151 (MQ=255)                    aTTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACca                                                                                                                            <  2:1100527/151‑1 (MQ=255)                       gCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACAt                                                                                                                         >  2:1161535/1‑151 (MQ=255)                        cccGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATc                                                                                                                        >  2:514555/1‑151 (MQ=255)                        cccGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATc                                                                                                                        >  2:715299/1‑151 (MQ=255)                         ccGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCa                                                                                                                       <  1:981893/151‑1 (MQ=255)                          cGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAg                                                                                                                      <  2:813543/151‑1 (MQ=255)                           gCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGa                                                                                                                     <  2:591654/151‑1 (MQ=255)                             aTGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGaaa                                                                                                                   <  2:904122/151‑1 (MQ=255)                                gTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCa                                                                                                                <  1:656769/151‑1 (MQ=255)                                    ccGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTc                                                                                                            >  2:1235571/1‑151 (MQ=255)                                      gggTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCt                                                                                                          <  2:146543/151‑1 (MQ=255)                                        gTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCt                                                                                                          >  2:723212/1‑149 (MQ=255)                                              ggcggcggGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGc                                                                                                   >  1:175960/1‑150 (MQ=255)                                                       aTGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGc                                                                                         <  2:75260/151‑1 (MQ=255)                                                       aTGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGc                                                                                         >  2:529991/1‑151 (MQ=255)                                                             aGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGAc                                                                                   <  2:1195502/151‑1 (MQ=255)                                                             aGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGAc                                                                                   <  1:299688/151‑1 (MQ=255)                                                                 aaGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCg                                                                               >  1:1187393/1‑151 (MQ=255)                                                                 aaGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCg                                                                               >  2:221497/1‑151 (MQ=255)                                                                     cGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTaaa                                                                           <  1:723184/151‑1 (MQ=255)                                                                       atatTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAAc                                                                          >  1:72215/1‑150 (MQ=255)                                                                              tcAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTg                                                                     <  1:332615/148‑1 (MQ=255)                                                                               cAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATa                                                                  <  1:794508/150‑1 (MQ=255)                                                                                aGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGc                                                                >  2:317083/1‑151 (MQ=255)                                                                                    gACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCAttt                                                            >  1:85383/1‑151 (MQ=255)                                                                                     aCTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCAtttt                                                           >  2:280715/1‑151 (MQ=255)                                                                                     aCTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCAtttt                                                           >  1:652273/1‑151 (MQ=255)                                                                                      cTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTg                                                          >  1:1152626/1‑151 (MQ=255)                                                                                         gACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGaa                                                        >  2:842571/1‑150 (MQ=255)                                                                                            aTTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGcc                                                    <  1:508182/151‑1 (MQ=255)                                                                                              tttttCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGt                                                  >  1:1026475/1‑151 (MQ=255)                                                                                                 ttCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTta                                                <  2:518872/150‑1 (MQ=255)                                                                                                      gAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTAtagt                                            >  1:744990/1‑149 (MQ=255)                                                                                                         cTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCg                                        >  1:930045/1‑150 (MQ=255)                                                                                                          ttCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGAc                                      >  1:182019/1‑151 (MQ=255)                                                                                                            cAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACtt                                    <  1:1102021/151‑1 (MQ=255)                                                                                                             aaCGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACtt                                    >  2:177463/1‑150 (MQ=255)                                                                                                              aCGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTca                                  <  2:594660/151‑1 (MQ=255)                                                                                                                    gcgATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCtt                            >  2:546123/1‑151 (MQ=255)                                                                                                                     cgATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTc                           <  1:961617/151‑1 (MQ=255)                                                                                                                      gATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTc                           >  1:1178688/1‑150 (MQ=255)                                                                                                                      gATTACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCa                          <  1:790219/151‑1 (MQ=255)                                                                                                                        ttACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCa                          >  1:784076/1‑149 (MQ=255)                                                                                                                        ttACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGc                        >  1:845875/1‑151 (MQ=255)                                                                                                                         tACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGc                        <  1:1169029/150‑1 (MQ=255)                                                                                                                         tACATCGAGAACATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCg                       >  1:536797/1‑151 (MQ=255)                                                                                                                                    cATCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACAt             >  1:433283/1‑150 (MQ=255)                                                                                                                                     aTCCAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATcc           >  2:1184591/1‑151 (MQ=255)                                                                                                                                       ccAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCa          <  1:699986/150‑1 (MQ=255)                                                                                                                                       ccAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCAc         >  2:423731/1‑151 (MQ=255)                                                                                                                                        cAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCAcc        >  2:897605/1‑151 (MQ=255)                                                                                                                                        cAGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCAcc        <  1:671580/151‑1 (MQ=255)                                                                                                                                         aGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCg       <  1:1235507/151‑1 (MQ=255)                                                                                                                                         aGAATTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCg       >  1:707781/1‑151 (MQ=255)                                                                                                                                            aTTTGTTTAATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGtc     >  2:1065732/1‑150 (MQ=255)                                                                                                                                              ttgtttaATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGtcat   <  2:869177/150‑1 (MQ=255)                                                                                                                                               tgtttaATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGtcatc  <  1:555043/150‑1 (MQ=255)                                                                                                                                                gtttaATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGtcatc  >  2:258539/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 AAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATC  >  NC_000913/1987546‑1987837

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