Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,081,094	A→T	6.5%	H250L (CAT→CTT)	PP_0938 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,081,094	0	A	T	6.5%	182.9 / 2.5	61	H250L (CAT→CTT)	PP_0938	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base T (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CCTGCCTGACGGCCAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCTGAATAATCTGGCTTT  >  NC_002947/1080948‑1081238                                                                                                                                                   |                                                                                                                                                 cctgcctgACGGCCAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCTGCCCAGGGTGTTGGGGCAGTGGCATg                                                                                                                                                >  7:262582/1‑149 (MQ=255) cctgcctgACGGCCAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATg                                                                                                                                                <  6:354208/149‑1 (MQ=255)         aCGGCCAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACAcc                                                                                                                                        >  5:399724/1‑149 (MQ=255)           ggCCAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCt                                                                                                                                      >  7:283598/1‑149 (MQ=255)            gCCAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATg                                                                                                                                                <  8:345557/138‑1 (MQ=255)            gCCAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATg                                                                                                                                                >  7:345557/1‑138 (MQ=255)             ccAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGGCTGCCGCAAAGTGACTGGCCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTgc                                                                                                                                    <  2:173872/149‑1 (MQ=255)              cAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTgcg                                                                                                                                   >  6:67884/1‑149 (MQ=255)              cAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTgcg                                                                                                                                   <  5:35919/149‑1 (MQ=255)               aGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGTCACCCTgcgg                                                                                                                                  <  3:360104/149‑1 (MQ=255)                   ccTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTgcgggcgg                                                                                                                              <  5:216336/149‑1 (MQ=255)                     tGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTgcgggcgggc                                                                                                                            >  1:286469/1‑149 (MQ=255)                        cGTTGAGCCTGACCAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGt                                                                                                                         <  4:75261/149‑1 (MQ=255)                             aGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCt                                                                                                                                              >  5:308977/1‑123 (MQ=255)                             aGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCt                                                                                                                                              <  6:308977/123‑1 (MQ=255)                              gCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTc                                                                                                                   <  5:67884/149‑1 (MQ=255)                               ccTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCt                                                                                                                  <  5:195257/149‑1 (MQ=255)                                  gAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGa                                                                                                               <  5:331482/149‑1 (MQ=255)                                    aCAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATc                                                                                                             <  3:347052/149‑1 (MQ=255)                                        ccGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCt                                                                                                         >  1:170314/1‑149 (MQ=255)                                        ccGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCt                                                                                                         >  4:272919/1‑149 (MQ=255)                                          gTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTgg                                                                                                       <  8:111427/149‑1 (MQ=255)                                               tCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCaa                                                                                                  <  5:161644/149‑1 (MQ=255)                                                  gCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAggg                                                                                               >  8:76576/1‑149 (MQ=255)                                                              cGTGATGGCGAGGTAGAATCCTATGTGTGCCTGCCGCAAAGTGACTGGGCAGATTGTTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGTTCGGCTGGGGCAAGGGGCAGTTGCAACg                                                                                   <  6:46070/149‑1 (MQ=255)                                                               gTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTgcgggc                                                                                                                                <  4:221605/103‑1 (MQ=255)                                                                  aTGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGcc                                                                               <  3:140094/149‑1 (MQ=255)                                                                     gCGTGGTAGAAGCCTATGTGAGCCTGGCGCAATGTGAGTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGcc                                                                               <  8:283598/146‑1 (MQ=255)                                                                      cGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCcgtcg                                                                           <  2:286469/149‑1 (MQ=255)                                                                         ggTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTaa                                                                        <  7:76576/149‑1 (MQ=255)                                                                                gCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCaa                                                                 >  1:53993/1‑149 (MQ=255)                                                                                gCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCaa                                                                 <  5:128598/149‑1 (MQ=255)                                                                                   tATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATg                                                               >  4:50158/1‑148 (MQ=255)                                                                                     tgtgAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCAc                                                            >  5:71875/1‑149 (MQ=255)                                                                                     tgtgAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCCGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCAc                                                            <  3:50158/149‑1 (MQ=255)                                                                                             tGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTg                                                    >  6:111024/1‑149 (MQ=255)                                                                                              gCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTg                                                    >  8:26224/1‑148 (MQ=255)                                                                                               ccGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTgg                                                   <  7:26224/148‑1 (MQ=255)                                                                                               ccGCAAAGGGCCTGGGCACATTGGTTGCCGCCCAGGGTGGTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAATGGGCAGTTGCAACTGGCCGTCGTAAGGCTCAATGCACCCAGGCTg                                                    <  6:399724/147‑1 (MQ=255)                                                                                                    aaGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAggggg                                             >  4:179831/1‑149 (MQ=255)                                                                                                         aCTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTAt                                        >  6:67275/1‑149 (MQ=255)                                                                                                           tGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGc                                      >  2:72526/1‑149 (MQ=255)                                                                                                             ggCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCTTGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCg                                    >  3:369402/1‑149 (MQ=255)                                                                                                              gCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCTTGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCgg                                   <  4:369402/149‑1 (MQ=255)                                                                                                                acaTTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCTTGCTGACACACTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGgc                                 >  6:250259/1‑149 (MQ=255)                                                                                                                   ttGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGgcgc                               >  7:93/1‑148 (MQ=255)                                                                                                                          ccgccCAGGGTGTTGGGGCAGTGGCTTGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGc                       <  5:250259/149‑1 (MQ=255)                                                                                                                             cccAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCaa                    <  6:83816/149‑1 (MQ=255)                                                                                                                              ccAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAAc                   <  2:53993/149‑1 (MQ=255)                                                                                                                               cAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGACCAGCACAAGGCTGACAAc                   >  6:269821/1‑148 (MQ=255)                                                                                                                                aGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCt                 <  5:66175/149‑1 (MQ=255)                                                                                                                                aGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCt                 <  2:170314/149‑1 (MQ=255)                                                                                                                                 gggTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTAGGGCAAGGGACAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCt                 >  2:260022/1‑148 (MQ=255)                                                                                                                                        ggggCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCTGAATAATc         <  8:262582/149‑1 (MQ=255)                                                                                                                                         gggCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCTGAATAATCt        >  1:354594/1‑149 (MQ=255)                                                                                                                                         gggCAGTGGCACGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCTGAATAATc         >  2:256941/1‑148 (MQ=255)                                                                                                                                         gggCAGTGGCACGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCTGAATAATc         <  1:256941/148‑1 (MQ=255)                                                                                                                                          ggCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCTGAATAATCTg       <  5:111024/149‑1 (MQ=255)                                                                                                                                             aGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCTGAATAATCTGGCt    >  7:162226/1‑149 (MQ=255)                                                                                                                                               tGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGAGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCTGAATAATCTGGCttt  >  1:87661/1‑149 (MQ=255)                                                                                                                                                 gCATGCTGACACCGTGCGGGCGGGGGGTCAGATCTCGATCCGCACGGCCAAGGGCCAGTTGCAACGGGCCGTCCTAAGCCTCAATGCACCCAGGCTGGAGGGCGGCTATGcccg                                   >  3:98880/1‑112 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CCTGCCTGACGGCCAGCCCCTGGCGTTGAGCCTGAACAGCCGTATGTCGGCCGAGGCCTGGCGTGATGGCGAGGTAGAAGCCTATGTGAGCCTGCCGCAAAGTGACTGGGCACATTGGTTGCCGCCCAGGGTGTTGGGGCAGTGGCATGCTGACACCCTGCGGGCGGGCGGTGAGTTCTGGATCGGCTGGGGCAAGGGGCAGTTGCAACGGGCCGTCGTAAGGCTCAATGCACCCAGGCTGGAGGGGGCCTATGCCGGGCGCAAGGCTGCCAACCTGAATAATCTGGCTTT  >  NC_002947/1080948‑1081238

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