Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,449,304	C→G	8.4%	L168V (CTC→GTC)	PP_3067 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,449,304	0	C	G	8.4%	216.5 / 9.4	71	L168V (CTC→GTC)	PP_3067	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base G (4/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

ATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAACACAAAGGTTCGCCAATTTT  >  NC_002947/3449165‑3449449                                                                                                                                            |                                                                                                                                                   aTCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGAGAGCGTCCCTGACGCTGTCGCtcacc                                                                                                                                           <  5:112415/149‑5 (MQ=255)      ttACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTc                                                                                                                                      <  8:53305/149‑1 (MQ=255)          tAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAgg                                                                                                                                  >  1:94696/1‑149 (MQ=255)            aaCATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGa                                                                                                                                    >  1:131807/1‑145 (MQ=255)            aaCATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGa                                                                                                                                    <  2:131807/145‑1 (MQ=255)               aTCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAg                                                                                                                             <  2:94696/149‑1 (MQ=255)                 cccTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGATCGAGGGGAAggt                                                                                                                           >  2:192311/1‑149 (MQ=255)                      ccTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAggtggg                                                                                                                        <  5:29658/147‑1 (MQ=255)                            tAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCTCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGCTGGGTGGGGaa                                                                                                                 >  1:49787/1‑148 (MQ=255)                              aGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAatat                                                                                                              <  2:232030/149‑1 (MQ=255)                                   ttGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTTTGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAggtgggtg                                                                                                                      >  1:360771/1‑136 (MQ=255)                                   ttGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAggtgggtg                                                                                                                      <  2:360771/136‑1 (MQ=255)                                        ttCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAaga                                                                                                    <  7:20706/149‑1 (MQ=255)                                        ttCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAaga                                                                                                    <  5:34361/149‑1 (MQ=255)                                         tCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGgtgggtgg                                                                                                                     <  4:356587/131‑1 (MQ=255)                                           aCGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGc                                                                                                 <  5:26608/149‑1 (MQ=255)                                            cGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCCCTGACGCTGTCGCtc                                                                                                                                              <  4:287808/103‑2 (MQ=255)                                            cGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCCCTGACGCTGTCGCtc                                                                                                                                              >  3:287808/1‑102 (MQ=255)                                              ttttCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTtg                                                                                              >  4:192199/1‑149 (MQ=255)                                                 tcttcAGGGGCCCAACAGGGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTtgat                                                                                            <  4:57945/148‑1 (MQ=255)                                                        gtgtCCAACAGTGTCATTGATTAACTTCAAGCGAAAGGTGGCGCTTATGATGAAATGGTGGTGAGCGAGAGCGTCAGGGACGCTCTCGCTGAATTTCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATGAt                                                                                    <  4:103049/149‑1 (MQ=255)                                                           tcgaaCAGTGTCATTGAGGAACTTCAAGCGAAAGGTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTaa                                                                                 <  2:49787/146‑1 (MQ=255)                                                           tCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGAGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTa                                                                                  <  4:79998/148‑1 (MQ=255)                                                              aaCAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTtgat                                                                                            <  5:244168/135‑1 (MQ=255)                                                              aaCAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTtgat                                                                                            >  6:244168/1‑135 (MQ=255)                                                               actgtgtCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATTATGAAATGGTTGTTAGCGACAGCGTCAGGGACGCTCTCGCTGAATTTCTCGAGGTGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTaaggaa                                                                             <  5:348893/146‑1 (MQ=255)                                                                 aGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCAtt                                                                                   >  7:202460/1‑141 (MQ=255)                                                                 aGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCAtt                                                                                   <  8:202460/141‑1 (MQ=255)                                                                  gtgtCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAAAGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGCGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAggg                                                                          >  7:186759/1‑149 (MQ=255)                                                                          gAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTc                                                                     >  5:207319/1‑146 (MQ=255)                                                                           aGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTc                                                                     >  5:77762/1‑145 (MQ=255)                                                                           aGGAACTTCAAGCGAAAGCTGCCGATTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTc                                                                     <  6:77762/145‑1 (MQ=255)                                                                            ggAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGAt                                                                <  3:192199/149‑1 (MQ=255)                                                                               aCTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACg                                                             >  6:243108/1‑149 (MQ=255)                                                                                cTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACg                                                             >  6:325830/1‑148 (MQ=255)                                                                                 ttCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGAt                                                                <  2:216016/144‑1 (MQ=255)                                                                                 ttCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGAt                                                                >  1:216016/1‑144 (MQ=255)                                                                                 ttCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                                            >  5:168201/1‑148 (MQ=255)                                                                                      gCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCa                                                      >  1:166136/1‑149 (MQ=255)                                                                                         aaaGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCatat                                                   >  6:322357/1‑149 (MQ=255)                                                                                          aaGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGa                                                                        <  8:279584/127‑1 (MQ=255)                                                                                          aaGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGa                                                                        >  7:279584/1‑127 (MQ=255)                                                                                           aGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATtc                                                 <  2:166136/149‑1 (MQ=255)                                                                                              tGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATtcct‑c                                             >  5:93275/1‑149 (MQ=255)                                                                                               gCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCa                                                                                     <  4:368111/109‑1 (MQ=255)                                                                                               gCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCa                                                                                     >  3:368111/1‑109 (MQ=255)                                                                                                ccGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATtcct‑c                                             <  5:322357/147‑1 (MQ=255)                                                                                                    ttATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGAt                                       >  2:236244/1‑149 (MQ=255)                                                                                                        gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGTTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                                            <  4:118389/125‑1 (MQ=255)                                                                                                        gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGTTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                                            >  3:118389/1‑125 (MQ=255)                                                                                                        gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATAc                                                              <  5:101315/123‑1 (MQ=255)                                                                                                        gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATAc                                                              >  6:101315/1‑123 (MQ=255)                                                                                                         atgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGaa                                  <  1:192311/149‑1 (MQ=255)                                                                                                          tgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTg                                    >  3:239327/1‑146 (MQ=255)                                                                                                           gaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTg                                    >  2:164645/1‑145 (MQ=255)                                                                                                               tggtggTGAGCGACATCTTCAGGGACGCTCTCGCTTAATTGCTCGAGGGGAAGGTGGGTGGGGAATATTATAAACAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGca                            <  4:239327/149‑1 (MQ=255)                                                                                                                  tggtgAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAATTCATATTCCT‑CATGGATGTGAAGAGGaaaa                          >  4:248954/1‑148 (MQ=255)                                                                                                                     tgAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGcaaacaaa                      <  1:164645/149‑1 (MQ=255)                                                                                                                        gCGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCAT‑CCTGGATTTGAAGATGCAaacaaacaa                   >  7:329372/1‑148 (MQ=255)                                                                                                                         cGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAacac                   >  8:5524/1‑148 (MQ=255)                                                                                                                         cGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAacac                   >  1:10216/1‑148 (MQ=255)                                                                                                                         cGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAacac                   >  5:60692/1‑148 (MQ=255)                                                                                                                         cGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAacac                   >  4:48555/1‑148 (MQ=255)                                                                                                                          gACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAACACaa                 <  6:93275/149‑1 (MQ=255)                                                                                                                                gTCAGGGACGCTCTCGCTGAATTTCTCGAGGTGAAGGTGGGTTTTTAATATGATAAAAAGAGGCTTTATGAAATCATTAAGGAAGGTGAGTATAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATTCAAACAATCACAAAGTTTc           <  6:168201/149‑1 (MQ=255)                                                                                                                                gTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGcaaac                         <  3:203495/135‑1 (MQ=255)                                                                                                                                gTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGcaaac                         >  4:203495/1‑135 (MQ=255)                                                                                                                                 tCATGGAAGCTCTTGCTTTATTTCTAGATGGGAATGTGGGTGGGGAATATGATAAAAATAGGCTTTATTAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCCTGGATTTGAAGATGCAAACAAACACAAAGGTTc           <  1:236244/149‑1 (MQ=255)                                                                                                                                 tCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGGTTTGAAGATGCAAACAAACACAAAGGTTCg          >  8:325355/1‑149 (MQ=255)                                                                                                                                     ggACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAACACAAAGGTTCGCCaa      >  1:252362/1‑149 (MQ=255)                                                                                                                                         gCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGAATTGAAGATCCAAACAAAAACAAACGTACGCCAAtttt  >  4:78767/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                   ATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAACACAAAGGTTCGCCAATTTT  >  NC_002947/3449165‑3449449

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