Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	902,619	A→C	7.0%	intergenic (+45/+31)	trxB → / ← nadC	thioredoxin reductase/quinolinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	902,619	0	A	C	7.0%	194.4 / 4.5	68	intergenic (+45/+31)	trxB/nadC	thioredoxin reductase/quinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/35);  new base C (2/3);  total (30/38)
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.56e-01

GACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTT  >  NC_002947/902473‑902760                                                                                                                                                   |                                                                                                                                              gACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAg                                                                                                                                             >  6:90379/1‑149 (MQ=255)  aCCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAgc                                                                                                                                            <  6:222577/149‑1 (MQ=255)    cACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAgcgc                                                                                                                                          >  7:90667/1‑149 (MQ=255)      cGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCg                                                                                                                                        <  7:279056/149‑1 (MQ=255)      cGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCg                                                                                                                                        <  2:3189/149‑1 (MQ=255)       gTTTACCGCCAGGCCATCACTACGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCAAAAAAAAAACCGCCAAACCCCCgg                                                                                                                                       >  2:322540/1‑149 (MQ=255)           aCCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                    <  2:25757/148‑1 (MQ=255)            ccGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGGCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttttt                                                                                                                                  <  4:69032/149‑1 (MQ=255)             cGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGtttttt                                                                                                                                 <  4:68750/149‑1 (MQ=255)                cAGGCCATCACTTCGGCCGGTGCCGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGAGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTAtt                                                                                                                              <  7:273837/149‑1 (MQ=255)                  ggCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGTTTTTTATTgg                                                                                                                            <  3:1568/149‑2 (MQ=39)                  ggCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGTTTTTTATTgg                                                                                                                            <  7:277636/149‑2 (MQ=39)                    ccATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCt                                                                                                                          <  6:13159/149‑1 (MQ=255)                     cATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGACGGTTGTTTATTGCCTg                                                                                                                         >  6:300240/1‑149 (MQ=255)                      aTCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGtttt                                                                                                                                   >  1:276383/1‑138 (MQ=255)                      aTCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGtttt                                                                                                                                   <  2:276383/138‑1 (MQ=255)                          cTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAAt                                                                                                                     >  7:309206/1‑148 (MQ=255)                           ttCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAAt                                                                                                                     >  7:28198/1‑147 (MQ=255)                             cGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAg                                                                                                                  <  5:90379/148‑1 (MQ=255)                              ggCCGGTGCGGGCTGCATGGCGGCGCTTGATGTCGAGCGCTACCTGGAAGGCTTGGCCAACGCCTCGTGCTGATCTGAGAAAAAGGGGGGCAAAGCGCCCCCAATAAAAAACCCGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCt                                                                                                                <  3:224867/149‑1 (MQ=255)                               gCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                    <  3:351359/128‑1 (MQ=255)                               gCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                    >  4:351359/1‑128 (MQ=255)                                 cGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTcac                                                                                                             <  1:186167/149‑1 (MQ=255)                                  ggTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCaa                                                                                                                                              >  7:42320/1‑115 (MQ=255)                                  ggTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCaa                                                                                                                                              <  8:42320/115‑1 (MQ=255)                                  ggTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTcaca                                                                                                            <  8:90667/149‑1 (MQ=255)                                  ggTGC‑GGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAg                                                                                                           <  7:67802/149‑1 (MQ=255)                                   gTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAg                                                                                                           >  3:169133/1‑149 (MQ=255)                                           tGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTc                                                                                                   <  2:142021/149‑1 (MQ=255)                                             cATGGCGGCGCTGGATGTCGCTCGCTACCTGGACGTTTTGGCCTACGCATACTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCGCCAATAATAAACAGGCCAAGCGCCTGTTATTTATTGCCTGCAATCAGCTCACAGGCGCAGTc                                                                                                   <  4:26612/147‑1 (MQ=255)                                             cATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGc                                                                                                 <  4:83189/149‑1 (MQ=255)                                             cATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTGTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCg                                                                                                  >  6:181932/1‑148 (MQ=255)                                                  cggcgCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCTTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGa                                                                                            <  3:201703/149‑1 (MQ=255)                                                      gcTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAgg                                                                                        >  5:265412/1‑149 (MQ=255)                                                      gcTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGAAATCAGCTCACAGGCTCAGTCGCATCGACAgg                                                                                        >  4:217300/1‑149 (MQ=255)                                                       cTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGt                                                                                       <  3:377002/149‑1 (MQ=255)                                                           aTGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCAc                                                                                   <  7:233926/149‑1 (MQ=255)                                                                     cTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTcac                                                                          >  8:104458/1‑148 (MQ=255)                                                                      tACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACAt                                                                        >  6:299449/1‑149 (MQ=255)                                                                           ggACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTg                                                                   >  5:176880/1‑149 (MQ=255)                                                                            gACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTgg                                                                  >  4:83411/1‑149 (MQ=255)                                                                            gACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTg                                                                   <  6:345304/148‑1 (MQ=255)                                                                              cGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTc                                                                <  8:28198/149‑1 (MQ=255)                                                                                gTTTGGCCAACGGCTTGGTCTGATCTGATAAAAACGGGGGTAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCAt                                                              <  1:322540/149‑1 (MQ=255)                                                                                       cAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGTTTTTTAtt                                                                                                                              >  5:177573/1‑78 (MQ=25)                                                                                       cAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGTTTTTTAtt                                                                                                                              <  6:177573/78‑1 (MQ=25)                                                                                         aCGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAAt                                                     >  2:214656/1‑149 (MQ=255)                                                                                            ccTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCga                                                  <  7:104458/149‑1 (MQ=255)                                                                                               cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           >  6:201355/1‑105 (MQ=255)                                                                                               cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           <  1:312367/105‑1 (MQ=255)                                                                                               cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           <  5:201355/105‑1 (MQ=255)                                                                                               cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           >  2:312367/1‑105 (MQ=255)                                                                                                gTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCgaga                                                <  5:97497/147‑1 (MQ=255)                                                                                                gTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGAt                                               <  5:299449/148‑1 (MQ=255)                                                                                                 ttCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATg                                              <  4:169133/148‑1 (MQ=255)                                                                                                    tgatctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGACACCAATCGAGATGTAGt                                          <  3:83411/149‑1 (MQ=255)                                                                                                     gatctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAAGCGAGATGTAGTc                                         >  6:223363/1‑149 (MQ=255)                                                                                                      atctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTc                                         >  2:287133/1‑148 (MQ=255)                                                                                                        ctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTTGTCCAc                                      >  1:282498/1‑149 (MQ=255)                                                                                                         tgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCAcc                                     >  8:152486/1‑149 (MQ=255)                                                                                                           agaAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATAAGATCACAGGCTCAGTCGCATCGAGAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCAcccc                                   >  5:274853/1‑149 (MQ=255)                                                                                                                      tGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCaa                        >  2:382275/1‑149 (MQ=255)                                                                                                                            gcgcCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCAccc                  <  7:152486/149‑1 (MQ=255)                                                                                                                                ccccAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAg              >  2:278035/1‑149 (MQ=255)                                                                                                                                aggcAATAAAAAACCGGCGCTTGGCCGGTTTGTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAg              >  1:380344/4‑149 (MQ=39)                                                                                                                                        aaaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgttt      >  4:43669/1‑149 (MQ=255)                                                                                                                                         aaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttc     <  5:97449/149‑1 (MQ=255)                                                                                                                                          aaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttcg    >  6:20695/1‑149 (MQ=255)                                                                                                                                            aCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCATGGtcgtttcgtt  <  3:153050/149‑1 (MQ=255)                                                                                                                                            aCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttcgtt  <  1:382275/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              GACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTT  >  NC_002947/902473‑902760

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