Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,149,440	C→A	8.0%	intergenic (‑426/+22)	PP_mr51 ← / ← PP_4535	C4 ncRNA/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,149,440	0	C	A	8.0%	276.2 / 11.1	87	intergenic (‑426/+22)	PP_mr51/PP_4535	C4 ncRNA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (43/37);  new base A (5/2);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

TATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGTTTC  >  NC_002947/5149292‑5149587                                                                                                                                                      |                                                                                                                                                      tATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATc                                                                                                                                                       <  5:401/149‑1 (MQ=255)       cgcgTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATccc‑cccc                                                                                                                                                <  6:157740/149‑1 (MQ=255)           tCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGt                                                                                                                                            <  3:266820/149‑1 (MQ=255)             ttGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑g                                                                                                                                         >  2:268916/1‑149 (MQ=255)               gAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGGGGGGATAGC‑CTCCTGGtt                                                                                                                                           <  8:68232/146‑1 (MQ=255)               gAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑g                                                                                                                                         <  4:233340/147‑1 (MQ=255)                aGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGa                                                                                                                                       <  8:249741/148‑1 (MQ=255)                 gTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAAc                                                                                                                                     >  5:26229/1‑149 (MQ=255)                 gTCGCTGTGCATTCTGAATGACCACGAGCAATCACACGCATAATACATAATGGGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGGCCAAGCCAATCCTCCTCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAAc                                                                                                                                     <  3:44270/149‑1 (MQ=255)                  tCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACA‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACt                                                                                                                                    >  1:169945/1‑149 (MQ=255)                    gCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTc                                                                                                                                  <  6:94097/149‑1 (MQ=255)                    gCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTc                                                                                                                                  <  7:46208/149‑1 (MQ=255)                          cATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTa                                                                                                                             >  5:137540/1‑148 (MQ=255)                            ttCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTc                                                                                                                          <  5:99064/149‑1 (MQ=255)                                 aaTGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGGGGGGATAGC‑CTCCTGGtt                                                                                                                                           <  6:100111/128‑1 (MQ=38)                                 aaTGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGGGGGGATAGC‑CTCCTGGtt                                                                                                                                           >  5:100111/1‑128 (MQ=38)                                  aTGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAAt                                                                                                                    <  6:31517/149‑1 (MQ=255)                                   tGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATc                                                                                                                   <  6:340776/149‑1 (MQ=255)                                    gACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCa                                                                                                                  >  3:32729/1‑149 (MQ=255)                                      ccACGAGCAATCAGCCGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑CGAACTTCAGTTACTCGGCAATCAGa                                                                                                                >  4:127709/1‑149 (MQ=255)                                      ccACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGa                                                                                                                <  2:139581/149‑1 (MQ=255)                                       cACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGAc                                                                                                               >  5:147351/1‑149 (MQ=255)                                                tCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGACTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGAc                                                                                                               >  1:197964/1‑140 (MQ=255)                                                tCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGACTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGAc                                                                                                               <  2:197964/140‑1 (MQ=255)                                                               aTAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACt                                                                                                                                    >  3:256240/1‑104 (MQ=255)                                                               aTAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACt                                                                                                                                    <  4:256240/104‑1 (MQ=255)                                                                 aaTGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGACCTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGAc                                                                                      >  7:189113/1‑148 (MQ=255)                                                                 aaTGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACAGCGCGTTCTCCACCGCGAGCGGAc                                                                                      >  8:104934/1‑148 (MQ=255)                                                                   tgtgTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCTGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGaa                                                                                   <  6:137540/149‑1 (MQ=255)                                                                     tgtCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAg                                                                                  <  4:32729/148‑1 (MQ=255)                                                                       tCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTcc                                                                               >  3:318385/1‑149 (MQ=255)                                                                       tCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTcc                                                                               <  3:127709/149‑1 (MQ=255)                                                                        ccACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCt                                                                              >  5:63847/1‑149 (MQ=255)                                                                         cACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTc                                                                             >  3:236411/1‑149 (MQ=255)                                                                              cgcccgccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGa                                                                        >  4:142877/1‑149 (MQ=255)                                                                                 ccgccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCgg                                                                                                                        >  3:90220/1‑98 (MQ=255)                                                                                 ccgccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCgg                                                                                                                        <  4:90220/98‑1 (MQ=255)                                                                                 ccgccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTc                                                                             <  3:85651/141‑1 (MQ=255)                                                                                 ccgccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTc                                                                             >  4:85651/1‑141 (MQ=255)                                                                                 ccgccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACa                                                                      >  7:95243/1‑148 (MQ=255)                                                                                 ccgccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACAt                                                                     <  2:150896/149‑1 (MQ=255)                                                                                  cgccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATa                                                                    <  6:147351/149‑1 (MQ=255)                                                                                    ccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGt                                                                  >  1:117444/1‑149 (MQ=255)                                                                                         tCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGcc                                                              <  3:60979/148‑1 (MQ=255)                                                                                           ttGGCTTCAACTTCAGTGTCCAAGCCAATCCTCTTCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTg                                                           >  3:35128/1‑149 (MQ=255)                                                                                            tGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTg                                                           >  5:112076/1‑148 (MQ=255)                                                                                             ggCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCt                                                            <  6:132763/146‑1 (MQ=255)                                                                                                ttCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCGAGGAGCCTAACCC‑CCCCGGGTT‑GGAACTTCAGTTAAACGGAAATCAAACCCCGCGTTGTCCACCGAGCGCGGACGAAGACGACCAAGACATAGTTGCCATGCAAc                                                       >  1:240810/1‑149 (MQ=255)                                                                                                ttCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACg                                                      >  5:99313/1‑149 (MQ=255)                                                                                                   aaCTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCCGCCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACgg                                                     >  1:263467/1‑148 (MQ=255)                                                                                                   aaCTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACgg                                                     >  3:15989/1‑147 (MQ=255)                                                                                                   aaCTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGc                                                    <  6:26229/148‑1 (MQ=255)                                                                                                     cTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAgg                                                 <  2:117444/149‑1 (MQ=255)                                                                                                         aGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCAt                                             <  4:318385/149‑1 (MQ=255)                                                                                                            gtCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATcc                                           <  4:236411/148‑1 (MQ=255)                                                                                                               cAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGt                                       >  6:157425/1‑149 (MQ=255)                                                                                                               cAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGt                                       >  1:409/1‑149 (MQ=255)                                                                                                               cAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGt                                       >  8:260944/1‑149 (MQ=255)                                                                                                                aaGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGt                                       >  4:41470/1‑148 (MQ=255)                                                                                                                    cAATCCTCTTCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCg                                  <  4:35128/149‑1 (MQ=255)                                                                                                                         cTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGa                              <  3:41470/148‑1 (MQ=255)                                                                                                                         cTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGa                              <  2:263467/148‑1 (MQ=255)                                                                                                                         cTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGAc                             <  6:112076/149‑1 (MQ=255)                                                                                                                           cATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACaa                           >  5:321736/1‑149 (MQ=255)                                                                                                                              ccAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGa                        >  6:260047/1‑149 (MQ=255)                                                                                                                                aaGTTCCAACT‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGccc                                                             >  6:105138/1‑110 (MQ=255)                                                                                                                                aaGTTCCAACT‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGccc                                                             <  5:105138/110‑1 (MQ=255)                                                                                                                                aaGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTTGGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGAt                       >  1:52008/1‑149 (MQ=255)                                                                                                                                aaGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGAtt                      >  8:174930/1‑149 (MQ=255)                                                                                                                                aaGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGAtt                      <  5:260047/149‑1 (MQ=255)                                                                                                                                aaGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGAtt                      <  2:52008/149‑1 (MQ=255)                                                                                                                                  gTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTg                    <  7:260944/149‑1 (MQ=255)                                                                                                                                   ttCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGc                   >  3:241063/1‑149 (MQ=255)                                                                                                                                   ttCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGc                   >  2:225537/1‑149 (MQ=255)                                                                                                                                     ccAACC‑AGGGGGGGATAGC‑CTCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCgg                 >  5:290200/1‑149 (MQ=39)                                                                                                                                     ccAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTACTCGAAGACATAGTTGCCCTGAAACGGCAGGCCATCCGCGTGATCGCACACAATGATTTGCgg                 >  6:122606/1‑149 (MQ=255)                                                                                                                                      cAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGc                >  1:65168/1‑149 (MQ=255)                                                                                                                                      cAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGc                <  4:15989/149‑1 (MQ=255)                                                                                                                                      cAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGc                >  6:222024/1‑149 (MQ=255)                                                                                                                                      cAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGAGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGc                >  6:35614/1‑149 (MQ=255)                                                                                                                                       aaCC‑AGGGGGGGATAGC‑CTCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGCCATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGACCGCAGACAATGAGTTGCGGcg               >  8:266581/1‑148 (MQ=255)                                                                                                                                       aaCC‑AGGGGGGGATAGC‑CTCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt               >  3:98788/1‑149 (MQ=255)                                                                                                                                       aaCC‑AGGGGGGGATAGC‑CTCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt               >  8:178619/1‑149 (MQ=255)                                                                                                                                        aCC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCtt              >  7:286894/1‑149 (MQ=255)                                                                                                                                        aCC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCtt              <  1:268916/149‑1 (MQ=255)                                                                                                                                               aggCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGcc        <  4:241063/149‑1 (MQ=255)                                                                                                                                                 gCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCgg      >  3:217042/1‑149 (MQ=255)                                                                                                                                                    aTCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGttt   >  3:35071/1‑149 (MQ=255)                                                                                                                                                     tCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGTTTc  <  2:169945/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACC‑AGGAGGCTATCCC‑CCCCTGGTT‑GGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGTTTC  >  NC_002947/5149292‑5149587

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