Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,575,203	T→C	9.1%	intergenic (‑189/+21)	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,203	0	T	C	9.1%	214.6 / 9.9	77	intergenic (‑189/+21)	pcaP/galP‑II	porin/porin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (37/33);  new base C (3/4);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGCAC  >  NC_002947/1575062‑1575349                                                                                                                                              |                                                                                                                                                   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                                                                                                                                      <  5:143972/105‑1 (MQ=255)                                                    tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCa                                                                                                                                      >  6:143972/1‑105 (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                                                                         >  5:373100/1‑123 (MQ=255)                                                                                               cAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                                         <  6:373100/123‑1 (MQ=255)                                                                                                      gCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgcgc                                        >  3:326563/1‑149 (MQ=255)                                                                                                        cTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCa                                      >  4:239512/1‑149 (MQ=255)                                                                                                         ttCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                                     >  8:329122/1‑149 (MQ=255)                                                                                                           cAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                                     <  3:268283/147‑1 (MQ=255)                                                                                                           cAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                                     <  7:202076/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                  >  6:209729/1‑140 (MQ=255)                                                                                                                                     gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGa                  <  5:209729/140‑1 (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)                                                                                                                                            tATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGca   <  1:232880/148‑1 (MQ=255)                                                                                                                                             ttCCGCGAAGAGGCCAATTTTGTCACAGAAGATTCAGCTGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGGATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGcac  <  6:57192/147‑1 (MQ=255)                                                                                                                                             aTCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGcac  <  7:285822/148‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   TCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGCAC  >  NC_002947/1575062‑1575349

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