Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,425	T→C	I294I (ATT→ATC)	ssrA →	chromosomal replication initiator protein DnaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,425	0	T	C	100.0%	235.2 / NA	90	I294I (ATT→ATC)	ssrA	chromosomal replication initiator protein DnaA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (82/8);  total (82/8)

TAGAAGACCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATTGAAGAAGAAAAATTAGATATTCCACCAGAAGCTTTAAATTATATAGCAAATCAAATTCAATCTAATATTCGTGAATTAGAAGGTGCA  >  NC_007793/1328‑1512                                                                                                  |                                                                                        tAGAAGACCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCga                                                                                       >  1:3527748/1‑100 (MQ=255) tAGAAGACCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCga                                                                                       >  1:2940894/1‑100 (MQ=255)  agaagaCCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCgaa                                                                                      >  1:2718115/1‑100 (MQ=255)   caagaCCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCgaag                                                                                     >  1:3076417/2‑100 (MQ=255)    aagaCCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCgaaga                                                                                    <  1:4079900/100‑1 (MQ=255)     agaCCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCgaagaa                                                                                   >  1:1569614/1‑100 (MQ=255)     agaCCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCgaagaa                                                                                   >  1:5648724/1‑100 (MQ=255)      gaCCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGACTGGCAATTTTGCAGAAGAAAATCgaagaag                                                                                  >  1:5406248/1‑100 (MQ=255)      gaCCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCgaagaag                                                                                  >  1:339893/1‑100 (MQ=255)      acccGATTACGTTCACGCTTTGAATGGGGGTTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCgaagaag                                                                                  >  1:1705342/3‑100 (MQ=255)           aTTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGaaaaa                                                                             >  1:3941930/1‑100 (MQ=255)            ttACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAAt                                                                            >  1:1325723/1‑100 (MQ=255)             tACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAAtt                                                                           >  1:2984302/1‑100 (MQ=255)              aCGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCAAAGAAGAAAAATTa                                                                          >  1:4885100/1‑100 (MQ=255)               ccctCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAg                                                                         >  1:4081158/4‑100 (MQ=255)               cGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAg                                                                         >  1:1526980/1‑100 (MQ=255)                tgtCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGa                                                                        >  1:3879645/3‑100 (MQ=255)                gcccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGa                                                                        >  1:5310459/4‑100 (MQ=255)                atccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGa                                                                        >  1:2990658/4‑100 (MQ=255)                 tgcACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGat                                                                       >  1:6410511/3‑100 (MQ=255)                 ctCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTTTGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGat                                                                       >  1:1134783/2‑100 (MQ=255)                 ccgaCGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAAAAGAAAAATTAGat                                                                       >  1:4965239/4‑100 (MQ=255)                 ccgaCGCTTTGAATGAGGGCTAGTTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGat                                                                       >  1:3207193/4‑100 (MQ=255)                 cccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAATTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGat                                                                       >  1:2050008/3‑100 (MQ=255)                 cccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAGTGGCAATTTTGCAGAAGAAAATCGAAGAAGaaaagatcgga                                                                       >  1:497359/3‑93 (MQ=255)                 cccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGat                                                                       >  1:3694847/3‑100 (MQ=255)                 cccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGat                                                                       >  1:2746903/3‑100 (MQ=255)                 cccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGat                                                                       >  1:307931/3‑100 (MQ=255)                 acgaCGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGat                                                                       >  1:2669272/4‑100 (MQ=255)                 acgaCGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGat                                                                       >  1:1846766/4‑100 (MQ=255)                  tctcGCTTTGAATGGGGGCTAATTATAGATATTGCGCCACCAGAGTATGAAACTCGAATGGTGATTTTGCAGACGAAAACCGAAGAAGAAAAATTAGata                                                                      >  1:964714/4‑100 (MQ=255)                  tCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGGAGAAAAATTAGata                                                                      >  1:4797294/1‑100 (MQ=255)                  tCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:3091376/1‑100 (MQ=255)                  tCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:3138572/1‑100 (MQ=255)                  tCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:4063477/1‑100 (MQ=255)                  tCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:6009933/1‑100 (MQ=255)                  tCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:1953550/1‑100 (MQ=255)                  tCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:3600331/1‑100 (MQ=255)                  tCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:5518927/1‑100 (MQ=255)                  gcACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:4691227/2‑100 (MQ=255)                  gcACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:3963387/2‑100 (MQ=255)                  ccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGTAATTTTGCAGAAGAAAATCGAAGAAG‑AAAATTAGatat                                                                     >  1:910445/2‑100 (MQ=255)                  ccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGCAAAATTAGata                                                                      >  1:561341/2‑100 (MQ=255)                  ccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:2583396/2‑100 (MQ=255)                  ccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:5749974/2‑100 (MQ=255)                  ccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:2775592/2‑100 (MQ=255)                  ccACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAACTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:5821546/2‑100 (MQ=255)                  acACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGTAGAAAAATTAGata                                                                      >  1:1848057/2‑100 (MQ=255)                  acACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGata                                                                      >  1:2862658/2‑100 (MQ=255)                   catgCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAGGCAAAATTAGatct                                                                     >  1:4953636/4‑98 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatgt                                                                     >  1:2097820/1‑98 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:6187741/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:6310902/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:712390/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:6057308/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:5649205/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:5976294/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:93831/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:3790281/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:1069426/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:2490759/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:4193976/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:2795305/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:2881942/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:4714239/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:3668596/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:339695/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:3237499/1‑100 (MQ=255)                   cACGCTTTGAATGGGGGCTAATTGATGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGTAGAAGAAAAATTAGatat                                                                     >  1:3597908/1‑100 (MQ=255)                   cACGCTTTGAATGGGAGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGACAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:2793715/1‑100 (MQ=255)                   accGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGatat                                                                     >  1:4265583/3‑100 (MQ=255)                    ccGCTTTGAATGGGGGCTAATTGTGGAGATAACGCCACCAGATTATGTAACTCGAATGGCAAATTTGCAGAAAAAAATCGAAGAAGAAAAAGTAGATAtt                                                                    >  1:5251203/2‑100 (MQ=255)                    aCGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATAtt                                                                    >  1:3879524/1‑100 (MQ=255)                     cGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTc                                                                   >  1:2994586/1‑100 (MQ=255)                         ttccctGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTtcacc                                                               <  1:2464884/95‑1 (MQ=255)                         ttGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTccacc                                                               <  1:101853/100‑1 (MQ=255)                         ttGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATGccacc                                                               <  1:6361966/100‑1 (MQ=255)                         ttGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAAAAGAAAAATTAGATATTccacc                                                               <  1:2401796/100‑1 (MQ=255)                          tGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTccacca                                                              >  1:6461220/1‑100 (MQ=255)                            aaTGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGa                                                            <  1:277987/100‑1 (MQ=255)                            aaTGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGa                                                            <  1:6368993/100‑1 (MQ=255)                            aaTGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAACGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGa                                                            >  1:1710423/1‑100 (MQ=255)                                                atTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGAAGCTTTAAATTATATAGcaa                                        >  1:5450632/1‑100 (MQ=255)                                                   aCGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGAAGCTTTAAATTATATAGcaaatc                                     >  1:5441384/1‑100 (MQ=255)                                                        accaGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAATATTAGATATTCCACCAGAAGCTTTAAATTATAGAGCAAATCAAAtt                                >  1:2701183/1‑100 (MQ=255)                                                        accaGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGAAGCTTTAAATTATGTAGCAAATCAAAtt                                >  1:4905235/1‑100 (MQ=255)                                                        accaGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGAAGCTTTAAATTATATAGcaagatcggaa                                >  1:6006795/1‑92 (MQ=255)                                                         ccaGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGAAGCTTTAAATTATATAGcaagatcggaag                               >  1:1682627/1‑91 (MQ=255)                                                               gaTGAAACTCGAATGGCAATTTTGCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGAAGCTTTAAATTATATAGCAAATCAAATTCAATCTa                         >  1:3872951/2‑100 (MQ=255)                                                                                      gCAGAAGAAAATCGAAGAAGAAAAATTAGATATTCCACCAGAAGCTTTAAATTATATAGCAAATCAAATTCAATCTAATATTCGTGAATTAGAAGGTGCa  <  1:4599869/100‑1 (MQ=255)                                                                                                  |                                                                                        TAGAAGACCGATTACGTTCACGCTTTGAATGGGGGCTAATTGTTGATATTACGCCACCAGATTATGAAACTCGAATGGCAATTTTGCAGAAGAAAATTGAAGAAGAAAAATTAGATATTCCACCAGAAGCTTTAAATTATATAGCAAATCAAATTCAATCTAATATTCGTGAATTAGAAGGTGCA  >  NC_007793/1328‑1512

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