Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,421,568	G→A	8.4%	G127R (GGG→AGG)	moaB‑I →	molybdenum cofactor biosynthesis protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,421,568	0	G	A	8.4%	162.1 / 3.2	59	G127R (GGG→AGG)	moaB‑I	molybdenum cofactor biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base A (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCTGTGA  >  NC_002947/2421420‑2421716                                                                                                                                                      |                                                                                                                                                      cggcggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCg                                                                                                                                                       <  8:189661/149‑1 (MQ=255) cggcggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCg                                                                                                                                                       <  6:14606/149‑1 (MQ=255)   gcggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCggg                                                                                                                                                     >  6:85495/1‑149 (MQ=255)    cggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGc                                                                                                                                                    <  6:179096/149‑1 (MQ=255)      gTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTg                                                                                                                                                  <  4:199906/149‑1 (MQ=255)       tACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTgg                                                                                                                                                 <  6:21826/149‑1 (MQ=255)                tACTGGTCGCGACAGCACCCCCGCGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaaggc                                                                                                                                          <  5:34388/148‑4 (MQ=255)                   tGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGc                                                                                                                                     >  3:114647/1‑149 (MQ=255)                     gTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCAA‑GGCACGcga                                                                                                                                   <  1:302963/149‑3 (MQ=255)                       cgcgACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGgt                                                                                                                                 <  1:344991/149‑1 (MQ=255)                        gcgACAGCACCCCCGAGGCGGTGGACTGCCTGCTGGACAATCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGAACggttgtt                                                                                                                                >  7:57737/1‑143 (MQ=255)                          gACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCTAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGtgtg                                                                                                                              >  5:374312/1‑149 (MQ=255)                                aCCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGctgcct                                                                                                                        >  7:5216/1‑149 (MQ=255)                                  ccccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGctgcctgc                                                                                                                      >  1:326330/1‑149 (MQ=255)                                   cccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaag                                                                                                                                            >  4:158030/1‑126 (MQ=255)                                    ccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCa                                                                                                                    <  8:55818/149‑1 (MQ=255)                                    ccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAAAGGGCCGCTGGTGTGctgcctgc                                                                                                                      >  3:195895/1‑147 (MQ=255)                                           ggTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTcca                                                                                                             >  5:213650/1‑149 (MQ=255)                                              ggCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGGCCGCTGGTGTGGTGCCTGCCCGGCTccacc                                                                                                           >  1:383089/1‑148 (MQ=255)                                                 ctgcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaaggca                                                                                                                                         <  1:209379/116‑5 (MQ=38)                                                 ctgcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaaggca                                                                                                                                         >  2:209379/1‑112 (MQ=38)                                                   gcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCt                                                                                                     >  4:231679/1‑149 (MQ=255)                                                       gctgGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGc                                                                                                 >  3:108052/1‑149 (MQ=255)                                                       gctgGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGc                                                                                                 >  4:98209/1‑149 (MQ=255)                                                           gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTg                                                                                                                                                  <  2:27396/96‑1 (MQ=255)                                                           gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTg                                                                                                                                                  >  1:27396/1‑96 (MQ=255)                                                           gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCATGCTCCACCAACGCTGTGCGTAc                                                                                             >  2:239925/1‑149 (MQ=255)                                                           gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTAc                                                                                             >  4:81055/1‑149 (MQ=255)                                                           gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTAc                                                                                             >  1:231193/1‑149 (MQ=255)                                                                gcaagTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGtt                                                                                        >  3:311270/1‑149 (MQ=255)                                                                    gTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGa                                                                                    >  8:222836/1‑149 (MQ=255)                                                                       gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGGC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                                  <  5:85495/148‑1 (MQ=255)                                                                       gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                                  <  7:222836/148‑1 (MQ=255)                                                                           gCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATc                                                                             <  2:326330/149‑1 (MQ=255)                                                                             tttGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCt                                                                           >  5:93215/1‑149 (MQ=255)                                                                                 gTGAGCTGTCCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGc                                                                        <  6:213650/148‑1 (MQ=255)                                                                                  tGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCg                                                                      <  4:108052/149‑1 (MQ=255)                                                                                    aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGag                                                                    <  5:210384/149‑1 (MQ=255)                                                                                    aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGa                                                                     <  4:124599/148‑1 (MQ=255)                                                                                         ttCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGagcagc                                                                <  2:385297/148‑1 (MQ=255)                                                                                         ttCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TGTGCCGGGCGGGCCAAACGCACCCGGGTGTGCTTCCAGCCCCGCTCCACCAACGCTGTGGGTTCCGGGTGGGGTTGCCTGCTGGCCGCGCCGCt                                                               >  1:300476/1‑149 (MQ=255)                                                                                           ccGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTc                                                              <  3:231679/148‑1 (MQ=255)                                                                                            cGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGa                                                            >  3:67840/1‑149 (MQ=255)                                                                                              gCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑g                                                         <  6:374312/149‑1 (MQ=255)                                                                                                 gATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGTCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑Gccc                                                      <  4:311270/149‑1 (MQ=255)                                                                                                   tttCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGt                                                    <  2:130834/149‑1 (MQ=255)                                                                                                     tCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTc                                                   <  4:114647/148‑1 (MQ=255)                                                                                                      ccGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCAt                                                 >  3:315297/1‑149 (MQ=255)                                                                                                            ccGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCAGCCATGCTCCACCAACCCTTTGCGTACCGGTTGGGATTGCATCCAGGCCGACAAGCACGAA‑GCCCGTCATCGccca                                           >  7:232975/1‑148 (MQ=255)                                                                                                                     gCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTc                                  <  4:195895/149‑1 (MQ=255)                                                                                                                         cTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTgg                              <  1:239925/149‑1 (MQ=255)                                                                                                                            cACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGgcg                            >  6:327271/1‑148 (MQ=255)                                                                                                                                gTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCt                       >  5:35296/1‑149 (MQ=255)                                                                                                                                   cAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGGTCCACCAACGCTGTGGGTTCCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTg                      >  6:287393/1‑147 (MQ=255)                                                                                                                                           tGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAgcaggcag            >  1:359623/1‑149 (MQ=255)                                                                                                                                            gCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGc           <  3:81055/149‑1 (MQ=255)                                                                                                                                             cc‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATCGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGc           <  1:280829/149‑1 (MQ=255)                                                                                                                                             cc‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCg          <  2:231193/149‑1 (MQ=255)                                                                                                                                             cc‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCg          <  5:327271/149‑1 (MQ=255)                                                                                                                                              c‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCgg         >  6:134092/1‑149 (MQ=255)                                                                                                                                              c‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCgg         >  3:116310/1‑149 (MQ=255)                                                                                                                                                  ggCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGc        <  2:359623/147‑1 (MQ=255)                                                                                                                                                      gggCTGGCCAACGTCACGCTGGTGTGCTGCCTGCCATGATCCACCAACGCTGTGCGTACCGGTTTGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTTATTGCTCTTGCAACTTTGTGGTGCATCTGAAGTATGTAGCGGCCTGTGa  <  8:57737/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCC‑TTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGAT‑GCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCTGTGA  >  NC_002947/2421420‑2421716

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