Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,449,297	G→T	8.4%	R165S (AGG→AGT) ‡	PP_3067 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,449,297	0	G	T	8.4%	214.0 / 8.5	72	R165S (AGG→AGT) ‡	PP_3067	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base T (4/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

TTTGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAACACAAAGGTTCGCCAA  >  NC_002947/3449150‑3449445                                                                                                                                                    |                                                                                                                                                      tttGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAgg                                                                                                                                                       <  3:358917/148‑1 (MQ=255)    gggggAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGAc                                                                                                                                                    >  7:53305/1‑148 (MQ=255)    gggggAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGAc                                                                                                                                                    <  3:162979/148‑1 (MQ=255)    gggggAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCAATCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGAc                                                                                                                                                    <  3:331847/148‑1 (MQ=255)                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                                                                                                <  5:34361/149‑1 (MQ=255)                                                       ttCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAaga                                                                                                <  7:20706/149‑1 (MQ=255)                                                        tCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGgtgggtgg                                                                                                                 <  4:356587/131‑1 (MQ=255)                                                          aCGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGc                                                                                             <  5:26608/149‑1 (MQ=255)                                                           cGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCCCTGACGCTGTCGCtc                                                                                                                                          >  3:287808/1‑102 (MQ=255)                                                           cGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCCCTGACGCTGTCGCtc                                                                                                                                          <  4:287808/103‑2 (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                                                        >  3:118389/1‑125 (MQ=255)                                                                                                                       gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGTTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                                        <  4:118389/125‑1 (MQ=255)                                                                                                                       gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATAc                                                          >  6:101315/1‑123 (MQ=255)                                                                                                                       gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATAc                                                          <  5:101315/123‑1 (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               >  1:10216/1‑148 (MQ=255)                                                                                                                                        cGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAacac               >  8:5524/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                     >  4:203495/1‑135 (MQ=255)                                                                                                                                               gTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGcaaac                     <  3:203495/135‑1 (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)                                                                                                                                                    |                                                                                                                                                      TTTGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCT‑CCTGGATTTGAAGATGCAAACAAACACAAAGGTTCGCCAA  >  NC_002947/3449150‑3449445

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