Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,740,471	0	A	T	26.2%	100.9 / 37.3	67	intergenic (+36/+35)	hutF/PP_5037	formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (19/29);  new base T (10/7);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGT  >  NC_002947/5740326‑5740617                                                                                                                                                  |                                                                                                                                                   tGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGcc                                                                                                                                                 >  6:413142/1‑149 (MQ=255)      gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCAcgcg                                                                                                                                             >  6:287884/1‑148 (MQ=39)        tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                          <  2:375299/149‑1 (MQ=39)          tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                                        <  3:324242/149‑1 (MQ=255)          tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                                        <  6:373585/149‑1 (MQ=255)          tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg                                                                                                                                         <  2:265670/148‑1 (MQ=255)            aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                      <  5:287884/149‑3 (MQ=38)            aTGGTGAATGGCCAATGGGTGGGGCGCCAGGGGCGGCATGCCCAGGAAGAGGAACGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTGGCGTGGCACCCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                      <  1:25304/149‑3 (MQ=255)             tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg                                                                                                                                         >  4:325395/1‑145 (MQ=255)             tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg                                                                                                                                         <  3:325395/145‑1 (MQ=255)             tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                     <  1:135330/149‑1 (MQ=255)             tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                     >  5:207412/1‑149 (MQ=255)                  aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc                                                                                                                                    <  6:7883/145‑1 (MQ=255)                  aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc                                                                                                                                    >  5:7883/1‑145 (MQ=255)                       tctttGGGTGGTGCGGCAGGGGCGGAATGCCAAGTAAGTTGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGAAACGTAAcc                                                                                                                           <  4:1422/145‑1 (MQ=255)                       gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGt                                                                                                                                       >  4:199848/1‑137 (MQ=255)                       gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGt                                                                                                                                       <  3:199848/137‑1 (MQ=255)                                 gtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCCTTTCAGCCCTTACCCTTGAAGCCCCCGAAAGGTTTCAACGCAACCCCAAAAAAc                                                                                                                  >  6:83079/1‑148 (MQ=255)                                  tgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTc                                                                                                                <  2:371033/149‑1 (MQ=255)                                     gcCAGGGGCGGCATGCCCAGGAAGATGAACGCGGGGGGGCGTTTGCCACGGGGTGGCCCCAGCTTTTGGGCTGCGCGTACCGGGGTTGCGTTCCCGCCCTTCGCGTGCACGCCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGt                                                                                                             <  6:368231/149‑1 (MQ=255)                                      ccAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCATACCCTTCAAACCCCCCCAGGGTTTGAACCCAACCCCaaaaaaataaatt                                                                                                            >  6:368349/1‑142 (MQ=255)                                              ggcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTa                                                                                                    <  5:413142/149‑1 (MQ=255)                                                     ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCATAGCGTGCAAGCCCGCCAAGGGTTTCAACCCCACCCCAAAAAAATAAGTGGAATTTTaaaaatt                                                                                              >  8:88515/1‑142 (MQ=255)                                                     ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                      >  1:348670/1‑106 (MQ=25)                                                     ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                      <  2:348670/108‑3 (MQ=25)                                                      cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCATAGCCGGGCAGCCCGCGAAGGGTTGGAACGCCACCCCAAAAAACTAAGTGCATTTTTACTATTTGc                                                                                            >  8:208911/1‑149 (MQ=255)                                                              tGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaac                                                                                                                                 <  1:339441/104‑1 (MQ=21)                                                              tGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaac                                                                                                                                 >  2:339441/1‑104 (MQ=21)                                                                                          ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTTCGTTTCAGCCCATAGCCTGGAAACCCCCGAAGGGTTGCAAAGCAACCCCAAAAAACAAAAGGCAATTTTTATATTTTATGATACGGCCCCCGCCACAAAAAACCGCAGGTGTc                                                         >  4:81942/1‑148 (MQ=255)                                                                                          ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCATAGCGTGCAAGCCCGCGCAGGGTTGGAAAGCAACCCCAAAAAAATCAATGCCTTTTTTATATTTTCTGATAGGACGCCCGCCAGATAAAGCGGCGGgttcga                                                        >  6:153291/1‑145 (MQ=255)                                                                                          ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                            <  5:153291/145‑1 (MQ=255)                                                                                               ccAGCTTTTGGGCTGAGCTTACCGGGGTTGCGTTCCAGCCCTTCGCGTCCAAGCCCGCGATGAGTTGCAACGCAACCCCAAAAAACt                                                                                                                 <  7:62198/87‑1 (MQ=255)                                                                                               ccAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                                   <  4:351168/149‑1 (MQ=255)                                                                                               ccAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                                   <  5:368349/149‑1 (MQ=255)                                                                                               ccAGCTTTTGGGCTGAGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                                   <  7:208911/149‑1 (MQ=255)                                                                                                 aGCTTTTGGGCTGAGCGTACCGGGGGTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCct                                                 <  7:88515/149‑1 (MQ=255)                                                                                                  gCTTTTGGGCTGAGCGTACCGGGGTTGCGTTCCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctg                                                <  3:339683/149‑1 (MQ=255)                                                                                                         ggCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                                         <  3:81942/149‑1 (MQ=255)                                                                                                              aGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                    <  1:394247/149‑1 (MQ=255)                                                                                                               gggTACCGGGGTTGCGTGGCAGCCCTTCGCGTGCAAGCCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                    <  4:210347/146‑1 (MQ=255)                                                                                                                ttttttggggTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCGAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                   >  4:262920/7‑148 (MQ=25)                                                                                                                 gTACCGGGGTTGCTTTTCAGCCCTTCGCGTGAAAGCCCGCTAAGGGTTGCAACGCAAACCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                 <  5:377564/149‑1 (MQ=255)                                                                                                                 gTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAACACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                 <  5:83079/149‑1 (MQ=255)                                                                                                                  tACCGGGGTTGCGTTGGAGCCCATAGCGTGCCCGCCCGCGCAGGGTTGGAACGCAACCCCCAAAAAAACAGTGCATTTTTACTATTTGCGGATAGGACCCCCCCCCCATCAAACCGCGGGTGGGAAAACCCAGCTGgcccgcccggcg                                 >  8:304135/1‑144 (MQ=255)                                                                                                                  tACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                    <  7:304135/145‑1 (MQ=255)                                                                                                                     cGGGGTTGCGTTGCAGCCCATAGCGTGCAAGCCCGCGCAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCAATTTTAATATTTGCTGATAGGACGCCCGCCAGATCAAGCGGCTGGTGGCATACCCCAGCTGCCGCGCCCAGGAcag                              >  2:85045/1‑148 (MQ=255)                                                                                                                      ggggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                            <  6:207412/149‑1 (MQ=255)                                                                                                                      ggggTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCGAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagc                             <  3:262920/148‑1 (MQ=37)                                                                                                                         gttgcgttgcAACCCTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                   >  5:339818/1‑107 (MQ=25)                                                                                                                         gttgcgttgcAACCCTTCGCGGGCTTCCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                   <  6:339818/107‑1 (MQ=25)                                                                                                                          ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt                        <  3:172929/149‑1 (MQ=255)                                                                                                                          ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGt                         >  4:287593/1‑148 (MQ=255)                                                                                                                              gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                              <  7:17947/123‑1 (MQ=255)                                                                                                                              gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                    >  1:400254/1‑149 (MQ=255)                                                                                                                              gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATATGCTGATCGGACGCCCGCCAGATCAAGCGCCTGGTGTCATAACCcagct                                              >  8:17947/1‑123 (MQ=255)                                                                                                                                tgcAACCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACg                  >  2:141901/1‑149 (MQ=255)                                                                                                                                  caacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt                >  7:351504/4‑149 (MQ=38)                                                                                                                                     cccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc              >  2:19463/1‑148 (MQ=39)                                                                                                                                       cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACg                  >  8:377823/1‑142 (MQ=255)                                                                                                                                       cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACg                  <  7:377823/142‑1 (MQ=255)                                                                                                                                       cTTCGCGGGCTTGCACGCGAAGCGCTGCAACGCAACACCAAAAAACTCAGTGCATTTTTACTTTTTGCTGATCTGACGCCCCCAAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCat           >  6:60565/1‑148 (MQ=255)                                                                                                                                        ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg          >  4:300465/1‑149 (MQ=39)                                                                                                                                        ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCAGTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCCCGCTCCCCgg          >  8:116933/1‑149 (MQ=38)                                                                                                                                           gcgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc        <  1:141901/148‑1 (MQ=255)                                                                                                                                            cgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg      <  1:85045/149‑1 (MQ=255)                                                                                                                                               gCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAg   >  2:234662/1‑149 (MQ=255)                                                                                                                                                cAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgt  >  4:197714/1‑149 (MQ=255)                                                                                                                                                 aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                 <  6:287437/117‑1 (MQ=255)                                                                                                                                                 aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                 >  5:287437/1‑117 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGT  >  NC_002947/5740326‑5740617

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