Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,449,288	C→G	6.1%	D162E (GAC→GAG)	PP_3067 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,449,288	0	C	G	6.1%	236.7 / 3.4	76	D162E (GAC→GAG)	PP_3067	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/37);  new base G (3/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

CGTCTCTTTTTTGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAACACAA  >  NC_002947/3449141‑3449434                                                                                                                                                    |                                                                                                                                                   agtgtcgtttttGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACa                                                                                                                                                   <  5:24319/142‑1 (MQ=255)  gTCTCTTTTTTGGGGTAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATTAAATGGTGGTGAGCGAc                                                                                                                                                    <  2:76816/147‑1 (MQ=255)  gTCTCTTTTTTGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGTTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAg                                                                                                                                                  >  3:103049/1‑149 (MQ=255)   tctcTTTTTTGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGc                                                                                                                                                 >  6:29658/1‑149 (MQ=255)       ttttttGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGt                                                                                                                                               <  6:131412/147‑1 (MQ=255)          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)                             aTTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTc                                                                                                                                              <  3:286420/126‑1 (MQ=255)                             aTTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTc                                                                                                                                              >  4:286420/1‑126 (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)                                                tCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTc                                                                                                                                              <  2:180155/107‑1 (MQ=255)                                                tCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTc                                                                                                                                              >  1:180155/1‑107 (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                                                                                                                              <  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                                                                   <  8:202460/141‑1 (MQ=255)                                                                                         aGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCAtt                                                                   >  7:202460/1‑141 (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                                                >  1:216016/1‑144 (MQ=255)                                                                                                         ttCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGAt                                                <  2:216016/144‑1 (MQ=255)                                                                                                         ttCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                            >  5:168201/1‑148 (MQ=255)                                                                                                              gCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCa                                      >  1:166136/1‑149 (MQ=255)                                                                                                                 aaaGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCatat                                   >  6:322357/1‑149 (MQ=255)                                                                                                                  aaGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGa                                                        >  7:279584/1‑127 (MQ=255)                                                                                                                  aaGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGa                                                        <  8:279584/127‑1 (MQ=255)                                                                                                                   aGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATtc                                 <  2:166136/149‑1 (MQ=255)                                                                                                                      tGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATtcctc                              >  5:93275/1‑149 (MQ=255)                                                                                                                       gCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCa                                                                     >  3:368111/1‑109 (MQ=255)                                                                                                                       gCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCa                                                                     <  4:368111/109‑1 (MQ=255)                                                                                                                        ccGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATtcctc                              <  5:322357/147‑1 (MQ=255)                                                                                                                            ttATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGAt                        >  2:236244/1‑149 (MQ=255)                                                                                                                                gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGTTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                            <  4:118389/125‑1 (MQ=255)                                                                                                                                gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGTTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                            >  3:118389/1‑125 (MQ=255)                                                                                                                                gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATAc                                              >  6:101315/1‑123 (MQ=255)                                                                                                                                gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATAc                                              <  5:101315/123‑1 (MQ=255)                                                                                                                                 atgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGaa                   <  1:192311/149‑1 (MQ=255)                                                                                                                                  tgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTg                     >  3:239327/1‑146 (MQ=255)                                                                                                                                   gaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTg                     >  2:164645/1‑145 (MQ=255)                                                                                                                                       tggtggTGAGCGACATCTTCAGGGACGCTCTCGCTTAATTGCTCGAGGGGAAGGTGGGTGGGGAATATTATAAACAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGca             <  4:239327/149‑1 (MQ=255)                                                                                                                                          tggtgAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAATTCATATTCCTCATGGATGTGAAGAGGaaaa           >  4:248954/1‑148 (MQ=255)                                                                                                                                             tgAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGcaaacaaa       <  1:164645/149‑1 (MQ=255)                                                                                                                                                gCGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCATCCTGGATTTGAAGATGCAaacaaacaa    >  7:329372/1‑148 (MQ=255)                                                                                                                                                 cGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAacac    >  1:10216/1‑148 (MQ=255)                                                                                                                                                 cGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAacac    >  8:5524/1‑148 (MQ=255)                                                                                                                                                 cGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAacac    >  4:48555/1‑148 (MQ=255)                                                                                                                                                 cGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAacac    >  5:60692/1‑148 (MQ=255)                                                                                                                                                  gACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAACACaa  <  6:93275/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CGTCTCTTTTTTGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAACACAA  >  NC_002947/3449141‑3449434

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