Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,518,539	A→T	5.4%	A14A (GCT→GCA)	PP_4850 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,518,539	0	A	T	5.4%	281.6 / 3.5	90	A14A (GCT→GCA)	PP_4850	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (43/41);  new base T (2/3);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTT  >  NC_002947/5518392‑5518687                                                                                                                                                      |                                                                                                                                                     gTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAg                                                                                                                                                     >  1:374579/1‑149 (MQ=255) gTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAg                                                                                                                                                     >  4:120324/1‑149 (MQ=255)    gACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCgg                                                                                                                                                  >  5:217419/1‑149 (MQ=255)     aCCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCgg                                                                                                                                                  <  4:303626/148‑1 (MQ=255)      ccAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCggg                                                                                                                                                 >  8:176080/1‑148 (MQ=255)       cAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGaa                                                                                                                                               >  6:145236/1‑149 (MQ=255)         gTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATg                                                                                                                                             <  6:290517/149‑1 (MQ=255)          ttCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CTGCggcattcc                                                                                                                                            >  1:6487/1‑143 (MQ=255)          ttCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGc                                                                                                                                            <  8:22582/149‑1 (MQ=255)             gAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGc                                                                                                                                         <  3:177619/149‑1 (MQ=255)              aaCACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCa                                                                                                                                        <  7:152616/149‑1 (MQ=255)               acacCAGCAGACCGGCCATCAGGTTTGGTAGATAGTTCAGGGGCAATGCCAGCAGGAGCCCGGCCATTACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACATCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCa                                                                                                                                        <  6:293725/148‑1 (MQ=255)               acacCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCg                                                                                                                                          >  7:151236/1‑146 (MQ=255)               aaacCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCACGGACAAACCCAGCAGAAGCACGGACATTACGATGATCAAGCTGGCCAGGGTCATCTGACGGGGGGTAAACGTCATACAGCCTCGA‑CGGCAGA‑CAGCGCGAATGCCGCa                                                                                                                                        >  5:293725/3‑148 (MQ=255)                 acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCa                                                                                                                                        <  8:151236/146‑1 (MQ=255)                  ccAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCt                                                                                                                                    <  5:195804/149‑1 (MQ=255)                   cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CTGCGGCATTCCCGCTGTCTg                                                                                                                                   <  6:17808/149‑1 (MQ=255)                   cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTg                                                                                                                                   >  1:123885/1‑149 (MQ=255)                    agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CTGCGGCATTCCCGCTGTCTGc                                                                                                                                  <  6:293968/149‑1 (MQ=255)                        gACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CTGCGGCATTCCCGCTGTCTGccgtc                                                                                                                              <  2:6487/149‑4 (MQ=255)                         aCCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGcc                                                                                                                             >  6:235847/1‑149 (MQ=255)                           cGGCCAGCAGGCTTGGCAGAAACTTCATGCGCAAAGCCAGCAGAAGCCCAGCAATGACGATGATCAGGCTGGCAAGTGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑AGCCAGA‑CAGCGGGAATGCAGAAGTCttcaaaca                                                                                                                             >  6:349878/1‑140 (MQ=255)                             gCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                         <  2:374579/149‑1 (MQ=255)                             gCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                         <  3:104885/149‑1 (MQ=255)                                    ggCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGTGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCa                                                                                                                  <  1:235738/149‑1 (MQ=255)                                     gCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                 <  6:235598/149‑1 (MQ=255)                                      cTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAgg                                                                                                                <  5:200257/149‑1 (MQ=255)                                       tGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGc                                                                                                               <  3:120324/149‑1 (MQ=255)                                           cagcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCCTACAGCCTCGA‑CGGCAGC‑CCGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAg                                                                                                           >  4:166445/1‑149 (MQ=255)                                           cagcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACCGGGGTAGCAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCa                                                                                                              >  1:341927/1‑148 (MQ=255)                                           cagcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAg                                                                                                           <  5:127781/149‑1 (MQ=255)                                           cagcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAg                                                                                                           <  5:235847/149‑1 (MQ=255)                                             gcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGgc                                                                                                         >  5:55423/1‑149 (MQ=255)                                                  ttCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGTCAGCCTTGGCGGCGCAGGCATAGGCGCAgg                                                                                                    >  5:208516/1‑149 (MQ=255)                                                     aGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCa                                                                                                 >  5:90211/1‑149 (MQ=255)                                                       gggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCa                                                                                                 <  5:127868/147‑1 (MQ=255)                                                       gggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCa                                                                                                 >  6:127868/1‑147 (MQ=255)                                                       gggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCAtt                                                                                               >  1:246347/1‑149 (MQ=255)                                                        ggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCAttt                                                                                              >  5:365616/1‑149 (MQ=255)                                                          cAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttct                                                                                            <  7:176080/149‑1 (MQ=255)                                                            aaGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttc                                                                                          <  2:123885/149‑1 (MQ=255)                                                             aGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGTGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCTCAGTCTGCCAGCCTTGGAGGCGCAGGCATAGGCGCAGGTCATttcttct                                                                                         <  5:349878/149‑1 (MQ=255)                                                                cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttctt                                                                                        >  7:319554/1‑147 (MQ=255)                                                                cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                      <  6:55423/149‑1 (MQ=255)                                                                 agcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGgcgc                                                                                                       <  3:189679/131‑1 (MQ=255)                                                                 agcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGgcgc                                                                                                       >  4:189679/1‑131 (MQ=255)                                                                      cagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCaa                                                                                 >  2:177942/1‑148 (MQ=255)                                                                                     gatgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCt                                                                                                                                    <  1:67453/82‑1 (MQ=255)                                                                                     gatgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCt                                                                                                                                    >  2:67453/1‑82 (MQ=255)                                                                                      atgatCAGGCTGGCCAGGGTGATCTGACGGGGTGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCtt                                                                <  1:177942/149‑1 (MQ=255)                                                                                      atgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCtt                                                                <  2:246347/149‑1 (MQ=255)                                                                                          tCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgcgc                                                            >  4:21735/1‑149 (MQ=255)                                                                                               cTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                             >  2:362177/1‑111 (MQ=255)                                                                                               cTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                             <  1:362177/111‑1 (MQ=255)                                                                                               cTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATc                                                       >  4:243503/1‑149 (MQ=255)                                                                                                      gggTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTc                                                <  3:21735/149‑1 (MQ=255)                                                                                                       ggTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGaa                                                                      <  6:321374/126‑1 (MQ=255)                                                                                                       ggTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGaa                                                                      >  5:321374/1‑126 (MQ=255)                                                                                                           aTCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGc                                           >  5:260784/1‑149 (MQ=255)                                                                                                               gACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCa                                                                          <  6:94349/114‑1 (MQ=255)                                                                                                               gACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCa                                                                          >  5:94349/1‑114 (MQ=255)                                                                                                                aCGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCAt                                                                         >  3:280972/1‑114 (MQ=255)                                                                                                                aCGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCAt                                                                         <  4:280972/114‑1 (MQ=255)                                                                                                                aCGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtg                                       <  3:243503/148‑1 (MQ=255)                                                                                                                       gTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGc                               >  4:148907/1‑149 (MQ=255)                                                                                                                       gTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGc                               <  6:217419/149‑1 (MQ=255)                                                                                                                       gTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGc                               >  7:52927/1‑149 (MQ=255)                                                                                                                       gTAAAGGTCATACAGCATCGA‑CGGAAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGTCGGCGCAGGAATAGGCGCAGGTAATTTATTCTTAAGGCAATCGCTCATGAAGGCCTTGCGCTCATCAACCTTCAGGGATTtggtggtgtc                               >  2:324581/1‑147 (MQ=255)                                                                                                                            ggTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAg                          <  7:147525/149‑1 (MQ=255)                                                                                                                             gTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGcc                                                                  <  3:144685/108‑1 (MQ=255)                                                                                                                             gTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGcc                                                                  >  4:144685/1‑108 (MQ=255)                                                                                                                             gTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAg                          >  8:51185/1‑148 (MQ=255)                                                                                                                                aTACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCAtt                      >  2:221648/1‑149 (MQ=255)                                                                                                                                 tACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                     >  3:330526/1‑149 (MQ=255)                                                                                                                                  aCAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGGGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                     >  2:334885/1‑148 (MQ=255)                                                                                                                                  aCAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                     >  6:271479/1‑148 (MQ=255)                                                                                                                                  aCAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                     >  2:332345/1‑148 (MQ=255)                                                                                                                                  aCAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                     >  2:376512/1‑148 (MQ=255)                                                                                                                                  aCAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                     >  4:340931/1‑148 (MQ=255)                                                                                                                                    aGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAg                  <  1:221648/149‑1 (MQ=255)                                                                                                                                        tCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCt              >  2:40061/1‑149 (MQ=255)                                                                                                                                         cGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCtt             >  6:34887/1‑149 (MQ=255)                                                                                                                                          gA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCtt             >  2:261152/1‑148 (MQ=255)                                                                                                                                          gA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCtt             <  6:260784/148‑1 (MQ=255)                                                                                                                                          gA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTc            <  8:319554/149‑1 (MQ=255)                                                                                                                                           a‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCa           <  6:90211/149‑1 (MQ=255)                                                                                                                                               gCAGA‑CTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAttt        >  4:20031/1‑149 (MQ=255)                                                                                                                                                  gA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAtttttt     <  1:261152/149‑1 (MQ=255)                                                                                                                                                   a‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTc    <  1:40061/149‑1 (MQ=255)                                                                                                                                                   a‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTc    <  4:330526/149‑1 (MQ=255)                                                                                                                                                      aGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTATTCATTTTTTCtt  >  2:380215/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGA‑CGGCAGA‑CAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTT  >  NC_002947/5518392‑5518687

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