Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,554,883	C→T	7.1%	E5K (GAA→AAA) ‡	PP_4886 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,554,883	0	C	T	7.1%	249.5 / 6.9	84	E5K (GAA→AAA) ‡	PP_4886	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (44/34);  new base T (3/3);  total (47/37)
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.55e-01

GCACCGAACTTGTAGCCGATGATTTCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTC  >  NC_002947/5554738‑5555030                                                                                                                                                  |                                                                                                                                                    gCACCGAACTTGTAGCCGATGATTTCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGa                                                                                                                                                  <  1:201348/149‑1 (MQ=255)   aCCGAACTTGTAGCCGATGATTTCGATCCAGTTGGCGTACTGCTTATCCAGGGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTAGGCTGATAACGACGATGTTCAGa                                                                                                                                                  >  1:258405/1‑147 (MQ=255)        aCTTGTAGCCGATGATTTCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTCTGAACATCGTc                                                                                                                                           <  6:136441/149‑1 (MQ=255)        aCTTGTAGCCGATGATTTCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTCTGAACATCGTc                                                                                                                                           >  5:136441/1‑149 (MQ=255)          ttGTAGCCGATGATTTCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTc                                                                                                                                           <  4:285376/147‑1 (MQ=255)          ttGTAGCCGATGATTTCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTc                                                                                                                                           >  3:285376/1‑147 (MQ=255)             tAGCCGATGATTTCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATcc                                                                                                                                       <  4:252179/148‑1 (MQ=255)                    tgatTTCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGtt                                                                                                                                >  7:303505/1‑148 (MQ=255)                       tttCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTc                                                                                                                            <  2:20285/149‑1 (MQ=255)                        ttCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTcc                                                                                                                           >  3:92204/1‑149 (MQ=255)                             tCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAg                                                                                                                      >  6:111016/1‑149 (MQ=255)                                aGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTTCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATc                                                                                                                                              >  8:110006/1‑122 (MQ=255)                                aGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTTCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATc                                                                                                                                              <  7:110006/122‑1 (MQ=255)                                aGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAg                                                                                                                      >  7:119858/1‑146 (MQ=255)                                aGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAg                                                                                                                      <  8:119858/146‑1 (MQ=255)                                   tGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAAc                                                                                                                >  6:246727/1‑149 (MQ=255)                                           tGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCg                                                                                                        >  6:293790/1‑149 (MQ=255)                                           tGCTTATACAGTGCTTCGCCAGTGATTTCGGCGATATGATCATAAGCATCGAAGGCCATGCTTACCGCGCCCTGGGAAAGGTTCGGCAGATACAGACGATGTTCAGACATCGTCATCCTGAGGTTGATCCCGAAGTAAACAgccggccg                                                                                                        >  8:298714/1‑149 (MQ=255)                                             cTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCgg                                                                                                       >  2:88714/1‑148 (MQ=255)                                              ttATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTa                                                                                                     >  3:154530/1‑149 (MQ=255)                                              ttATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTa                                                                                                     >  4:3815/1‑149 (MQ=255)                                               tATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTCTGAACATCGTcg                                                                                                                                          <  6:96090/111‑2 (MQ=255)                                               tATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTCTGAACATCGTcg                                                                                                                                          >  5:96090/1‑110 (MQ=255)                                                aTCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTa                                                                                                     <  4:92204/147‑1 (MQ=255)                                                 tCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTacc                                                                                                   >  6:123989/1‑148 (MQ=255)                                                  ccAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTaccga                                                                                                 >  7:196108/1‑149 (MQ=255)                                                       gCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTCTGAACATCGTc                                                                                                                                           <  7:161146/102‑1 (MQ=255)                                                       gCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTCTGAACATCGTc                                                                                                                                           >  8:161146/1‑102 (MQ=255)                                                           cGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTaccga                                                                                                 <  6:119031/140‑1 (MQ=255)                                                           cGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTaccga                                                                                                 >  5:119031/1‑140 (MQ=255)                                                           cGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCtt                                                                                          >  2:170886/1‑147 (MQ=255)                                                            gCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAg                                                                                       <  8:303505/149‑1 (MQ=255)                                                             ccAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAgg                                                                                      >  1:146790/1‑149 (MQ=255)                                                                    tttCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATa                                                                               >  5:255577/1‑149 (MQ=255)                                                                        cgcgATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATacaa                                                                            >  6:206520/1‑148 (MQ=255)                                                                         gcgATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgttt                                                                                        >  5:229048/1‑131 (MQ=255)                                                                         gcgATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgttt                                                                                        <  6:229048/135‑5 (MQ=255)                                                                          cgATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAacac                                                                         <  5:206520/149‑1 (MQ=255)                                                                             tatGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTaa                                                                      <  6:255577/149‑1 (MQ=255)                                                                               tGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAAcc                                                                    <  8:196108/149‑1 (MQ=255)                                                                                 aaCATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgtt                                                                                         <  4:229185/126‑4 (MQ=255)                                                                                 aaCATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgtt                                                                                         >  3:229185/1‑123 (MQ=255)                                                                                  aCATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCt                                                                 <  5:246727/149‑1 (MQ=255)                                                                                    ttAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTAt                                                               <  3:269315/148‑1 (MQ=255)                                                                                      aaGCATCGAAAGCCATTCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATTTTCAGACATCGTAATCCTGAGGTTGCTACCGAAGTAAAAAGGAGGCCGGTACCGACCGCCTTCAGGCTCAATACAACAATAAACCCTATga                                                             >  4:269315/1‑148 (MQ=255)                                                                                      aaGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGCCATCGTCATCCAGAGGTTGCTCCCGAAGTAAAGCGGCGGCCGGACCCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATg                                                              >  4:278926/1‑148 (MQ=255)                                                                                      aaGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCa                                                                                        >  6:40752/1‑122 (MQ=255)                                                                                      aaGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCa                                                                                        <  5:40752/122‑1 (MQ=255)                                                                                      aaGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATg                                                              >  4:195296/1‑148 (MQ=255)                                                                                      aaGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGt                                                             >  3:7797/1‑149 (MQ=255)                                                                                        gCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGt                                                             >  5:67228/1‑147 (MQ=255)                                                                                          aTCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGc                                                         <  5:293790/149‑1 (MQ=255)                                                                                           tCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGc                                                         <  5:185349/148‑1 (MQ=255)                                                                                           tCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGc                                                         >  6:185349/1‑148 (MQ=255)                                                                                               aaGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTg                                                    <  1:88714/149‑1 (MQ=255)                                                                                                  ccATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCaa                                                 <  3:3815/149‑1 (MQ=255)                                                                                                   cATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAAc                                                <  3:195296/149‑1 (MQ=255)                                                                                                    aTGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAAGCTCAATACAACACAAAACCATATGTCAGACCTTGAAAc                                                >  2:335236/1‑148 (MQ=255)                                                                                                          aCCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGc                                                                                     >  2:117304/1‑105 (MQ=255)                                                                                                          aCCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGc                                                                                     <  1:117304/105‑1 (MQ=255)                                                                                                          aCCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACa                                               <  7:262652/143‑1 (MQ=255)                                                                                                          aCCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACa                                               >  8:262652/1‑143 (MQ=255)                                                                                                            ctctCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGc                                       >  2:129003/1‑149 (MQ=255)                                                                                                                cccTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTa                                   >  1:24387/1‑149 (MQ=255)                                                                                                                 ccTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTa                                   >  5:12997/1‑148 (MQ=255)                                                                                                                 ccTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTAt                                  <  2:146790/149‑1 (MQ=255)                                                                                                                     gTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAg                              <  6:67228/149‑1 (MQ=255)                                                                                                                               ggCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTgg                    <  8:262868/149‑1 (MQ=255)                                                                                                                                 cTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTggg                   >  5:277292/1‑148 (MQ=255)                                                                                                                                  tGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAg                              <  2:114630/136‑1 (MQ=255)                                                                                                                                  tGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGa                  >  3:190532/1‑148 (MQ=255)                                                                                                                                  tGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATATCAGCCCTTGCAACAGCAAACGCCCTATCAAg                              >  1:114630/1‑136 (MQ=255)                                                                                                                                   gATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGaaa                >  4:123181/1‑149 (MQ=255)                                                                                                                                       acgacgATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGc                                             >  4:277624/1‑116 (MQ=255)                                                                                                                                       acgacgATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGc                                             <  3:277624/116‑1 (MQ=255)                                                                                                                                        cgacgaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAATGAATATTTGGGAAATCTcc           <  4:7797/149‑1 (MQ=255)                                                                                                                                         gacgaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTcc           >  5:172802/1‑148 (MQ=255)                                                                                                                                          acgaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTa         <  2:24387/149‑1 (MQ=255)                                                                                                                                          acgaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTa         <  6:12997/149‑1 (MQ=255)                                                                                                                                           cgaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTAc        <  1:170886/149‑1 (MQ=255)                                                                                                                                            gaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTa         >  7:289011/1‑147 (MQ=255)                                                                                                                                            gaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACa       <  5:123989/149‑1 (MQ=255)                                                                                                                                                 tCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGTCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTc  >  5:155464/1‑149 (MQ=255)                                                                                                                                                 tCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTc  >  3:99022/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GCACCGAACTTGTAGCCGATGATTTCGATCCAGTTGGCGTACTGCTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTC  >  NC_002947/5554738‑5555030

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