Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	346,119	T→A	5.9%	N576I (AAC→ATC)	PP_0287 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	346,119	0	T	A	5.9%	227.1 / 2.7	78	N576I (AAC→ATC)	PP_0287	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base A (3/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

GACGCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCGG  >  NC_002947/345971‑346267                                                                                                                                                     |                                                                                                                                                     gACGCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGt                                                                                                                                                      >  2:128597/1‑149 (MQ=255) gACGCCATTGCGTACCACCAGGCAGCCGCGCAGCTCCAGGAAGGCGGTGTCCATGCCCCGTGGTTCGCTGCTCAGCCTCTTGCGATTGAGCGTGGCGATCGCTTGGCACAACTGCTGTTCCAGGTTGGCAATGACCAGCACGCCATCGt                                                                                                                                                      >  5:148585/1‑149 (MQ=255)       aTTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATgg                                                                                                                                                >  3:246514/1‑149 (MQ=255)        ttGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGc                                                                                                                                               <  1:92881/149‑1 (MQ=255)         tGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCAACGATGGc                                                                                                                                               <  8:12617/148‑1 (MQ=255)           cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCAACGATGGCgtg                                                                                                                                            <  5:264837/149‑3 (MQ=255)            gTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAg                                                                                                                                           >  8:142971/1‑149 (MQ=255)            gTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTACTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAg                                                                                                                                           >  3:25704/1‑149 (MQ=255)               ccaccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGttt                                                                                                                                        <  1:343099/149‑1 (MQ=255)                 accaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgc                                                                                                                                      >  5:31702/1‑149 (MQ=255)                     ggCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCATCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTg                                                                                                                                  <  6:275898/149‑1 (MQ=255)                          ccGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGt                                                                                                                              >  8:92292/1‑148 (MQ=255)                              gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCGTCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                         <  7:335232/149‑1 (MQ=255)                              gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                         <  2:268975/149‑1 (MQ=255)                              gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                         <  4:246514/149‑1 (MQ=255)                                  cTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGgtgt                                                                                                                     <  7:215738/149‑1 (MQ=255)                                   tCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTc                                                                                                                    >  8:218967/1‑149 (MQ=255)                                     ctggttgggggTGTGCTTGGCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACATCTGCTGTTCCATGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTc                                                                                                                    <  6:148585/141‑1 (MQ=255)                                     cTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTc                                                                                                                    <  8:96257/147‑1 (MQ=255)                                      tGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                 >  3:101511/1‑149 (MQ=255)                                       ggAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                 >  2:35041/1‑148 (MQ=255)                                        gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgc                                                                                                                                      >  7:309968/1‑126 (MQ=255)                                        gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgc                                                                                                                                      <  8:309968/126‑1 (MQ=255)                                        gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                 >  5:289476/1‑147 (MQ=255)                                        gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                 >  4:362807/1‑147 (MQ=255)                                           gggggTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGAGTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTg                                                                                                            >  1:271287/1‑149 (MQ=255)                                             gggTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCt                                                                                                          <  7:92292/149‑1 (MQ=255)                                                tgtCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTc                                                                                                       <  2:315241/149‑1 (MQ=255)                                                  tCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTg                                                                                                     <  3:371875/149‑1 (MQ=255)                                                    cTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGAATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGc                                                                                                   <  6:273084/149‑1 (MQ=255)                                                       gCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                                  <  3:362807/147‑1 (MQ=255)                                                          gcgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTAc                                                                                             >  1:225363/1‑149 (MQ=255)                                                                ttCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTgg                                                                                       >  7:100006/1‑149 (MQ=255)                                                                  cGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGt                                                                                                             >  3:167489/1‑125 (MQ=255)                                                                  cGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGt                                                                                                             <  4:167489/125‑1 (MQ=255)                                                                    ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACATCTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGATGTTAGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                   >  2:178995/1‑149 (MQ=255)                                                                      gctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTg                                                                                                            <  7:143047/122‑1 (MQ=255)                                                                      gctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTg                                                                                                            >  8:143047/1‑122 (MQ=255)                                                                      gctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                   >  4:273060/1‑147 (MQ=255)                                                                       ctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACTATTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTc                                                                                  <  5:331649/147‑1 (MQ=255)                                                                         cAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGa                                                                               >  5:283237/1‑148 (MQ=255)                                                                          aGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGAt                                                                              <  7:132773/148‑1 (MQ=255)                                                                            cGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTTCTGTTCCAGGTTGGCATTGACGAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATc                                                                             <  4:25704/147‑1 (MQ=255)                                                                            cGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCAc                                                                           <  1:128597/149‑1 (MQ=255)                                                                                     aTTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcaggg                                                                  <  4:168868/149‑1 (MQ=255)                                                                                          gCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagggtcag                                                              >  7:301835/1‑148 (MQ=255)                                                                                          gCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagggtcag                                                              <  8:301835/148‑1 (MQ=255)                                                                                           cGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTcagca                                                            <  4:101511/149‑1 (MQ=255)                                                                                                      cTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGAcc                                                 >  5:340108/1‑149 (MQ=255)                                                                                                                ctgctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCt                                       <  6:283237/149‑1 (MQ=255)                                                                                                                  gctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTc                                                                                  >  6:108907/1‑104 (MQ=255)                                                                                                                  gctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTc                                                                                  <  5:108907/104‑1 (MQ=255)                                                                                                                      ttCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCt                                 >  5:112737/1‑149 (MQ=255)                                                                                                                        ccAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCtt                                <  2:225363/148‑1 (MQ=255)                                                                                                                         cAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCACTTGACCGGCGGCGTCTGTTCCTTGc                              >  7:124464/1‑149 (MQ=255)                                                                                                                            gTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCtt                           >  4:73986/1‑149 (MQ=255)                                                                                                                             ttGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTc                          <  6:289476/149‑1 (MQ=255)                                                                                                                              tGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCg                         >  3:160943/1‑149 (MQ=255)                                                                                                                               ggCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCg                         <  5:216400/148‑1 (MQ=255)                                                                                                                                gCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCCGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCgg                        <  3:53781/148‑1 (MQ=255)                                                                                                                                 cATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGGTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTc                      >  8:141866/1‑149 (MQ=255)                                                                                                                                  aTTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCt                     >  4:187052/1‑149 (MQ=255)                                                                                                                                  aTTGACCAGCACGCCAACGATGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTCTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGAGCGGCGGCGTCTGTTGCTTGCCTTCGGTCt                     >  8:180954/1‑149 (MQ=255)                                                                                                                                      aCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat                 >  8:336680/1‑149 (MQ=255)                                                                                                                                      aCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat                 >  3:332959/1‑149 (MQ=255)                                                                                                                                      aCCAGCACGCCATCGTTGAAGGCGAAGTTTGCGCAGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat                 >  6:74339/1‑149 (MQ=255)                                                                                                                                       ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgatg                <  8:100006/149‑1 (MQ=255)                                                                                                                                       ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgatg                <  1:35041/149‑1 (MQ=255)                                                                                                                                       ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgatg                <  7:235503/149‑1 (MQ=255)                                                                                                                                       ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat                 >  8:265789/1‑148 (MQ=255)                                                                                                                                        cAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCtt                                >  6:62737/1‑132 (MQ=255)                                                                                                                                        cAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCtt                                <  5:62737/132‑1 (MQ=255)                                                                                                                                         aGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGaa              <  7:265789/149‑1 (MQ=255)                                                                                                                                           ttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTGTGTTCCTTGCCTTCGGTCTTGATGAAgt            >  7:318066/3‑149 (MQ=255)                                                                                                                                           ttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgt            >  1:47482/3‑149 (MQ=255)                                                                                                                                              gCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGtt         <  6:31702/149‑1 (MQ=255)                                                                                                                                               ccATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTc        <  2:271287/149‑1 (MQ=255)                                                                                                                                                  tCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCAc     <  8:124464/149‑1 (MQ=255)                                                                                                                                                    gTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCg   <  2:17184/149‑1 (MQ=255)                                                                                                                                                    gTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCg   >  8:228188/1‑149 (MQ=255)                                                                                                                                                     ttGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCgg  >  8:232093/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GACGCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCGG  >  NC_002947/345971‑346267

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