Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,575,200	T→A	9.2%	intergenic (‑186/+24)	pcaP ← / ← galP‑II	porin/porin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,575,200	0	T	A	9.2%	219.6 / 11.0	76	intergenic (‑186/+24)	pcaP/galP‑II	porin/porin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base A (3/4);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGAC  >  NC_002947/1575054‑1575344                                                                                                                                                   |                                                                                                                                                 gATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTa                                                                                                                                                >  8:244954/1‑149 (MQ=255)         tCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGaa                                                                                                                                        >  7:81714/1‑149 (MQ=255)           aGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGaa                                                                                                                                        >  6:263786/1‑147 (MQ=255)            ggATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAatag                                                                                                                                     <  8:81714/149‑1 (MQ=255)             gATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAgg                                                                                                                                    <  5:263786/149‑1 (MQ=255)              aTAGCCTGCAGTCACCACCTTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAATAGGc                                                                                                                                   <  4:159112/149‑1 (MQ=255)                  ccTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGACTATTCGCTGGTTTATCCGCGAAGAGGCCAAt                                                                                                                               >  5:345741/1‑149 (MQ=255)                   cTGCAGTCACCTCCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGTCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAATAGGCCaga                                                                                                                               <  7:329875/148‑3 (MQ=38)                     gCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAAtt                                                                                                                              >  2:175748/1‑147 (MQ=255)                      cAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAg                                                                                                                           <  3:187237/149‑1 (MQ=255)                                       agctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCg                                                                                                          <  7:244954/149‑1 (MQ=255)                                        gctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGATGGCGTCTGGCCTATTCGCGGGTTTCTCCGCGAAGAGGCCAATTGAGTCAAAGAAGATGCAGCg                                                                                                          >  8:297023/1‑148 (MQ=255)                                         ctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGc                                                                                                                                   <  2:281609/122‑1 (MQ=255)                                         ctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGc                                                                                                                                   >  1:281609/1‑122 (MQ=255)                                          tcagcGAATGGCGGCAGCTCGCACTTGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCCTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                            <  6:127575/128‑1 (MQ=255)                                            agcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATaa                                                                                                     <  2:338284/149‑1 (MQ=255)                                            agcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATa                                                                                                      >  8:347415/1‑148 (MQ=255)                                             gcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAAc                                                                                                    <  1:175748/149‑1 (MQ=255)                                              cGAATGGCGGCAGCTCGCACTGGTCAATTAATGTCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAATAGGGAGGGCGTCTTGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACt                                                                                                   <  6:345741/149‑1 (MQ=255)                                              cGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAAc                                                                                                    >  8:202076/1‑148 (MQ=255)                                                    gcggcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACaa                                                                                              <  3:136924/148‑1 (MQ=255)                                                       gcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAg                                                                                                            >  6:232319/1‑131 (MQ=255)                                                       gcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAg                                                                                                            <  5:232319/131‑1 (MQ=255)                                                       gcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCaa                                                                                          <  3:188096/149‑1 (MQ=255)                                                          gCTCGCACTGGTCAATTCATTGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCTAAAGGGAGGGCGTCTGGCCTATTCGCTGGTTTATGCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCg                                                                                       <  7:309749/149‑1 (MQ=255)                                                            tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCa                                                                                                                                 >  6:143972/1‑105 (MQ=255)                                                            tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCa                                                                                                                                 <  5:143972/105‑1 (MQ=255)                                                            tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGt                                                                                     >  5:360044/1‑149 (MQ=255)                                                            gcGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGt                                                                                     <  7:297023/148‑1 (MQ=255)                                                               cACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTtct                                                                                  <  8:49491/149‑1 (MQ=255)                                                                 cTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCa                                                                                <  7:347415/149‑1 (MQ=255)                                                                  tGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCAt                                                                               >  6:313828/1‑149 (MQ=255)                                                                  tGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCAt                                                                               >  8:40233/1‑149 (MQ=255)                                                                       aaTTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAAc                                                                          >  4:138098/1‑149 (MQ=255)                                                                          tCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAAcgcg                                                                       >  7:41044/1‑149 (MQ=255)                                                                             tGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                                    >  7:23274/1‑149 (MQ=255)                                                                             tGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                                    <  8:23274/149‑1 (MQ=255)                                                                               gCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAAGTCACAATCAAGCGGTTCTCATCGACCGCGTTgtt                                                                  >  1:175577/1‑147 (MQ=255)                                                                                 aGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAg                                                                                                            >  8:273543/1‑105 (MQ=255)                                                                                 aGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAg                                                                                                            <  7:273543/105‑1 (MQ=255)                                                                                 aGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                           <  4:254337/106‑1 (MQ=255)                                                                                 aGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                           >  3:254337/1‑106 (MQ=255)                                                                                 aGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGc                                                                >  8:297939/1‑149 (MQ=255)                                                                                             aCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGa                                                     <  6:360044/148‑1 (MQ=255)                                                                                               gCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGc                                                  <  5:313828/149‑1 (MQ=255)                                                                                                  gTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATg                                               >  5:57192/1‑149 (MQ=255)                                                                                                   ttGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATg                                               <  1:357015/148‑1 (MQ=255)                                                                                                     gCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTc                                            >  4:268283/1‑149 (MQ=255)                                                                                                      ccAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTc                                            >  5:67317/1‑148 (MQ=255)                                                                                                       cAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                                    <  6:373100/123‑1 (MQ=255)                                                                                                       cAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                                    >  5:373100/1‑123 (MQ=255)                                                                                                              gCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgcgc                                   >  3:326563/1‑149 (MQ=255)                                                                                                                cTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCa                                 >  4:239512/1‑149 (MQ=255)                                                                                                                 ttCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                                >  8:329122/1‑149 (MQ=255)                                                                                                                   cAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                                <  7:202076/147‑1 (MQ=255)                                                                                                                   cAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                                <  3:268283/147‑1 (MQ=255)                                                                                                                         ggagggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgc                                     <  7:308702/136‑1 (MQ=255)                                                                                                                         ggagggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgc                                     >  8:308702/1‑136 (MQ=255)                                                                                                                                gTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTcc                 >  6:275831/1‑149 (MQ=255)                                                                                                                                  ttGGCCTCTTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCtt               >  8:208694/2‑149 (MQ=39)                                                                                                                                  cTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCtt               >  2:144201/1‑149 (MQ=255)                                                                                                                                    ggCCTCTTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACATGTCCtt               <  3:21013/147‑1 (MQ=255)                                                                                                                                    ggCCTCTTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCtt               >  4:21013/1‑147 (MQ=255)                                                                                                                                     gCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCtt               <  7:40233/146‑1 (MQ=255)                                                                                                                                     gCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGa             <  7:297939/148‑1 (MQ=255)                                                                                                                                     gCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGa             <  6:67317/148‑1 (MQ=255)                                                                                                                                       cTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTAAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGATGAATTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGc          >  8:181237/1‑149 (MQ=255)                                                                                                                                        tATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCt         >  6:289849/1‑149 (MQ=255)                                                                                                                                             gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCt                <  8:271334/137‑1 (MQ=255)                                                                                                                                             gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCt                >  7:271334/1‑137 (MQ=255)                                                                                                                                             gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGa             <  5:209729/140‑1 (MQ=255)                                                                                                                                             gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGa             >  6:209729/1‑140 (MQ=255)                                                                                                                                             gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGccc     >  8:285822/1‑148 (MQ=255)                                                                                                                                             gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGccc     >  2:232880/1‑148 (MQ=255)                                                                                                                                              cgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGa   >  8:248312/1‑149 (MQ=255)                                                                                                                                               gggTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGAc  >  2:204625/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGAC  >  NC_002947/1575054‑1575344

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: ‑