Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,224,724	T→C	6.0%	E213G (GAG→GGG)	minD ←	inhibitor of FtsZ ring polymerization; chromosome‑membrane tethering protein; membrane ATPase of the MinCDEE system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,224,724	0	T	C	6.0%	210.2 / 4.1	67	E213G (GAG→GGG)	minD	inhibitor of FtsZ ring polymerization; chromosome‑membrane tethering protein; membrane ATPase of the MinCDEE system
Reads supporting (aligned to +/- strand):  ref base T (35/28);  new base C (2/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-02

CTTTCTTCTCTTCTTCAATGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCATTTTCTGCGCGGCGTGA  >  NC_000913/1224578‑1224872                                                                                                                                                   |                                                                                                                                                     cTTTCTTCTCTTCTTCAATGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCctct                                                                                                                                                    >  3:517831/1‑149 (MQ=255)   ttcttcTCTTCTTCAATGAAGCGGAAAGGACGTTCTTCTCCTAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTgg                                                                                                                                                  <  2:527054/149‑1 (MQ=255)   ttcttcTCTTCTTCAATGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTgg                                                                                                                                                  <  1:237513/149‑1 (MQ=255)   ttcttcTCTTCTTCAATGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTgg                                                                                                                                                  <  1:399789/149‑1 (MQ=255)            tcttcAATGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCCCTGGGATCACGcc                                                                                                                                         >  3:444264/1‑149 (MQ=255)            tcttcAATGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCCCTGGGATCACGcc                                                                                                                                         >  3:383735/1‑149 (MQ=255)            tcttcAATGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGCCGCGCAATACTGATTGATCCTCTGGGATCACGcc                                                                                                                                         >  3:519874/1‑149 (MQ=255)               tcAATGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCcgac                                                                                                                                      >  1:508540/1‑149 (MQ=255)                 aaTGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCcgacga                                                                                                                                    >  1:406894/1‑149 (MQ=255)                   tGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTTATTGATCCTCTGGGATCACGCCGACGAGt                                                                                                                                  <  4:489253/149‑1 (MQ=255)                    gAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGtt                                                                                                                                 >  1:416177/1‑149 (MQ=255)                    gAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGtt                                                                                                                                 <  3:43240/149‑1 (MQ=255)                    gAAGCGGAAAGGACGTTCTTCTCCAAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGTGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGtt                                                                                                                                 >  1:527776/1‑149 (MQ=255)                       gccgAAAGGACGTTGTTCTCCCTTCGGACGTTCTACGGTATCTGCGTAGGCTTTCCCTGGATAGGATTTAATGTCTTGAATGATCGGTTCTTCCTAGTTGGTTGCGCGCAATACTGATTTATCCTCTGGGATCACGCGGGCGAGTTTGa                                                                                                                              <  4:205074/146‑1 (MQ=255)                       gCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGa                                                                                                                              <  2:533558/149‑1 (MQ=255)                        cGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGAt                                                                                                                             <  2:508540/149‑1 (MQ=255)                          gAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATgc                                                                                                                           >  1:475164/1‑149 (MQ=255)                              ggACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAg                                                                                                                       <  3:108298/149‑1 (MQ=255)                               gACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAgg                                                                                                                      >  4:376677/1‑149 (MQ=255)                                aCGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGa                                                                                                                     >  3:457934/1‑149 (MQ=255)                                     cttctCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTcc                                                                                                                <  4:467415/149‑1 (MQ=255)                                       tctcCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAg                                                                                                              <  3:150783/149‑1 (MQ=255)                                       tctcCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAg                                                                                                              <  4:534114/149‑1 (MQ=255)                                         tcCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGca                                                                                                            >  4:22781/1‑149 (MQ=255)                                            cAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACAt                                                                                                         >  3:210695/1‑149 (MQ=255)                                            cAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACAt                                                                                                         >  1:439873/1‑149 (MQ=255)                                             aaCAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATc                                                                                                        >  1:419426/1‑149 (MQ=255)                                             aaCAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATc                                                                                                        >  1:90850/1‑149 (MQ=255)                                              aCAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCt                                                                                                       >  1:470166/1‑149 (MQ=255)                                                 gACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCcgacga                                                                                                                                    <  1:248948/117‑1 (MQ=255)                                                 gACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCcgacga                                                                                                                                    >  2:248948/1‑117 (MQ=255)                                                     ttCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGc                                                                                                >  1:300006/1‑149 (MQ=255)                                                       cTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTc                                                                                              <  2:225903/149‑1 (MQ=255)                                                       cTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTc                                                                                              <  2:44537/149‑1 (MQ=255)                                                          cGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGc                                                                                           <  2:470166/149‑1 (MQ=255)                                                             tATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATg                                                                                        >  3:444464/1‑149 (MQ=255)                                                                      aGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTg                                                                               >  3:181178/1‑149 (MQ=255)                                                                               ccGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCgg                                                                      >  1:480087/1‑149 (MQ=255)                                                                               ccGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCgg                                                                      >  3:112765/1‑149 (MQ=255)                                                                                cGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGCCGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGc                                                                     >  1:426610/1‑149 (MQ=255)                                                                                 gCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGcc                                                                    <  1:273912/149‑1 (MQ=255)                                                                                     cGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTggg                                                                <  4:489841/149‑1 (MQ=255)                                                                                       gCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGTTGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGtt                                                              >  3:535659/1‑149 (MQ=255)                                                                                         gTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTtat                                                            >  3:296596/1‑149 (MQ=255)                                                                                                tCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGt                                                     >  4:195066/1‑149 (MQ=255)                                                                                                    gaATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAAc                                                 >  3:493316/1‑149 (MQ=255)                                                                                                     aaTGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAAca                                                <  2:214616/149‑1 (MQ=255)                                                                                                        gACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAAcagca                                             >  3:483853/1‑149 (MQ=255)                                                                                                             gTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATTCCCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGc                                        <  4:383735/149‑1 (MQ=255)                                                                                                             gTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCCCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGc                                        <  4:444264/149‑1 (MQ=255)                                                                                                                  cccTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCttt                                   <  2:527776/149‑1 (MQ=255)                                                                                                                  cccTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCttt                                   <  2:416177/149‑1 (MQ=255)                                                                                                                   ccTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTaa                                                  >  2:162119/1‑133 (MQ=255)                                                                                                                   ccTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTaa                                                  <  1:162119/133‑1 (MQ=255)                                                                                                                   ccTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTa                                  <  1:537391/149‑1 (MQ=255)                                                                                                                   ccTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTa                                  >  2:100090/1‑149 (MQ=255)                                                                                                                      ggTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATa                               <  4:444464/149‑1 (MQ=255)                                                                                                                            agGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGctct                         >  4:185337/1‑149 (MQ=255)                                                                                                                            agGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGctct                         >  4:267234/1‑149 (MQ=255)                                                                                                                              gcgcgcAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTc                       >  1:541100/1‑149 (MQ=255)                                                                                                                              gcgcgcAATACTGATTGATCCTCTGGGATCACGCCGACGACTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTc                       <  4:493316/149‑1 (MQ=255)                                                                                                                                      tACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCAtttt               >  4:91603/1‑149 (MQ=255)                                                                                                                                        cTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCATTTTCt             <  2:409936/149‑1 (MQ=255)                                                                                                                                          gattgatCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCATTTTCTgc           <  4:457934/149‑1 (MQ=255)                                                                                                                                          gattgatCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCATTTTCTgc           <  1:436312/149‑1 (MQ=255)                                                                                                                                              gatCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCATTTTCTGCgcgg       >  4:478060/1‑149 (MQ=255)                                                                                                                                              gatCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCATTTTCTGCgcgg       >  2:320275/1‑149 (MQ=255)                                                                                                                                                 ccTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCATTTTCTGCGCGGCGt    >  1:265627/1‑149 (MQ=255)                                                                                                                                                  ctctGGGATCACGCCGACGAGTTTGATGCGCAGGATATCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATGGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCTATTTTCTGCGCGGCGTg   >  4:478657/1‑149 (MQ=255)                                                                                                                                                   tctGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCATTTTCTGCGCGGCGTGa  >  2:419528/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CTTTCTTCTCTTCTTCAATGAAGCGGAAAGGACGTTCTTCTCCCAACAGACGTTCTACGGTATCTGCGTAGGCTTTACCCGCATCGGCGTTAATGTCGAGAATGACCGGTTCACCCTGGTTAGAGGCGCGCAATACTGATTGATCCTCTGGGATCACGCCGACGAGTTTGATGCGCAGGATCTCCAGCACATCTTCCATGCTCAGCATGTCACCTCTGCTTACGCGGCCTGGGTTATAGCGCGTTAACAGCAGGTGCTCTTTAATAGGCTCTTCGCCATTTTCTGCGCGGCGTGA  >  NC_000913/1224578‑1224872

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