Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,338,404	T→A	5.9%	noncoding (128/395 nt)	PP_mr53 →	CrcZ ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,338,404	0	T	A	5.9%	207.8 / 2.8	68	noncoding (128/395 nt)	PP_mr53	CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base A (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

ACGGCACCTGCTATATGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGA  >  NC_002947/5338256‑5338548                                                                                                                                                     |                                                                                                                                                  accgcacCTGCTTAATGATTGGTACAAGAAAAATAACAAGCACTGCAACTAAGAATAATAACAATACGAAACGGCTCACGCACAATAAAAACAACATGGGTGAGGCTAAGCAAACTGATTCATTTGGAGAGGAAGCGTGTTTGGGGCtt                                                                                                                                                   >  2:118869/4‑149 (MQ=255)  cggcacCTGCTATATGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑g                                                                                                                                                 >  3:2718/1‑149 (MQ=255)    gcacCTGCTATATGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑Gcc                                                                                                                                               <  3:59096/149‑1 (MQ=255)         tGCTATATGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCcac                                                                                                                                           >  7:227093/1‑148 (MQ=255)           cTATATGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTGTGGGGCAA‑Gcccaaacac                                                                                                                                        <  7:175589/149‑7 (MQ=255)            tataTGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAAc                                                                                                                                        >  4:187338/1‑148 (MQ=255)              tatGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAAcc                                                                                                                                       <  3:120949/147‑1 (MQ=255)              tatGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAg                                                                                                                                     >  7:99011/1‑149 (MQ=255)              tatGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAg                                                                                                                                     >  1:305395/1‑149 (MQ=255)                  tttGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAgg                                                                                                                                    >  2:109295/1‑146 (MQ=255)                  tttGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAAAAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAgg                                                                                                                                    <  1:109295/146‑1 (MQ=255)                    tgctACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGTGGCTT‑GCCCCACAACCAGGCaga                                                                                                                                <  4:106943/145‑1 (MQ=255)                       tACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAAGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGaa                                                                                                                             >  3:311628/1‑148 (MQ=255)                        aCAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTCGGGGCTT‑GCCCCACAACCAGGCAGAGaac                                                                                                                            >  5:226248/1‑148 (MQ=255)                           aGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGaacaac                                                                                                                         >  4:286087/1‑148 (MQ=255)                            gAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAacaaca                                                                                                                        >  1:205016/1‑148 (MQ=255)                              aCAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACaaa                                                                                                                      >  1:214126/1‑148 (MQ=255)                                aaTAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTCGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAAc                                                                                                                   <  6:226248/149‑1 (MQ=255)                                 aTAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACt                                                                                                                  <  2:205016/149‑1 (MQ=255)                                    aCAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCa                                                                                                               <  2:305395/149‑1 (MQ=255)                                    aCAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCa                                                                                                               <  8:227093/149‑1 (MQ=255)                                     cAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCAc                                                                                                              >  3:88256/1‑149 (MQ=255)                                      aaGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCAc                                                                                                              >  8:204542/1‑148 (MQ=255)                                          aCTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTaa                                                                                                           >  3:162208/1‑147 (MQ=255)                                               aaCTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAagcagc                                                                                                    <  4:7064/149‑1 (MQ=255)                                                 cTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAgcgc                                                                                                  <  3:7486/149‑1 (MQ=255)                                                      aaTAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGaa                                                                                             >  2:42403/1‑149 (MQ=255)                                                       aTAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCAA‑GCCCcaaac                                                                                                                                         >  4:179397/1‑101 (MQ=255)                                                       aTAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCAA‑GCCCcaaac                                                                                                                                         <  3:179397/104‑4 (MQ=255)                                                          agaagaagaCGACACGGCTCACGCACAATATAAACAAGATGGCGGAGGCGCAGCTAACTCATTCTTTTGGAGTTGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGATCAACAAAAACTGCACTAAAAGAAGCGCCTGAACTg                                                                                          <  6:175932/148‑1 (MQ=255)                                                           gAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGaacaac                                                                                                                         >  5:203408/1‑116 (MQ=255)                                                           gAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGaacaac                                                                                                                         <  6:203408/116‑1 (MQ=255)                                                           gAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGt                                                                                        <  7:34290/149‑1 (MQ=255)                                                            aaCAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAAGCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGtt                                                                                       <  1:81968/149‑1 (MQ=255)                                                            aaCAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGtt                                                                                       <  4:162208/149‑1 (MQ=255)                                                             aCAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCTACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTg                                                                                      <  1:268074/149‑1 (MQ=255)                                                             aCAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTg                                                                                      >  3:327175/1‑149 (MQ=255)                                                             aCAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTg                                                                                      <  2:214126/149‑1 (MQ=255)                                                               aaGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGa                                                                                    >  7:171496/1‑149 (MQ=255)                                                                aGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGAt                                                                                   >  7:43813/1‑149 (MQ=255)                                                                 gACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATc                                                                                  >  7:200912/1‑149 (MQ=255)                                                                 gACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATc                                                                                  <  3:286087/149‑1 (MQ=255)                                                                   cGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTg                                                                                      >  7:212531/1‑143 (MQ=255)                                                                   cGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTg                                                                                      <  8:212531/143‑1 (MQ=255)                                                                     aaaCGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACa                                                                              <  3:88912/149‑1 (MQ=255)                                                                         ggCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATc                                                                          <  8:43813/149‑1 (MQ=255)                                                                              aCGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACg                                                                      <  2:343155/148‑1 (MQ=255)                                                                                  acaATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACa                                                                 >  5:300362/1‑149 (MQ=255)                                                                                                 tGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGATCAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCa                                                  <  3:187338/149‑1 (MQ=255)                                                                                                    cGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATc                                               <  6:262583/149‑1 (MQ=255)                                                                                                     ggAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATcc                                              <  8:200912/149‑1 (MQ=255)                                                                                                         gcgcAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGttt                                          <  4:88256/149‑1 (MQ=255)                                                                                                             aGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGctc                                      >  6:280107/1‑149 (MQ=255)                                                                                                              gCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGctct                                     <  4:327175/149‑1 (MQ=255)                                                                                                                 aaCTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTa                                  >  5:186645/1‑149 (MQ=255)                                                                                                                  aCTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTAc                                 >  4:185343/1‑149 (MQ=255)                                                                                                                   cTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTAcc                                <  4:2718/149‑1 (MQ=255)                                                                                                                    tGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTAccc                               >  2:159015/1‑149 (MQ=255)                                                                                                                     gATTCTTTTGGAGAGGATGCGTGTTTGGGGCTTGGCCCAACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGATCTTTAccc                               >  1:221483/1‑149 (MQ=255)                                                                                                                      aTTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATGAGCGGCCAAAGCAATCCGTTTGCTCTTTAccc                               <  1:42403/147‑1 (MQ=255)                                                                                                                      aTTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTAcccc                              <  2:221483/148‑1 (MQ=255)                                                                                                                       ttCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCg                            <  6:186645/149‑1 (MQ=255)                                                                                                                            ttGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGt                          <  5:280107/146‑1 (MQ=255)                                                                                                                                  agGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAg                                                    <  6:336534/114‑1 (MQ=255)                                                                                                                                  agGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAg                                                    >  5:336534/1‑114 (MQ=255)                                                                                                                                         gtgtTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGAACTAAAAGCAGCTCCTGAACTGGTAGGATCGACAGATCAACTTGACATCAGCGGCCAAAGCAATCCGTTAGCTCTTTACCCCAGGTTTGGGGGTAGTCaac           >  8:150739/1‑148 (MQ=255)                                                                                                                                          gttgtGGGGCAA‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGTCATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCcac           >  4:274984/5‑147 (MQ=255)                                                                                                                                                 ggCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCga  <  8:99011/149‑1 (MQ=255)                                                                                                                                                  gCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTc    >  8:183311/1‑146 (MQ=255)                                                                                                                                                     |                                                                                                                                                  ACGGCACCTGCTATATGTTTGGTACAAGAACAATAACAAGCACTGCAACTCAGAATAAGAACAAGACGAAACGGCTCACGCACAATAAAAACAAGATGGCGGAGGCGCAGCTAACTGATTCTTTTGGAGAGGATGCGTGTTTGGGGCTT‑GCCCCACAACCAGGCAGAGAACAACAAAAACTGCACTAAAAGCAGCGCCTGAACTGGTTGGATCGACAGATCAACGTGACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGA  >  NC_002947/5338256‑5338548

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