Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	77,500	A→T	9.0%	L15H (CTC→CAC)	rsmB ←	ribosomal RNA small subunit methyltransferase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	77,500	0	A	T	9.0%	196.7 / 10.0	67	L15H (CTC→CAC)	rsmB	ribosomal RNA small subunit methyltransferase B
Reads supporting (aligned to +/- strand):  ref base A (37/24);  new base T (3/3);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

GCTTTTGCAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCT  >  NC_002947/77351‑77644                                                                                                                                                       |                                                                                                                                                 gCTTTTGCAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGa                                                                                                                                                  <  1:18157/150‑1 (MQ=255)  cTTTTGCAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAg                                                                                                                                                 >  1:268208/1‑150 (MQ=255)        cagcagCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGc                                                                                                                                           >  1:373664/1‑150 (MQ=255)         agcagcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCg                                                                                                                                          >  2:1027008/1‑150 (MQ=255)         agcagcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCg                                                                                                                                          >  2:1206978/1‑150 (MQ=255)             gcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGTGCTCAGCGGcaa                                                                                                                                      <  1:1309490/150‑3 (MQ=255)                 tgCGGCCAGCAGGTCGACGCGTGGCAGACAGCGTGTGGTACGGAACGCGAAGGGCTGGGTGAAGCCGAGGTGGCGTTCATCTTCTTTATCTAGTTGCGCCGGCCGCGAGCTGTTCAGAGAAGCCTT‑GCCGCTGAGCACAGCGGCCAgcgc                                                                                                                                  >  2:521405/1‑150 (MQ=255)                     gCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACgg                                                                                                                              <  1:676620/150‑1 (MQ=255)                      ccAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGTGCTCAGCGGCAAg                                                                                                                                     >  2:539662/1‑142 (MQ=255)                      ccAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGTGCTCAGCGGCAAg                                                                                                                                     <  1:539662/142‑1 (MQ=255)                         gcagGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGc                                                                                                                                        <  2:1177118/136‑1 (MQ=255)                         gcagGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGc                                                                                                                                        >  1:1177118/1‑136 (MQ=255)                         gcagGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGggcgg                                                                                                                          >  2:139504/1‑150 (MQ=255)                                       ggCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGTCGACGCCACGGTCATCTACCTTGTCCAGTTGAACCTGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTggg                                                                                                            <  1:1385987/150‑1 (MQ=255)                                          tGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTc                                                                                                         >  2:150515/1‑150 (MQ=255)                                           cctggCGTTCGGTACCGCACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGTGCTCAGCGGcaa                                                                                                                                      <  2:1405772/116‑3 (MQ=255)                                           cctggCGTTCGGTACCGCACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGTGCTCAGCGGcaa                                                                                                                                      >  1:1405772/5‑118 (MQ=255)                                             cAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCAtt                                                                                                      >  1:183580/1‑150 (MQ=255)                                             cAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCAtt                                                                                                      >  2:684517/1‑150 (MQ=255)                                             cAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCAtt                                                                                                      <  1:1376344/150‑1 (MQ=255)                                              aGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTg                                                                                                     <  2:268208/150‑1 (MQ=255)                                               gCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTgg                                                                                                    <  1:150515/150‑1 (MQ=255)                                                     ggTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGc                                                                                              >  2:778926/1‑150 (MQ=255)                                                       tACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGt                                                                                                               <  2:812542/131‑1 (MQ=255)                                                       tACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTACAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGt                                                                                                               >  1:812542/1‑131 (MQ=255)                                                            aaCGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGc                                                                                       >  2:1077653/1‑150 (MQ=255)                                                             aCGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGcc                                                                                      <  1:213348/150‑1 (MQ=255)                                                                  aGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCcggcgg                                                                                 <  1:684517/150‑1 (MQ=255)                                                                     tCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGa                                                                              >  2:390006/1‑150 (MQ=255)                                                                      ccTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGaa                                                                             <  1:584928/150‑1 (MQ=255)                                                                       cTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAGGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGTCCCAGCACCTTGCCGGCGGCGAAc                                                                            >  1:1203069/1‑150 (MQ=255)                                                                          ggTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTtc                                                                         >  1:171732/1‑150 (MQ=255)                                                                                gCCGCGGTCGCGTTCATCTACCTTTTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGCACTTCTCgcgg                                                                   >  1:57277/1‑150 (MQ=255)                                                                                  cgcgGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCtt                                                                                         <  1:608382/126‑1 (MQ=255)                                                                                  cgcgGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCtt                                                                                         >  2:608382/1‑126 (MQ=255)                                                                                   gcgGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGcggcgg                                                                >  1:815826/1‑150 (MQ=255)                                                                                      gTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTg                                                             >  2:1352984/1‑150 (MQ=255)                                                                                         gcgTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTg                                                          <  2:171732/150‑1 (MQ=255)                                                                                           gTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                            >  2:1089327/1‑130 (MQ=255)                                                                                           gTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                            <  1:1089327/130‑1 (MQ=255)                                                                                            tgtaTCTACCTTGTGCAGTTGCTCCTGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAAc                                                       >  2:1090560/4‑150 (MQ=255)                                                                                                tCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCTTCGGCTGTTAAGCAGGt                                                   >  2:793003/1‑150 (MQ=255)                                                                                                tCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGt                                                   <  2:815826/150‑1 (MQ=255)                                                                                                 cTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTc                                                  >  2:82842/1‑150 (MQ=255)                                                                                                   aCCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCAc                                                >  2:537076/1‑150 (MQ=255)                                                                                                    ccTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACt                                               <  1:96624/150‑1 (MQ=255)                                                                                                    ccTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACt                                               >  1:986301/1‑150 (MQ=255)                                                                                                      ttGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGa                                             >  1:1090990/1‑150 (MQ=255)                                                                                                       tGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGaa                                            <  1:139504/150‑1 (MQ=255)                                                                                                       tGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGaa                                            <  2:1090990/150‑1 (MQ=255)                                                                                                        gTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGACGTGTTCGCCG‑GCAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAg                                           <  1:1027008/150‑1 (MQ=255)                                                                                                              ttGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAAt                                     <  1:1206978/150‑1 (MQ=255)                                                                                                                      gCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGACGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTgcc                             >  1:110703/1‑150 (MQ=255)                                                                                                                          cGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGCACAGGTCACTGAAGCTCAATGCCTTGCCGCCa                         <  1:741672/150‑1 (MQ=255)                                                                                                                                gTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTgcc                             <  1:492243/140‑1 (MQ=255)                                                                                                                                gTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTgcc                             >  2:492243/1‑140 (MQ=255)                                                                                                                                 ttCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGtt                  >  1:947976/1‑150 (MQ=255)                                                                                                                                 ttCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACCTCTTTCCTGTGCCGGCTAACTGGTGAGGGAAGCTCAATGCCTTGCCGCCAGGCAGtt                  <  1:603991/150‑1 (MQ=255)                                                                                                                                   cAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGc                <  1:476299/150‑1 (MQ=255)                                                                                                                                     gCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAg              >  2:438920/1‑150 (MQ=255)                                                                                                                                      cGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAgg             >  1:504652/1‑150 (MQ=255)                                                                                                                                        aaGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCg           >  1:136105/1‑150 (MQ=255)                                                                                                                                            cTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTc       >  1:744995/1‑150 (MQ=255)                                                                                                                                             tt‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa      >  2:901623/1‑150 (MQ=255)                                                                                                                                              tggCCGCTGTGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGAGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa      >  2:392049/3‑150 (MQ=255)                                                                                                                                                 ccGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAgg                       >  2:1154629/1‑130 (MQ=255)                                                                                                                                                 ccGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGCTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAgg                       <  1:1154629/130‑1 (MQ=255)                                                                                                                                                  cGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCt  >  1:893811/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GCTTTTGCAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTT‑GCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCT  >  NC_002947/77351‑77644

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