Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,741,951	A→G	7.6%	F87L (TTC→CTC)	PP_5039 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,741,951	0	A	G	7.6%	220.6 / 6.6	79	F87L (TTC→CTC)	PP_5039	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (3/3);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

CGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGCGTCG  >  NC_002947/5741803‑5742099                                                                                                                                                     |                                                                                                                                                      cGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                      >  1:710333/1‑150 (MQ=255)    tCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcg                                                                                                                                                   <  2:1063186/150‑1 (MQ=255)      cACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcgcg                                                                                                                                                 >  1:193841/1‑150 (MQ=255)      cACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcgcg                                                                                                                                                 >  1:1166953/1‑150 (MQ=255)       aCGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcgcgc                                                                                                                                                <  1:1606690/150‑1 (MQ=255)         gCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCgcgg                                                                                                                                              >  2:1666962/1‑150 (MQ=255)          cacaGGTTCTGTCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGgcgaccggc                                                                                                                                             <  2:710307/150‑7 (MQ=255)          cacaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGTATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGgcgaccggc                                                                                                                                             <  2:163840/150‑7 (MQ=255)           acaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGgcgaccggcg                                                                                                                                            >  1:836278/1‑143 (MQ=255)               gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                      >  1:892605/1‑136 (MQ=255)               gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                      <  2:892605/136‑1 (MQ=255)               gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTACTCGTATGCCGTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTc                                                                                                                                        >  1:1772373/1‑150 (MQ=255)                  cTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTcg‑cg                                                                                                                                    >  2:1519927/1‑150 (MQ=255)                      cAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAg                                                                                                                                >  1:1200002/1‑150 (MQ=255)                         tCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATc                                                                                                                             >  1:542535/1‑150 (MQ=255)                         cctaGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATc                                                                                                                             >  2:505077/4‑150 (MQ=255)                           cAGCACTACGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCCTCTGCAGGCCCTGCATGGCGACCACCTGTGCACGCCGGTCCCCCTCGAAGGCGCGCGGCGTGTCG‑CGGCAGATccc                                                                                                                           >  2:90827/1‑148 (MQ=255)                            atcattgcGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATg                                                                                                                          <  2:367075/143‑1 (MQ=255)                             gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGc                                                                                                                         <  2:193841/150‑1 (MQ=255)                             gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGc                                                                                                                         <  2:542535/150‑1 (MQ=255)                                 ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcg                                                                                                                                                   >  2:1615552/1‑121 (MQ=255)                                 ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcg                                                                                                                                                   <  1:1615552/121‑1 (MQ=255)                                   cGGCATGCTCCAGTACCAGCAGCACGCCCACGCGCTGGCACGGCAGCAGCTCGTAGTCCCAGCCCTTCTGCCGCTGCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCa                                                                                                                   <  1:1067265/150‑1 (MQ=255)                                     gCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGc                                                                                                                 >  1:1811049/1‑150 (MQ=255)                                       attcTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTg                                                                                                               >  1:274651/4‑150 (MQ=255)                                       aTGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTg                                                                                                               <  2:1488370/150‑1 (MQ=255)                                        tGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTgg                                                                                                              <  2:1440719/150‑1 (MQ=255)                                          cTCCAGTACCAGCAGCACGAACACGCACAGGATCGTCAGCAGCTAATAGTCCCAGCCCTGCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGCAGCCGCCCGGCGTGTCG‑CGCTAGATCATGCCTGGCAGCTGGTc                                                                                                            <  2:1552445/150‑1 (MQ=255)                                            ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                          <  1:1844158/150‑1 (MQ=255)                                            ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                          <  1:1844185/150‑1 (MQ=255)                                            ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATCGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                          <  1:233370/150‑1 (MQ=255)                                                   cagcagCACGAACACGCGCTGGATCGGCAGCAGCTAGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGGAAGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCaaaaa                                                                                                   >  2:61714/1‑149 (MQ=255)                                                      cagcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcag                                                                                                <  1:395913/150‑1 (MQ=255)                                                          aCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAgg                                                                                            <  1:1197836/150‑1 (MQ=255)                                                                acGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCcaagc                                                                                      <  1:902291/150‑1 (MQ=255)                                                                        gATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTg                                                                              >  1:166024/1‑150 (MQ=255)                                                                            ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                          <  1:248282/150‑1 (MQ=255)                                                                            ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                          <  1:1519927/150‑1 (MQ=255)                                                                             gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                                         >  2:522324/1‑150 (MQ=255)                                                                             gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                                         >  2:1280966/1‑150 (MQ=255)                                                                                gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                                      <  1:68905/150‑1 (MQ=255)                                                                                gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                                      <  2:1166953/150‑1 (MQ=255)                                                                                gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                                      <  2:274651/150‑1 (MQ=255)                                                                                   gcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAGGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                                                   >  1:1100788/1‑150 (MQ=255)                                                                                   gcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgat                                                                   <  2:935835/150‑3 (MQ=255)                                                                                    cTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgctg                                                                  >  2:1001662/1‑150 (MQ=255)                                                                                           tCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGGCGACCGGCgtg                                                                                                                                          >  2:1519865/1‑72 (MQ=18)                                                                                           tCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGGCGACCGGCgtg                                                                                                                                          <  1:1519865/72‑1 (MQ=18)                                                                                             ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                                         >  1:953992/1‑150 (MQ=255)                                                                                             ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                                         >  1:672591/1‑150 (MQ=255)                                                                                               aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                          <  2:1510477/99‑1 (MQ=255)                                                                                               aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                          >  1:1510477/1‑99 (MQ=255)                                                                                                 cccTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGc                                                     <  2:672591/150‑1 (MQ=255)                                                                                                  ccTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGcc                                                    >  1:1083037/1‑150 (MQ=255)                                                                                                    ttCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCa                                                  <  2:1772373/150‑1 (MQ=255)                                                                                                       tGCAGGCCCTGCATGGCGACCTCCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCAccc                                               <  1:1513856/150‑1 (MQ=255)                                                                                                        gCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCg                                              >  1:1694029/1‑150 (MQ=255)                                                                                                         cAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGc                                             >  2:951147/1‑150 (MQ=255)                                                                                                            gCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAg                                          >  2:1737263/1‑150 (MQ=255)                                                                                                             cccTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGAG‑TAGTCCAGGCGTCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAgg                                         <  2:50324/150‑1 (MQ=255)                                                                                                              ccTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAggg                                        >  2:138713/1‑150 (MQ=255)                                                                                                               cTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATca                                                                                                 <  2:1633754/92‑1 (MQ=255)                                                                                                               cTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATca                                                                                                 >  1:1633754/1‑92 (MQ=255)                                                                                                               cTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGc                                       <  1:138713/150‑1 (MQ=255)                                                                                                                  cATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGccc                                                                                <  1:321309/106‑1 (MQ=255)                                                                                                                  cATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGccc                                                                                >  2:321309/1‑106 (MQ=255)                                                                                                                   aTGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGc                                   <  1:1001662/150‑1 (MQ=255)                                                                                                                   aTGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGc                                   <  1:929307/150‑1 (MQ=255)                                                                                                                           caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAg                           <  1:1737263/150‑1 (MQ=255)                                                                                                                            accTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAgg                          >  1:123564/1‑150 (MQ=255)                                                                                                                            accTGGGCACGCCGGTCGCCCTCGAAGGCGCGAGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAgg                          >  1:439834/1‑150 (MQ=255)                                                                                                                              cTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTc                        <  2:300604/150‑1 (MQ=255)                                                                                                                                      gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                                                   >  1:1443282/1‑99 (MQ=255)                                                                                                                                      gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                                                   <  2:1443282/99‑1 (MQ=255)                                                                                                                                      gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACa                >  1:240759/1‑150 (MQ=255)                                                                                                                                               ccTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATggg       >  2:533197/1‑150 (MQ=255)                                                                                                                                               ccTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATggg       >  2:342536/1‑150 (MQ=255)                                                                                                                                                 tCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGCg     >  1:1815192/1‑150 (MQ=255)                                                                                                                                                 tCGAAGGCGAGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGTCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGCg     >  1:1735860/1‑150 (MQ=255)                                                                                                                                                    aaGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTAGCCAAACAGCGCATGGGCGTCg  <  1:522324/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGCGTCG  >  NC_002947/5741803‑5742099

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