Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,030,350	C→A	7.0%	intergenic (+35/+32)	PP_5729 → / ← rpmG	hypothetical protein/50S ribosomal protein L33

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,030,350	0	C	A	7.0%	211.4 / 2.8	69	intergenic (+35/+32)	PP_5729/rpmG	hypothetical protein/50S ribosomal protein L33
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base A (3/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCACCCGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCGGTGGTG  >  NC_002947/6030203‑6030487                                                                                                                                                      |                                                                                                                                          gCCGCACCCGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGCCCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCg                                                                                                                                          >  8:73842/1‑149 (MQ=255)        ccGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCgg                                                                                                                                   >  6:257724/1‑149 (MQ=255)        ccGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGCATGCgg                                                                                                                                   >  4:131541/1‑149 (MQ=255)          gggTGGTCGGCCACAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGt                                                                                                                                 >  4:186673/1‑149 (MQ=255)           ggtggtCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGtt                                                                                                                                >  4:243669/1‑149 (MQ=255)             tggtCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGtttt                                                                                                                              <  8:26386/149‑1 (MQ=255)               gtCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGttttt                                                                                                                             <  8:230216/148‑1 (MQ=255)                tCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCCAAAAAAAAA‑ACCCGCATCCGAAAGGAGGCGGGtttttt                                                                                                                            <  7:174994/149‑1 (MQ=38)                gcGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGttttttt                                                                                                                           <  4:79949/148‑1 (MQ=255)                  ggCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCATCCGAAAGGAGGCGGGttttttttt                                                                                                                         <  2:292964/149‑3 (MQ=38)                   gCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAg                                                                                                                        <  6:219139/149‑1 (MQ=255)                       cAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCATCCGAAAGGAGGCGGGtttttttttggctg                                                                                                                    >  7:29925/1‑142 (MQ=38)                        aaCGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGc                                                                                                                   <  6:164519/149‑1 (MQ=255)                           gAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGa                                                                                                                <  8:169640/149‑1 (MQ=255)                             aaTGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATc                                                                                                              <  3:285053/149‑1 (MQ=255)                                   cACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTAc                                                                                                        >  6:222067/1‑149 (MQ=255)                                      ccGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCgg                                                                                                                                   >  2:16974/1‑119 (MQ=255)                                      ccGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCgg                                                                                                                                   <  1:16974/119‑1 (MQ=255)                                      ccGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTa                                                                                                         <  2:174165/145‑1 (MQ=255)                                      ccGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTa                                                                                                         >  1:174165/1‑145 (MQ=255)                                         cgAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCA‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATc                                                                                                  >  1:281389/1‑149 (MQ=255)                                           aGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCtt                                                                                                <  7:164986/149‑1 (MQ=255)                                               cgacgTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCt                                                                                            <  1:162164/149‑1 (MQ=255)                                                   gTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCt                                                                                        <  3:243669/149‑1 (MQ=255)                                                    tCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCtt                                                                                       <  4:250232/149‑1 (MQ=255)                                                       aGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTa                                                                                    <  3:265083/149‑1 (MQ=255)                                                          gcATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCtt                                                                                       <  6:109397/143‑1 (MQ=255)                                                          gcATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCtt                                                                                       >  5:109397/1‑143 (MQ=255)                                                                    cAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTAAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTg                                                                       >  1:167934/1‑149 (MQ=255)                                                                     aaGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGc                                                                      <  7:152315/149‑1 (MQ=255)                                                                      aGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCg                                                                     <  3:186673/149‑1 (MQ=255)                                                                       ggCCGGTGCCGACCTGCTTATCCAGGGCTTCGACTTCAACTGAACCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATTCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGa                                                                    <  5:39078/149‑1 (MQ=255)                                                                       ggCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGa                                                                    <  5:222067/149‑1 (MQ=255)                                                                              gCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCg                                                             >  8:181481/1‑149 (MQ=255)                                                                               ccGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAAGACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCg                                                             >  1:84304/1‑149 (MQ=255)                                                                               ccGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑‑AAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAAcc                                                              >  3:221417/1‑146 (MQ=255)                                                                               ccGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑‑AAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAAcc                                                              <  4:221417/146‑1 (MQ=255)                                                                                cGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGa                                                           >  7:213000/1‑149 (MQ=255)                                                                                 gACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCAGCCTCCTTTAACATGCCGGTTTTTTTAATGTCGCCGATAAATTAAATAAAAATGGCTTCATTGTACACCACTTGATTGCGAACAACCAGa                                                           >  4:100968/1‑148 (MQ=255)                                                                                  aCCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGAt                                                          >  2:146715/1‑148 (MQ=255)                                                                                       ctgATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATAtt                                                    >  1:302228/1‑149 (MQ=255)                                                                                          aTCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCt                                                 <  8:213000/149‑1 (MQ=255)                                                                                           tCCAGTGCTTAGACTTCAACTGACCCTGCAACAGCA‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCtt                                                <  2:167934/149‑1 (MQ=255)                                                                                             cAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTg                                               <  3:131541/148‑1 (MQ=255)                                                                                             cAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTg                                               <  5:257724/148‑1 (MQ=255)                                                                                              aGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCCCCTCGTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTAATTGATCTTGGCTTCCTTGTAACCCACGTGCTTGCGAACAACAGGATCATATTTGTTGa                                              >  1:159670/1‑148 (MQ=255)                                                                                              aGGGCTTAGACTTAAACTGACCCTGAAACAGAC‑AAAAAAAA‑ACACGCCTCCATTCTGATGCCGGTATTTTTCAGGTCGCCGATCAATTACTTGATCTTTGCTTCCTAGTACACCCCGTGCTTCCTAAAAACCGGAACAAATATCTTGa                                              >  4:59717/1‑148 (MQ=255)                                                                                                   ttCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGa                                        >  2:297407/1‑149 (MQ=255)                                                                                                   ttCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCATTCTGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGa                                        <  2:159670/149‑1 (MQ=255)                                                                                                       aCTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGAtttt                                    <  3:154509/149‑1 (MQ=255)                                                                                                       aCTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGAtttt                                    >  4:6928/1‑149 (MQ=255)                                                                                                           cAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGt                                  <  2:281389/147‑1 (MQ=255)                                                                                                            aaCTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCg                               >  3:280289/1‑149 (MQ=255)                                                                                                                  cccTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGggtg                          <  2:302228/148‑1 (MQ=255)                                                                                                                    cTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGgtggtg                       >  6:255239/1‑149 (MQ=255)                                                                                                                       cAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTgc                      <  4:280289/147‑1 (MQ=255)                                                                                                                         aCAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTCGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCtt                  <  3:232895/149‑1 (MQ=255)                                                                                                                            gCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGtt               <  7:73842/149‑1 (MQ=255)                                                                                                                                aaaaaaaa‑aCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCggg                             <  7:2641/132‑1 (MQ=255)                                                                                                                                aaaaaaaa‑aCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCggg                             >  8:2641/1‑132 (MQ=255)                                                                                                                                aaaaaaaa‑aCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCtt            <  3:6928/149‑1 (MQ=255)                                                                                                                                 gaaaaaa‑aCCCGCATCCGAAAGGAGGCGGGATTTTTTCAGGTCGCCGATCAATTACTTGATCTTGCCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATGACGATTTTCTCCGGGGTGGTCCGCTTGTTCTTg           >  6:93248/2‑149 (MQ=255)                                                                                                                                   aaaaa‑aCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTc         >  7:283282/1‑149 (MQ=255)                                                                                                                                    aaaa‑aCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCg        >  5:173136/1‑149 (MQ=255)                                                                                                                                      aa‑aCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATc                                                         <  6:80689/98‑1 (MQ=255)                                                                                                                                      aa‑aCCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATc                                                         >  5:80689/1‑98 (MQ=255)                                                                                                                                          ccctgctcgtTTCGGTTGCGGGTTTTTTTCTGGTCGCCGATCAATTAGTTGATCTTGGCTTCCTTGTACACCACGTGCTTGGGAACAAGCGGTTCATATTTCTTGTTCTCGTTTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCggtggt   <  3:59717/144‑1 (MQ=255)                                                                                                                                          cccGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCggtggt   <  1:297407/149‑1 (MQ=255)                                                                                                                                              cctcctTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCGgtggtg  >  3:192424/1‑146 (MQ=255)                                                                                                                                                      |                                                                                                                                          GCCGCACCCGGGTGGTCGGCCCCAACGAAATGTGCACCCGCGAGTACGACGTCGAGCGCATCCAGATCAAGGCCGGTGCCGACCTGCTGATCCAGGGCTTCGACTTCAACTGACCCTGCAACAGCC‑AAAAAAAA‑ACCCGCCTCCTTTCGGATGCGGGTTTTTTTCAGGTCGCCGATCAATTACTTGATCTTGGCTTCCTTGTACACCACGTGCTTGCGAACAACCGGATCATATTTCTTGATCTCGATTTTGTCCGGGGTGGTGCGCTTGTTCTTGTCGGTGGTG  >  NC_002947/6030203‑6030487

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑