Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,622,740	G→A	6.2%	P206L (CCG→CTG)	PP_1422 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,622,740	0	G	A	6.2%	190.1 / 3.3	64	P206L (CCG→CTG)	PP_1422	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base A (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGC  >  NC_002947/1622595‑1622888                                                                                                                                                   |                                                                                                                                                     gCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCgg                                                                                                                                                     >  8:126913/1‑147 (MQ=255)  ccAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAg                                                                                                                                                  >  4:42318/1‑149 (MQ=255)  ccAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAg                                                                                                                                                  <  2:140419/149‑1 (MQ=255)   cAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAgg                                                                                                                                                 <  2:294002/149‑1 (MQ=255)   cAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAgg                                                                                                                                                 >  1:51695/1‑149 (MQ=255)   cAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAgg                                                                                                                                                 >  6:208106/1‑149 (MQ=255)    aaaCGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGt                                                                                                                                                >  4:245840/1‑149 (MQ=255)         gCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCATGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGAc                                                                                                                                            >  5:225092/1‑148 (MQ=255)             ccccAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACt                                                                                                                                       >  6:282153/1‑149 (MQ=255)                        ggTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAAc                                                                                                                            <  3:245840/149‑1 (MQ=255)                            cGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTccc                                                                                                                        <  3:357965/149‑1 (MQ=255)                               cGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGAc                                                                                                                                            <  7:159506/126‑1 (MQ=255)                               cGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGAc                                                                                                                                            >  8:159506/1‑126 (MQ=255)                                ggCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGc                                                                                                                    >  1:356659/1‑149 (MQ=255)                                     cgcagcCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCATGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGa                                                                                                               <  6:225092/149‑1 (MQ=255)                                     cgcagcCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCtt                                                                                                                 >  3:8265/1‑147 (MQ=255)                                     cgcagcCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGa                                                                                                               <  5:208106/149‑1 (MQ=255)                                      gcagcCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAg                                                                                                                      >  2:127995/1‑141 (MQ=255)                                      gcagcCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAg                                                                                                              <  4:8265/149‑1 (MQ=255)                                            aaTGGCCGGTGACGGCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCtt                                                                                                        <  3:40514/149‑1 (MQ=255)                                            aaTGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCtt                                                                                                        <  1:332314/149‑1 (MQ=255)                                                gCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTcc                                                                                                    <  5:63428/149‑1 (MQ=255)                                                         gCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCAcc                                                                                           <  5:282153/149‑1 (MQ=255)                                                               gcTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTGTGCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTg                                                                                      >  1:101593/1‑149 (MQ=255)                                                                ctgactgaTGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGa                                                                                    <  2:356659/149‑1 (MQ=255)                                                                     tgaTGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGc                                                                               >  4:165416/1‑149 (MQ=255)                                                                     tgaTGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGc                                                                               >  8:240589/1‑149 (MQ=255)                                                                     tgaTGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGc                                                                               >  2:369952/1‑149 (MQ=255)                                                                      gaTGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGc                                                                               <  1:289115/148‑1 (MQ=255)                                                                          cAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGa                                                                           <  7:126913/148‑1 (MQ=255)                                                                          cAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGa                                                                           <  6:87079/148‑1 (MQ=255)                                                                                  gATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAAc                                                                                                                            >  1:182160/1‑91 (MQ=255)                                                                                  gATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAAc                                                                                                                            <  2:182160/91‑1 (MQ=255)                                                                                        cACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTTCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAg                                                              <  5:97979/147‑1 (MQ=255)                                                                                        cACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAgg                                                             >  6:266614/1‑148 (MQ=255)                                                                                        cACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGc                                                            <  1:18909/149‑1 (MQ=255)                                                                                                 gatgtAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCATGCCTTGAGGATCTTGTCCCACATCAACGTCTGGAGTTGCTCGATAACATAATGCAGGAAAAACggct                                                   >  5:301636/1‑147 (MQ=255)                                                                                                 gatgtAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCAGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGt                                                   >  4:287616/1‑149 (MQ=255)                                                                                                 gatgtAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCAGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGt                                                   >  4:282576/1‑149 (MQ=255)                                                                                                  atgtAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTc                                                  <  1:369952/149‑1 (MQ=255)                                                                                                     tAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGa                                               <  4:105487/149‑1 (MQ=255)                                                                                                      aaTCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGa                                               >  4:337593/1‑148 (MQ=255)                                                                                                       aTCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGac                                              >  7:215986/1‑148 (MQ=255)                                                                                                        tCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGaca                                             >  5:147730/1‑148 (MQ=255)                                                                                                            tCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGt                                        >  2:269811/1‑149 (MQ=255)                                                                                                              ttGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTc                                      >  2:334332/1‑149 (MQ=255)                                                                                                                gTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGacaaca                                          <  2:51695/143‑1 (MQ=255)                                                                                                                   aGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTtg                                  <  3:42318/148‑1 (MQ=255)                                                                                                                          gCGGTCCAGGTCATGAACCTGGCTCAGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCtt                           <  3:282576/149‑1 (MQ=255)                                                                                                                          gCGGTCCAGGTCATGAACCTGGCTCAGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCtt                           <  3:287616/149‑1 (MQ=255)                                                                                                                            ggTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCg                         <  8:215986/149‑1 (MQ=255)                                                                                                                             gTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCgg                        >  6:97035/1‑149 (MQ=255)                                                                                                                                 aGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATa                    >  4:78312/1‑149 (MQ=255)                                                                                                                                    tCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGAt                 <  3:165416/149‑1 (MQ=255)                                                                                                                                    tCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGAt                 <  3:337593/149‑1 (MQ=255)                                                                                                                                      aTGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTTGAGTTGCTCGATCAGATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATg                <  5:266614/148‑1 (MQ=255)                                                                                                                                      aTGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGAt                 >  3:302574/1‑147 (MQ=255)                                                                                                                                      aTGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGa               <  3:78312/149‑1 (MQ=255)                                                                                                                                      aTGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGa               >  7:96850/1‑149 (MQ=255)                                                                                                                                       tGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAg              >  6:334880/1‑149 (MQ=255)                                                                                                                                        gAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGAt                                                                          >  5:253701/1‑88 (MQ=255)                                                                                                                                        gAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGAt                                                                          <  6:253701/88‑1 (MQ=255)                                                                                                                                                 tCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCTCCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGAt                 <  3:309802/136‑1 (MQ=255)                                                                                                                                                 tCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTATTCGGCATAGAt                 >  4:309802/1‑136 (MQ=255)                                                                                                                                                  cGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATACATGAGTACCGGCAgcg   >  1:155445/1‑149 (MQ=255)                                                                                                                                                   ggCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGCGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAgcgc  >  3:222479/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGATGCAGCGGGCGATCTGCACGATGTCGATGTAATCGATCTTGTCAGTCTG‑GCGGTCCAGGTCATGAACCTGGCTCGGCAGGTGGACCAACTGCTCCGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGC  >  NC_002947/1622595‑1622888

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