Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,815,591	A→G	P575P (CCA→CCG)	katE →	catalase HPII, heme d‑containing

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,815,591	0	A	G	100.0%	233.6 / NA	64	P575P (CCA→CCG)	katE	catalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/30);  total (34/30)

AGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTG  >  NC_000913/1815448‑1815735                                                                                                                                                |                                                                                                                                                 aGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGAATCGAACTGACTGCCGACCAGCTGAATATCACCCCGCCTCCgg                                                                                                                                           >  1:1061949/1‑151 (MQ=255)     ttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGCCTGACGACCAGCTGAATATCCCCCCGCCTCCGGACg                                                                                                                                        <  1:1126047/150‑1 (MQ=255)     ttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGt                                                                                                                                       >  2:1199871/1‑151 (MQ=255)      tGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTc                                                                                                                                      >  2:367469/1‑151 (MQ=255)       gAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCa                                                                                                                                     <  1:392966/151‑1 (MQ=255)           tAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACgg                                                                                                                                 >  1:972619/1‑151 (MQ=255)            aaGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACgg                                                                                                                                 <  1:691267/150‑1 (MQ=255)            aaGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACgg                                                                                                                                 >  1:1165646/1‑150 (MQ=255)            aaGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACgg                                                                                                                                 <  1:288790/150‑1 (MQ=255)               cAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTg                                                                                                                             <  1:421759/151‑1 (MQ=255)                 aaGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGaa                                                                                                                           <  1:67630/151‑1 (MQ=255)                 aaGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGaa                                                                                                                           <  2:63283/151‑1 (MQ=255)                    tcgttcgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGaaaa                                                                                                                         >  2:692748/1‑150 (MQ=255)                    tcgttcgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGaaaa                                                                                                                         >  1:860006/1‑150 (MQ=255)                        tcgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGat                                                                                                                    >  2:827755/1‑151 (MQ=255)                        tcgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGa                                                                                                                     <  1:593705/150‑1 (MQ=255)                         cgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGatc                                                                                                                   >  2:850593/1‑151 (MQ=255)                            ccGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGatcca                                                                                                                 <  2:452606/150‑1 (MQ=255)                             cGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGatccatc                                                                                                               >  1:119019/1‑151 (MQ=255)                                atatTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCtt                                                                                                            >  1:1045542/1‑151 (MQ=255)                                atatTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACGCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCtt                                                                                                            <  1:1172463/151‑1 (MQ=255)                                 tatTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCtt                                                                                                            <  1:705817/150‑1 (MQ=255)                                   ttCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTa                                                                                                           >  1:105085/1‑149 (MQ=255)                                     cGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGtt                                                                                                       >  2:33884/1‑151 (MQ=255)                                        gAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTg                                                                                                     <  1:1013310/150‑1 (MQ=255)                                        gAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGt                                                                                                    >  2:772797/1‑151 (MQ=255)                                         aGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGt                                                                                                    >  1:943578/1‑150 (MQ=255)                                           cgcgTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGCCAAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTa                                                                                                   <  1:621126/149‑1 (MQ=255)                                             cgTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGcc                                                                                               >  2:874976/1‑151 (MQ=255)                                                   tgACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTcc                                                                                          <  2:176917/150‑1 (MQ=255)                                                    gACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCt                                                                                         <  1:802325/150‑1 (MQ=255)                                                    gACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCt                                                                                         <  1:1104581/150‑1 (MQ=255)                                                    gACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTg                                                                                        >  2:115689/1‑151 (MQ=255)                                                      ccAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGAc                                                                                      <  2:1195169/151‑1 (MQ=255)                                                         gCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGt                                                                                   <  2:1024942/151‑1 (MQ=255)                                                                cATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTTGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGAtgtg                                                                             >  1:956215/1‑150 (MQ=255)                                                                cATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGAtgtg                                                                             <  2:102652/150‑1 (MQ=255)                                                                  tatTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGaa                                                                           <  2:709315/150‑1 (MQ=255)                                                                  tatTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGaa                                                                           >  1:782809/1‑150 (MQ=255)                                                                         ctcACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTACCGAACTGGCTGACGACCAGCTGAATATCCCCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTcg                                                                     <  2:559097/149‑1 (MQ=255)                                                                               ctGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGt                                                               >  2:260289/1‑149 (MQ=255)                                                                                 ggCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAg                                                             <  2:490297/149‑1 (MQ=255)                                                                                  gCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGa                                                          <  2:512807/151‑1 (MQ=255)                                                                                    ccAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAtt                                                        >  1:545526/1‑151 (MQ=255)                                                                                      aGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAtttt                                                      >  1:733229/1‑151 (MQ=255)                                                                                      aGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAtttt                                                      >  1:370326/1‑151 (MQ=255)                                                                                        gCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTAc                                                    >  1:781412/1‑151 (MQ=255)                                                                                          ggTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACtt                                                  >  2:1023198/1‑151 (MQ=255)                                                                                             ggCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAt                                                         <  1:722565/141‑1 (MQ=255)                                                                                             ggCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAt                                                         >  2:722575/1‑141 (MQ=255)                                                                                             ggCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAat                                               >  2:172863/1‑151 (MQ=255)                                                                                             ggCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAat                                               <  1:844962/151‑1 (MQ=255)                                                                                                        tcGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTg                                      >  2:949475/1‑149 (MQ=255)                                                                                                         cGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTgag                                    <  1:946115/150‑1 (MQ=255)                                                                                                             aTCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTg                                      >  1:461330/1‑144 (MQ=255)                                                                                                             aTCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCg                                <  1:1030131/150‑1 (MQ=255)                                                                                                                      actgacGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTc                      >  2:1136575/1‑151 (MQ=255)                                                                                                                          acgacCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTgg                   <  1:757279/150‑1 (MQ=255)                                                                                                                            gacCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGcc                 >  2:115402/1‑150 (MQ=255)                                                                                                                                    gAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCaa         <  2:980184/150‑1 (MQ=255)                                                                                                                                     aaTATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAgg       >  1:30281/1‑151 (MQ=255)                                                                                                                                      atatCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgc      >  1:741025/1‑151 (MQ=255)                                                                                                                                        atCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgcg     <  1:155903/150‑1 (MQ=255)                                                                                                                                           aCCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTg  >  2:347582/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                 AGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTG  >  NC_000913/1815448‑1815735

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