Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,741,941	T→G	6.0%	D90A (GAC→GCC) ‡	PP_5039 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,741,941	0	T	G	6.0%	221.5 / 5.9	84	D90A (GAC→GCC) ‡	PP_5039	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/40);  new base G (2/3);  total (41/43)
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.51e-01

CATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACA  >  NC_002947/5741793‑5742085                                                                                                                                                     |                                                                                                                                                  caTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTc                                                                                                                                                  <  1:845279/150‑1 (MQ=255) caTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTc                                                                                                                                                  <  1:155592/150‑1 (MQ=255) caTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTCCATGGCGCCAACCTGAGCACGACCGTc                                                                                                                                                  <  2:504468/150‑1 (MQ=255)  aTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCg                                                                                                                                                 >  2:902291/1‑150 (MQ=255)   ttCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTATACAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGc                                                                                                                                                >  1:1452901/1‑150 (MQ=255)      gATAACGCTCCACGCACAGGTTCTGCGGCGCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCt                                                                                                                                             >  1:50324/1‑150 (MQ=255)      gATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCCCGCCGGTCGCCCt                                                                                                                                             >  1:1908941/1‑150 (MQ=255)           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                                                                                                                                        <  2:892605/136‑1 (MQ=255)                         gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                        >  1:892605/1‑136 (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)                                 aGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGc                                                                                                                                                <  2:1478169/120‑1 (MQ=255)                                 aGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGc                                                                                                                                                >  1:1478169/1‑120 (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:1844185/150‑1 (MQ=255)                                                      ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                            <  1:1844158/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:1519927/150‑1 (MQ=255)                                                                                      ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                            <  1:248282/150‑1 (MQ=255)                                                                                       gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                           >  2:522324/1‑150 (MQ=255)                                                                                       gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                           >  2:1280966/1‑150 (MQ=255)                                                                                          gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                        <  2:1166953/150‑1 (MQ=255)                                                                                          gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                        <  2:274651/150‑1 (MQ=255)                                                                                          gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                        <  1:68905/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                                                     <  2:1443282/99‑1 (MQ=255)                                                                                                                                                gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                                     >  1:1443282/1‑99 (MQ=255)                                                                                                                                                gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACa  >  1:240759/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACA  >  NC_002947/5741793‑5742085

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