Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,490,825	A→T	5.1%	Y214N (TAC→AAC)	PP_5715 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,490,825	0	A	T	5.1%	232.0 / 2.3	79	Y214N (TAC→AAC)	PP_5715	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base T (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

AACAGGATGTCGTCACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCTGCCACAGCCAATCCATCGCCA  >  NC_002947/5490679‑5490971                                                                                                                                                   |                                                                                                                                                   aaCAGGATGTCGTCACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACt                                                                                                                                                  >  7:167695/1‑149 (MQ=255)      gATGTCGTCACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAAt                                                                                                                                             >  8:70026/1‑149 (MQ=255)      gATGTCGTCACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAAt                                                                                                                                             >  8:204083/1‑149 (MQ=255)      gATCTCGTCACGGGCTTCGCGACTAAAACCTGGCAAGAATCACAGAAGGCTGCTGTGGCCGAGTTTGCCACTGGCGGCACATTGATTGAACTCACTCAATTTCGACTCCTGAGGTAGTGGCAAGTTAGACCTTGCTGCGTTACTGCAAt                                                                                                                                             >  8:204085/1‑149 (MQ=255)       aTGTCGTCACGGGCTTCGCTTCTACCACCTGTCAAGAATGTCATCACGCGTCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCCCTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCTAtt                                                                                                                                            <  7:204085/149‑1 (MQ=255)       aTGTCGTCACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAAtt                                                                                                                                            >  8:73502/1‑149 (MQ=255)       aTGTCGTCACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAAtt                                                                                                                                            <  7:204083/149‑1 (MQ=255)           cgtcACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCtt                                                                                                                                        >  1:239355/1‑149 (MQ=255)            gtcACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTg                                                                                                                                       <  4:96835/149‑1 (MQ=255)             tcACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTGCTGAGGTAGTGGAAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGa                                                                                                                                      >  2:267315/1‑149 (MQ=255)              cACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGAc                                                                                                                                     >  4:281744/1‑149 (MQ=255)               aCGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGAc                                                                                                                                     <  7:278040/148‑1 (MQ=255)                    cTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGAt                                                                                                                               <  7:305363/149‑1 (MQ=255)                        gcgACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGc                                                                                                                            <  2:239355/148‑1 (MQ=255)                           aCTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCa                                                                                                                        >  1:6706/1‑149 (MQ=255)                                 acCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCTACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAg                                                                                                                  <  4:64529/149‑1 (MQ=255)                                 acCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAg                                                                                                                  >  3:306764/1‑149 (MQ=255)                                   cTGGAAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAAATCACTCATTTTCGACTCCTGAGGTAGTTGCAATTTTGACATGGCTGCGGTACTGGAATTGGTTGACATTGATAGCGTCATGCAATcc                                                                                                                >  7:303480/1‑149 (MQ=255)                                     ggCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACATTGATTGAACTCACTAATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATTGCGTCATGCAAGCCCg                                                                                                              >  2:88867/1‑149 (MQ=255)                                      gCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGGAAGCCCg                                                                                                              <  8:167695/148‑1 (MQ=255)                                      gCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGccc                                                                                                               <  3:281744/147‑1 (MQ=255)                                      gCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCg                                                                                                              >  2:30061/1‑148 (MQ=255)                                       cAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGa                                                                                                            <  7:73502/149‑1 (MQ=255)                                         aGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCt                                                                                                                                         >  7:14893/1‑118 (MQ=255)                                         aGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCt                                                                                                                                         <  8:14893/118‑1 (MQ=255)                                            aTGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCTCTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTTATAGGGTCATGCAAGCCCGGACtaa                                                                                                        >  8:134156/1‑146 (MQ=255)                                               aCAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCt                                                                                                                                         <  5:256542/112‑1 (MQ=255)                                               aCAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCt                                                                                                                                         >  6:256542/1‑112 (MQ=255)                                                 agcagGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAAc                                                                                                  <  1:305433/149‑1 (MQ=255)                                                          ctgTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACt                                                                                          <  2:6706/148‑1 (MQ=255)                                                             tgGACGAGTTTGCCACTGGCGGCAGCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACttttt                                                                                      <  2:217480/149‑1 (MQ=255)                                                              ggACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACttttt                                                                                      >  5:222176/1‑148 (MQ=255)                                                                  gAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATgg                                                                                 >  2:23142/1‑149 (MQ=255)                                                                  gAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATgg                                                                                 <  7:70026/149‑1 (MQ=255)                                                                         ccACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGaaa                                                                           >  5:300747/1‑148 (MQ=255)                                                                         ccACTGGCGGCACCATGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAaat                                                                          <  8:167652/149‑1 (MQ=255)                                                                          cACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTTCTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAaat                                                                          >  8:220783/1‑148 (MQ=255)                                                                          cACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTTCTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAaat                                                                          >  8:225034/1‑148 (MQ=255)                                                                          cACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAaat                                                                          >  7:137670/1‑148 (MQ=255)                                                                               gcggcACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTgg                                                                    >  2:310641/1‑149 (MQ=255)                                                                               gcggcACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTg                                                                     >  3:162152/1‑148 (MQ=255)                                                                                   cACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATc                                                                <  4:306764/149‑1 (MQ=255)                                                                                   cACCTTGATTGAACTCACTCATTTTAGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATc                                                                <  3:203970/149‑1 (MQ=255)                                                                                      cTTGATTGAACTCACTGATTTTCGACTCCTGATGGAGTTTCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATATCGTCATGCAAGCCCTGACTTATCGAACTTATCACTTTTTCATGGTTGAAATAATTTGCATCGaa                                                             <  8:303480/149‑1 (MQ=255)                                                                                      cTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGaa                                                             >  1:131495/1‑149 (MQ=255)                                                                                       ttgattgaACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAAATGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGCTGAAATAATTGGCATCGAAc                                                            >  5:206949/1‑149 (MQ=255)                                                                                                ctcactcaTTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTg                                                   >  5:296800/1‑149 (MQ=255)                                                                                                 tcactcaTTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTgg                                                  >  5:301074/1‑149 (MQ=255)                                                                                                   actcaTTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCAt                                                                 >  8:244862/1‑132 (MQ=255)                                                                                                   actcaTTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCAt                                                                 <  7:244862/132‑1 (MQ=255)                                                                                                   actcaTTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGcc                                                >  8:132822/1‑149 (MQ=255)                                                                                                       aTTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTggg                                            >  1:29015/1‑149 (MQ=255)                                                                                                       aTTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTggg                                            <  6:301074/149‑1 (MQ=255)                                                                                                       aTTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTg                                              <  1:267315/147‑1 (MQ=255)                                                                                                           tCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTTCTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCAc                                         <  7:225034/148‑1 (MQ=255)                                                                                                           tCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTTCTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCAc                                         <  7:220783/148‑1 (MQ=255)                                                                                                               cTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGaa                                    <  1:310641/149‑1 (MQ=255)                                                                                                                tCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGaaa                                   <  4:162152/149‑1 (MQ=255)                                                                                                                tCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGa                                     >  1:277544/1‑147 (MQ=255)                                                                                                                  cTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGAt                                                                                              <  5:115234/88‑1 (MQ=255)                                                                                                                  cTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGAt                                                                                              >  6:115234/1‑88 (MQ=255)                                                                                                                  cTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAg                                  <  2:29015/148‑1 (MQ=255)                                                                                                                   tGAGGTAGTGGCAAGTTTGACCTGGCTGAGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGACCTTCTGCGTGGCCTGGGCACTTGAAAgct                                >  8:150247/1‑147 (MQ=255)                                                                                                                     aGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCAt                                                                                   >  7:45019/1‑96 (MQ=255)                                                                                                                     aGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCAt                                                                                   <  8:45019/96‑1 (MQ=255)                                                                                                                         aGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCg                          >  7:223876/1‑149 (MQ=255)                                                                                                                            ggCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATTCAAGCCCGGACTTATAGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGc                        <  7:134156/148‑1 (MQ=255)                                                                                                                             gCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCt                       <  8:137670/148‑1 (MQ=255)                                                                                                                             gCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCTg                      >  4:312318/1‑149 (MQ=255)                                                                                                                                   ttGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCgg                         >  7:270420/1‑140 (MQ=255)                                                                                                                                   ttGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCgg                         <  8:270420/140‑1 (MQ=255)                                                                                                                                    tGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCTGCCACAg                <  6:296800/148‑1 (MQ=255)                                                                                                                                      aCCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCTGCCACAGcc              >  2:139517/1‑148 (MQ=255)                                                                                                                                          ggCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACt                                        >  6:203987/1‑118 (MQ=255)                                                                                                                                          ggCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACt                                        <  5:203987/118‑1 (MQ=255)                                                                                                                                           gCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCTGCCACAGCCAATc          <  8:223876/147‑1 (MQ=255)                                                                                                                                              gCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCTGCCACAGCCATTccctag     >  8:289334/1‑145 (MQ=255)                                                                                                                                                 gTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCTGCCACAGCCAATCCATCGc    >  8:77394/1‑147 (MQ=255)                                                                                                                                                 gTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCTGCCACAGCCAATCCATCGCCa  <  2:277544/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AACAGGATGTCGTCACGGGCTTCGCGACTAACACCTGGCAAGAATGACAGCAGGCTGCTGTGGACGAGTTTGCCACTGGCGGCACCTTGATTGAACTCACTCATTTTCGACTCCTGAGGTAGTGGCAAGTTTGACCTGGCTGCGGTACTGCAATTGCTTGACCTTGATAGCGTCATGCAAGCCCGGACTTATCGAACTGATCACTTTTTCATGGGTGAAATAATTGGCATCGAACTGCTGAGTGGCCTGGGCACTTGAAAGATTCAGCGGCTGCCACAGCCAATCCATCGCCA  >  NC_002947/5490679‑5490971

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