Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,356,518	C→A	7.6%	R748L (CGC→CTC)	infB ←	translation initiation factor IF‑2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,356,518	0	C	A	7.6%	208.7 / 6.3	66	R748L (CGC→CTC)	infB	translation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (3/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-01

GCCTTCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGAAGCCGAACAGTAC  >  NC_002947/5356370‑5356666                                                                                                                                                     |                                                                                                                                                     gCCTTCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTcgc                                                                                                                                                      <  6:80107/149‑1 (MQ=255) gCCTTCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTcgc                                                                                                                                                      <  2:226599/149‑1 (MQ=255)    ttCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGaa                                                                                                                                                   >  4:46042/1‑149 (MQ=255)    ttCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGaa                                                                                                                                                   >  8:5392/1‑149 (MQ=255)    ttCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGa                                                                                                                                                    >  3:202488/1‑148 (MQ=255)      cGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGaa                                                                                                                                                   <  5:102526/147‑1 (MQ=255)       gAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACAt                                                                                                                                                <  5:36905/149‑1 (MQ=255)         aGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCg                                                                                                                                              <  4:270413/149‑1 (MQ=255)            tCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTg                                                                                                                                           <  4:199476/149‑1 (MQ=255)               caACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGcc                                                                                                                                         >  6:165563/1‑148 (MQ=255)                aaCGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGa                                                                                                                                       <  4:202488/149‑1 (MQ=255)                aaCGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGa                                                                                                                                       <  3:46042/149‑1 (MQ=255)                aaCGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGa                                                                                                                                       <  3:310651/149‑1 (MQ=255)                aaCGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGa                                                                                                                                       <  3:254994/149‑1 (MQ=255)                 aCGTCATCACGCAGTACGCGAATCGGCCGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAAAGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAg                                                                                                                                      <  8:70030/149‑1 (MQ=255)                    tcatcaCGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCAt                                                                                                                                   <  3:251665/149‑1 (MQ=255)                      atcaCGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGc                                                                                                                                 >  5:320598/1‑149 (MQ=255)                         aCGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCgg                                                                                                                              >  1:214645/1‑149 (MQ=255)                         aCGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCgg                                                                                                                              >  1:209920/1‑149 (MQ=255)                                    gAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCcttc                                                                                                                   <  6:284792/149‑1 (MQ=255)                                      aTCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATCCCGGTCAGGGCcttctt                                                                                                                 <  8:116339/149‑1 (MQ=255)                                                   ggTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCgatg                                                                                                    <  6:320598/149‑1 (MQ=255)                                                   ggTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCgat                                                                                                     >  5:305434/1‑148 (MQ=255)                                                   ggTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCgat                                                                                                     >  1:266133/1‑148 (MQ=255)                                                    gTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCGCGAGAACATcctg                                                                                                                                            >  5:104564/1‑105 (MQ=255)                                                    gTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCGCGAGAACATcctg                                                                                                                                            <  6:104564/108‑4 (MQ=255)                                                    gTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCGCGAGAACATcctg                                                                                                                                            >  5:109395/1‑105 (MQ=255)                                                    gTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCGCGAGAACATcctg                                                                                                                                            <  6:109395/108‑4 (MQ=255)                                                          ggTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTACGGCGACACCCAGGATGTTCTCGCGAACATCCCTGCCGAGCATGCCGGTCACGGCCTTCTTGACGTCTTCGATGATGTCg                                                                                              >  8:283102/1‑148 (MQ=255)                                                            tACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCg                                                                                              >  3:22284/1‑146 (MQ=255)                                                                         cATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCgatgat                                                                                                  <  5:181677/129‑1 (MQ=255)                                                                         cATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCgatgat                                                                                                  >  6:181677/1‑129 (MQ=255)                                                                            aCAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAAc                                                                           >  3:120809/1‑149 (MQ=255)                                                                             cagccagcGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAAc                                                                           <  7:141369/148‑1 (MQ=255)                                                                                   gcGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATc                                                                    <  6:305434/149‑1 (MQ=255)                                                                                        ttcgcCGAAGTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGAc                                                               >  3:217848/3‑149 (MQ=255)                                                                                         gcgcCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGAc                                                               <  5:165563/148‑1 (MQ=255)                                                                                              gAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCt                                                         <  7:5392/149‑1 (MQ=255)                                                                                                aCTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGc                                                       <  2:266133/149‑1 (MQ=255)                                                                                                  ttCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAAc                                                                           >  1:16213/1‑127 (MQ=255)                                                                                                  ttCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAAc                                                                           <  2:16213/127‑1 (MQ=255)                                                                                                  ttCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTc                                                     >  7:99810/1‑149 (MQ=255)                                                                                                       cGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTcgacg                                                 >  8:132590/1‑148 (MQ=255)                                                                                                        gAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGAt                                               <  4:120809/149‑1 (MQ=255)                                                                                                           cGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCtt                                            <  7:146657/149‑1 (MQ=255)                                                                                                           cGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCtt                                            <  8:99810/149‑1 (MQ=255)                                                                                                           cGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCATGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCtt                                            >  7:235971/1‑149 (MQ=255)                                                                                                            ggAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTg                                           <  4:22284/149‑1 (MQ=255)                                                                                                             gAACACGTCACGTACTTCGGCGACACCCATGATGTTCTAGCGAACATCTCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTgc                                          <  2:214645/149‑1 (MQ=255)                                                                                                             gAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATACCGCATATCCAGACCTTCCTGCTCGACGATCTTgc                                          <  2:209920/149‑1 (MQ=255)                                                                                                             gAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTgc                                          >  4:298720/1‑149 (MQ=255)                                                                                                              aaCACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTgcg                                         <  1:171270/149‑1 (MQ=255)                                                                                                                      acgtACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCg                                  <  4:217848/148‑1 (MQ=255)                                                                                                                        gtACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCa                               >  4:182253/1‑149 (MQ=255)                                                                                                                        gtACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCa                               >  2:265674/1‑149 (MQ=255)                                                                                                                          aCTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATc                             >  8:258922/1‑149 (MQ=255)                                                                                                                           cTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCg                            >  5:212343/1‑149 (MQ=255)                                                                                                                            ttCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCgg                           >  8:66412/1‑149 (MQ=255)                                                                                                                                tcGACACCCATGATGTTCTCGCTAACATCGCTGCCGAGAATGCCGGTCAGGGCCTTCATGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGgcacg                       <  8:235971/148‑1 (MQ=255)                                                                                                                                     acCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGt                  >  6:5287/1‑149 (MQ=255)                                                                                                                                       ccAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGtt                 >  5:38964/1‑148 (MQ=255)                                                                                                                                       ccAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGtt                 >  1:277972/1‑148 (MQ=255)                                                                                                                                       ccAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGtt                 >  5:213361/1‑148 (MQ=255)                                                                                                                                       cagcgATGTTCGCGAGAACATC‑CTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTg                >  4:325490/5‑148 (MQ=255)                                                                                                                                        cAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTg                >  2:196539/1‑148 (MQ=255)                                                                                                                                        cAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGa               <  6:212343/149‑1 (MQ=255)                                                                                                                                               ttCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGAAGCCGaa        <  3:298720/149‑1 (MQ=255)                                                                                                                                                  tcGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGAAGCCGAACAg     >  1:296631/1‑149 (MQ=255)                                                                                                                                                     cgAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGAAGCCGAACAGTAc  >  2:1668/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCCTTCGAAGATCACAACGTCATCACGCAGTACGCGAATCGGACGGTTGCGGTACACGGTACCCTCGATGACCATACAGCCAGCGATGGCGCCGAACTTCGGCGAACGGAACACGTCACGTACTTCGGCGACACCCAGGATGTTCTCGCGAACATCGCTGCCGAGCATGCCGGTCAGGGCCTTCTTGACGTCTTCGATGATGTCGTAGATCACGTTGTAATAACGCATATCCAGACCTTCCTGCTCGACGATCTTGCGCGCACCGGCATCGGCACGCACGTTGAAGCCGAACAGTAC  >  NC_002947/5356370‑5356666

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