Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,163,435	T→A	6.1%	I99N (ATT→AAT)	PP_5482 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,163,435	0	T	A	6.1%	196.2 / 2.5	65	I99N (ATT→AAT)	PP_5482	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base A (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

GGGGTAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAATTTTGG  >  NC_002947/2163287‑2163580                                                                                                                                                     |                                                                                                                                                   ggggTAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGcatcat                                                                                                                                                    <  1:379062/149‑1 (MQ=255)  gggTAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAAGGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCAtt                                                                                                                                                   >  1:285389/1‑149 (MQ=255)  gggTAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCAtt                                                                                                                                                   <  1:187732/149‑1 (MQ=255)    gTAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGa                                                                                                                                                 <  6:27017/149‑1 (MQ=255)     tAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCATCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGAt                                                                                                                                                <  6:5469/149‑1 (MQ=255)     tAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGAt                                                                                                                                                <  7:130353/149‑1 (MQ=255)           ggCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGGAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGCAACTCGAATTCGAGTATGAAGGCATCATTGATGCCttt                                                                                                                                          <  8:64982/149‑1 (MQ=255)            gCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTc                                                                                                                                         <  4:283718/149‑1 (MQ=255)             cAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCa                                                                                                                                        <  1:225497/149‑1 (MQ=255)             cAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCAATGATGCCttcat                                                                                                                                        <  4:153274/149‑4 (MQ=255)               aCTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCAATGATGCCttcatac                                                                                                                                      <  1:103186/149‑6 (MQ=255)                   gTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCtata                                                                                                                                  <  3:177245/149‑1 (MQ=255)                    tGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGTATCATTGATGCCTTTCAGCTatat                                                                                                                                 <  3:16860/149‑1 (MQ=255)                           aCTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGAAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGAc                                                                                                                           >  5:362646/1‑148 (MQ=255)                                 aTGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGc                                                                                                                     <  3:12541/148‑1 (MQ=255)                                      ttAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATgg                                                                                                                      >  6:148579/1‑142 (MQ=255)                                      ttAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATgg                                                                                                                      <  5:148579/142‑1 (MQ=255)                                       tAAATCGGTAGTCATTGAGCGTAATTACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTTGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTc                                                                                                                <  7:249899/147‑1 (MQ=255)                                          aaCGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCAtt                                                                                                           >  6:324941/1‑149 (MQ=255)                                                  tCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTg                                                                                                   <  7:37994/149‑1 (MQ=255)                                                    aTTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑a                                                                                                >  8:117306/1‑149 (MQ=255)                                                    aTTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑a                                                                                                >  5:383296/1‑149 (MQ=255)                                                    aTTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGAATAATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑a                                                                                                >  3:214287/1‑149 (MQ=255)                                                    aTTGAGCGAAATAACGATGAGCATGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑a                                                                                                >  3:155530/1‑149 (MQ=255)                                                       gAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATAc                                                                                             >  2:50283/1‑149 (MQ=255)                                                          cGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGCCCATGGCGACTCGCATTCGTTGCTGG‑ATACCAg                                                                                          >  8:118817/1‑149 (MQ=255)                                                              atatCGATGAGCTTTTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCg                                                                                      <  6:383296/149‑1 (MQ=255)                                                                     tGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGATTATGAAGGCATCATTGATGACTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTg                                                                                                   <  2:64245/130‑1 (MQ=255)                                                                     tGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTg                                                                                                   >  1:64245/1‑130 (MQ=255)                                                                          gtgtTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCg                                                                            >  4:352153/1‑147 (MQ=255)                                                                             ttATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCg                                                                                      >  6:193935/1‑134 (MQ=255)                                                                             ttATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCg                                                                                      <  5:193935/134‑1 (MQ=255)                                                                             ttATAGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGcc                                                                       <  8:325076/149‑1 (MQ=255)                                                                                  gCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCAc                                                                   <  8:215021/148‑1 (MQ=255)                                                                                      ccATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTaaa                                                              >  6:56531/1‑149 (MQ=255)                                                                                        aTCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAAtt                                                            <  1:50283/149‑1 (MQ=255)                                                                                          cAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCCCCTAAATTCg                                                          >  8:128493/1‑149 (MQ=255)                                                                                                tATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCgg                                                    >  4:159690/1‑149 (MQ=255)                                                                                                     ttGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCAttt                                               >  8:33163/1‑149 (MQ=255)                                                                                                       gACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTg                                              >  6:193870/1‑148 (MQ=255)                                                                                                        aCGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGa                                             <  4:214287/148‑1 (MQ=255)                                                                                                         cGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGaaa                                           <  8:335069/149‑1 (MQ=255)                                                                                                         agAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGaaa                                           <  6:362646/148‑1 (MQ=255)                                                                                                          gAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAAc                                          >  1:21992/1‑149 (MQ=255)                                                                                                          gAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTCTGCACAGCAGACGCACCTAAAATCTCACCGGAATTTTaaa                                           >  4:102918/1‑148 (MQ=255)                                                                                                            ggCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTa                                        >  6:142055/1‑149 (MQ=255)                                                                                                               tCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCATGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACgg                                     >  4:323484/1‑149 (MQ=255)                                                                                                               tCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACgg                                     >  5:75785/1‑149 (MQ=255)                                                                                                               tCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACgg                                     >  5:32115/1‑149 (MQ=255)                                                                                                                 gTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACgg                                     <  7:117306/147‑1 (MQ=255)                                                                                                                 gTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGtt                                   <  6:344419/149‑1 (MQ=255)                                                                                                                   aGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAg                                                                         >  2:317262/1‑109 (MQ=255)                                                                                                                   aGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAg                                                                         <  1:317262/109‑1 (MQ=255)                                                                                                                     gCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGaaa                               <  2:311093/149‑1 (MQ=255)                                                                                                                       aGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGt                             >  8:144552/1‑149 (MQ=255)                                                                                                                       aGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGt                             >  2:190824/1‑149 (MQ=255)                                                                                                                        ggAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGtt                            >  7:106768/1‑149 (MQ=255)                                                                                                                               gAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCAc                                                                              >  4:237987/1‑92 (MQ=255)                                                                                                                               gAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCAc                                                                              <  3:237987/92‑1 (MQ=255)                                                                                                                               gAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGc                     >  6:294503/1‑149 (MQ=255)                                                                                                                                  ttCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGC‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTa                  <  8:196397/149‑1 (MQ=255)                                                                                                                                     gAGTATGAAGGCATCATTGATGCCTTTGAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAACGTTTCAACGCTTacac               <  7:33163/149‑1 (MQ=255)                                                                                                                                     gAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTacac               <  5:294503/149‑1 (MQ=255)                                                                                                                                     gAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTacac               <  6:75785/149‑1 (MQ=255)                                                                                                                                       gTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCa             <  3:352153/149‑1 (MQ=255)                                                                                                                                        ctgaaaGGCATCAATGATGCCTCTCAGCTATATGTGGACCATGGCGCCTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAg            >  7:351491/5‑149 (MQ=255)                                                                                                                                         tgaaaGGCATCAATGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAg            >  6:172916/4‑148 (MQ=255)                                                                                                                                               gCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGAATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAAtt      >  1:103435/1‑149 (MQ=255)                                                                                                                                                 atcatTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAATTTTg   <  6:32115/149‑1 (MQ=255)                                                                                                                                                 atcatTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAATTTTg   >  6:181901/1‑149 (MQ=255)                                                                                                                                                 atcatTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAATTTTg   >  3:3762/1‑149 (MQ=255)                                                                                                                                                  tcatTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTGACACCAGTAAATTTTgg  >  6:258523/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGGGTAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGG‑ATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCAGCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAATTTTGG  >  NC_002947/2163287‑2163580

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