Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	11,040	T→A	6.0%	L500Q (CTG→CAG)	dnaA →	chromosomal replication initiator protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	11,040	0	T	A	6.0%	209.9 / 2.9	67	L500Q (CTG→CAG)	dnaA	chromosomal replication initiator protein
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base A (2/2);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

AAGGAGTTGACCAACCACAGTCTGCCGGAAATCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAGCGCCGTCAGACCTTGCCGG  >  NC_002947/10901‑11187                                                                                                                                             |                                                                                                                                                    aaGGAGTTGACCAACCACAGTCTGCCGGAAATCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGAc                                                                                                                                            >  4:362567/1‑149 (MQ=255)       ttGACCAACCACAGTCTGCCGGAAATCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTgac                                                                                                                                      >  8:69274/1‑149 (MQ=255)        tGACCAACCACAGTCTGCCGGAAATCGGCGACATGTTCGGTGGTCGCGACCATACGAGCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTgac                                                                                                                                      <  7:69274/148‑1 (MQ=255)          accaaccaCAGTCTGCCGGAAATCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTgacgac                                                                                                                                   <  5:313459/149‑1 (MQ=255)           ccaaccaCAGTCTGCCGGAAATCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGAcc                                                                                                                                  >  5:92605/1‑149 (MQ=255)              accaCAGTCTGCCGGAAATCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTg                                                                                                                                        <  3:264298/140‑1 (MQ=255)              accaCAGTCTGCCGGAAATCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTg                                                                                                                                        >  4:264298/1‑140 (MQ=255)                           ggAAATCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTgcg                                                                                                                   >  3:108767/1‑148 (MQ=255)                              aaTCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGc                                                                                                               >  4:353987/1‑149 (MQ=255)                              aaTCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGc                                                                                                               >  6:27179/1‑149 (MQ=255)                                      aCATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGaa                                                                                                       <  3:362567/149‑1 (MQ=255)                                                  gtCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTa                                                                                           <  3:353987/149‑1 (MQ=255)                                                  gtCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTa                                                                                           <  5:27179/149‑1 (MQ=255)                                                     gcgACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGAc                                                                                        >  3:164875/1‑149 (MQ=255)                                                                aCCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCa                                                                             <  6:92605/149‑1 (MQ=255)                                                                      ctgcACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCAtt                                                                        >  4:367474/1‑148 (MQ=255)                                                                         cACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAAc                                                                    >  4:309906/1‑149 (MQ=255)                                                                          aCGCCTGCCGCAAAATCAATGAACTTAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAAcg                                                                   <  7:193260/149‑1 (MQ=255)                                                                               tGCCGAAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAg                                                              >  5:291703/1‑149 (MQ=255)                                                                                   gCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCaa                                                                     >  7:110830/1‑138 (MQ=255)                                                                                   gCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCaa                                                                     <  8:110830/138‑1 (MQ=255)                                                                                   gCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCt                                                          <  4:108767/149‑1 (MQ=255)                                                                                     aaaaTCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGt                                                        >  1:9467/1‑149 (MQ=255)                                                                                         tCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTAACCATTCAACGCGaa                                                               <  6:241116/138‑1 (MQ=255)                                                                                         tCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTAACCATTCAACGCGaa                                                               >  5:241116/1‑138 (MQ=255)                                                                                          cAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACCATGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGa                                                     >  1:308405/1‑148 (MQ=255)                                                                                            aTGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAAcc                                                 >  8:41371/1‑149 (MQ=255)                                                                                            aTGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCAGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAAca                                                 >  2:299864/1‑148 (MQ=255)                                                                                             tGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAAccc                                                >  3:178693/1‑149 (MQ=255)                                                                                               aaCTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAAcccc                                               <  2:9467/148‑1 (MQ=255)                                                                                                aCTGAATGAATCCGACGTTTACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCt                                              <  6:291703/148‑1 (MQ=255)                                                                                                aCTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTGGAAAccc                                                <  6:64850/146‑1 (MQ=255)                                                                                                aCTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTCGAAAccc                                                >  5:64850/1‑146 (MQ=255)                                                                                                  tGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGc                                            >  7:225200/1‑148 (MQ=255)                                                                                                   gAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAAcg                                                                   >  4:295808/1‑124 (MQ=255)                                                                                                   gAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAAcg                                                                   <  3:295808/124‑1 (MQ=255)                                                                                                    aaGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCa                                           <  4:178693/147‑1 (MQ=255)                                                                                                    aaGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAAc                                         >  8:162956/1‑149 (MQ=255)                                                                                                     aGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGaaa                                                   >  6:4867/1‑138 (MQ=255)                                                                                                     aGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGaaa                                                   <  5:4867/138‑1 (MQ=255)                                                                                                      ggAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTg                                       <  3:309906/149‑1 (MQ=255)                                                                                                          tCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCg                                   >  6:28349/1‑149 (MQ=255)                                                                                                             gACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTCCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCg                                >  2:51498/1‑149 (MQ=255)                                                                                                                gcgGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGgtg                             <  7:123857/149‑1 (MQ=255)                                                                                                                gcgGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTCGTCGCCGgtg                             >  8:123857/1‑149 (MQ=255)                                                                                                                     cATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTgtcgtc                        <  7:162956/149‑1 (MQ=255)                                                                                                                      aTCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGAcc                                                                                       >  8:374584/1‑85 (MQ=255)                                                                                                                      aTCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGAcc                                                                                       <  7:374584/85‑1 (MQ=255)                                                                                                                        ccGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCg                                   <  1:202201/135‑1 (MQ=255)                                                                                                                        ccGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCg                                   >  2:202201/1‑135 (MQ=255)                                                                                                                        ccGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAg                     >  3:200106/1‑149 (MQ=255)                                                                                                                         cgcgAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAg                     <  3:367474/148‑1 (MQ=255)                                                                                                                              ggACTACAAGAACC‑AGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGaa                                                               <  8:251442/101‑1 (MQ=255)                                                                                                                              ggACTACAAGAACC‑AGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGaa                                                               <  8:257077/101‑1 (MQ=255)                                                                                                                              ggACTACAAGAACC‑AGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGaa                                                               >  7:251442/1‑101 (MQ=255)                                                                                                                              ggACTACAAGAACC‑AGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGaa                                                               >  7:257077/1‑101 (MQ=255)                                                                                                                                  tACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAGCGCCGTCAGa           >  8:92017/1‑149 (MQ=255)                                                                                                                                  tACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAGCGCCGTCAGa           >  8:352331/1‑149 (MQ=255)                                                                                                                                    cAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCa                                                                             >  5:217548/1‑81 (MQ=255)                                                                                                                                    cAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCa                                                                             <  6:217548/81‑1 (MQ=255)                                                                                                                                    cAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGt                                                        >  4:359796/1‑102 (MQ=255)                                                                                                                                    cAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAGCGCCGTCAGAc          >  6:343422/1‑148 (MQ=255)                                                                                                                                    cAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAGCGCCGTCAGAc          <  5:343422/148‑1 (MQ=255)                                                                                                                                     aaGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCt                                              <  4:262349/111‑1 (MQ=255)                                                                                                                                     aaGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCt                                              >  3:262349/1‑111 (MQ=255)                                                                                                                                     aaGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAGCGCCGTCAGACCt        >  1:84384/1‑149 (MQ=255)                                                                                                                                       gAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAGCGCCGTCAGACCTTg      >  5:285284/1‑149 (MQ=255)                                                                                                                                           c‑tgctgCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAGCGCCGTCAGACCTTGCCgg  <  5:28349/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                    AAGGAGTTGACCAACCACAGTCTGCCGGAAATCGGCGACATGTTCGGTGGTCGCGACCATACGACCGTGCTGCACGCCTGCCGCAAAATCAATGAACTGAAGGAATCCGACGCGGACATCCGCGAGGACTACAAGAACC‑TGCTGCGGACGCTGACGACCTGATGGCCGTCTGCGCAGCTAATTGAAGGCAAGGGACTAGACCATGCATTTCACCATTCAACGCGAAGCCCTGTTGAAACCCCTGCAACTGGTCGCCGGTGTCGTCGAGCGCCGTCAGACCTTGCCGG  >  NC_002947/10901‑11187

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