Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	58,439	G→C	8.7%	intergenic (‑350/+338)	PP_5424 ← / ← PP_0050	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	58,439	0	G	C	8.7%	211.8 / 8.8	69	intergenic (‑350/+338)	PP_5424/PP_0050	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (3/3);  total (36/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 = 2.04e-01

CCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAA  >  NC_002947/58293‑58584                                                                                                                                                   |                                                                                                                                                    ccGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGtt                                                                                                                                                   >  3:79636/1‑149 (MQ=255)      aGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCctct                                                                                                                                              >  3:148012/1‑149 (MQ=255)         ccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                           <  5:126090/149‑3 (MQ=255)         ccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                           <  5:283076/149‑3 (MQ=255)          cggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAag                                                                                                                                          <  5:196435/149‑1 (MQ=255)                   tCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcgg                                                                                                                                 >  5:39443/1‑149 (MQ=255)                      tGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCGTTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcggcgg                                                                                                                              >  8:334669/1‑149 (MQ=255)                         aTCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCa                                                                                                                           >  3:276945/1‑149 (MQ=255)                            aGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGACCCGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCAtt                                                                                                                         <  2:294279/148‑1 (MQ=255)                             gATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGc                                                                                                                       <  7:29608/149‑1 (MQ=255)                                   gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTc                                                                                                                  <  6:11544/148‑1 (MQ=255)                                        acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGAtt                                                                                                             <  7:269066/148‑1 (MQ=255)                                          acGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCt                                                                                                           >  7:90343/1‑148 (MQ=255)                                            gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                           >  4:119314/1‑114 (MQ=255)                                            gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                           <  3:119314/114‑1 (MQ=255)                                            gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGATGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                        >  7:72512/1‑149 (MQ=255)                                            gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                        <  2:22693/149‑1 (MQ=255)                                            gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACAGCGCATTTGTCGGGCCTGGGTACCAGCAGGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCAGTGCCAGTCGGATtgtgca                                                                                                        >  6:301687/1‑147 (MQ=255)                                             cgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAAc                                                                                                       <  8:236230/149‑1 (MQ=255)                                                 ccGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCaa                                                                                                   <  4:276945/149‑1 (MQ=255)                                                           tAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGATCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGa                                                                                        <  6:14044/149‑1 (MQ=255)                                                              ccTTGCAATTGAGGGCTGCAACTGCGCGTTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTa                                                                                     <  7:334669/149‑1 (MQ=255)                                                              ccTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCAtt                                                                                                                         <  6:233362/114‑1 (MQ=255)                                                               cTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCt                                                                                                           >  3:240263/1‑127 (MQ=255)                                                               cTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCt                                                                                                           <  4:240263/127‑1 (MQ=255)                                                                ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑g                                                                                   >  8:189081/1‑148 (MQ=255)                                                                ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑g                                                                                   >  8:6631/1‑148 (MQ=255)                                                                  gCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAg                                                                                >  3:222602/1‑149 (MQ=255)                                                                   cAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAg                                                                                <  4:274332/148‑1 (MQ=255)                                                                    aaTTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGc                                                                               >  6:224025/1‑148 (MQ=255)                                                                                cAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCt                                                                   <  2:205352/148‑1 (MQ=255)                                                                                 aaCTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGATGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGc                                                                 <  8:72512/149‑1 (MQ=255)                                                                                 aaCTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGc                                                                 <  7:189081/149‑1 (MQ=255)                                                                                 aaCTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGc                                                                 <  4:222602/149‑1 (MQ=255)                                                                                   cTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGa                                                               >  1:350693/1‑149 (MQ=255)                                                                                                 ggCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATccc                                                 >  8:186639/1‑149 (MQ=255)                                                                                                  gCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCGGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCt                                                >  5:345912/1‑149 (MQ=255)                                                                                                  gCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCt                                                >  2:292623/1‑149 (MQ=255)                                                                                                       ggTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTAcc                                            <  7:186639/148‑1 (MQ=255)                                                                                                        gTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCg                                           >  5:199745/1‑148 (MQ=255)                                                                                                          aCCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAAc                                                                                                       >  2:33305/1‑88 (MQ=255)                                                                                                          aCCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAAc                                                                                                       <  1:33305/88‑1 (MQ=255)                                                                                                           ccAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGaa                                         >  5:116197/1‑147 (MQ=255)                                                                                                             aGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGt                                     >  2:315976/1‑149 (MQ=255)                                                                                                              gCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGaaaa                                       <  4:325091/146‑1 (MQ=255)                                                                                                              gCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGaaaa                                       >  3:325091/1‑146 (MQ=255)                                                                                                              gCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGt                                     >  6:330164/1‑148 (MQ=255)                                                                                                              gCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTc                                    <  8:90343/149‑1 (MQ=255)                                                                                                                       cccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                                  <  1:207411/142‑1 (MQ=255)                                                                                                                       cccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAAACCAAGC‑TCGCCAGACTA‑GGACCTAAAGCATCGCAGCGACTACTCCTCAACCCTTACCGAAAAGTCAc                                  >  2:207411/1‑142 (MQ=255)                                                                                                                        ccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcggt                          >  6:270382/1‑149 (MQ=255)                                                                                                                         cGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcggtc                         >  2:41978/1‑149 (MQ=255)                                                                                                                         cGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcggtc                         <  7:6631/149‑1 (MQ=255)                                                                                                                         cGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGACAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcggtc                         >  3:76274/1‑149 (MQ=255)                                                                                                                          gCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcggtcg                        <  5:270382/149‑1 (MQ=255)                                                                                                                                 tGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACg                  <  5:330164/148‑1 (MQ=255)                                                                                                                                   gcagcaGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGc                <  6:199745/148‑1 (MQ=255)                                                                                                                                     ttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCt                                                                              <  1:14132/82‑1 (MQ=39)                                                                                                                                     ttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCt                                                                              >  2:14132/3‑84 (MQ=39)                                                                                                                                     ttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCt             >  7:41257/3‑149 (MQ=255)                                                                                                                                     ttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCt             >  1:88737/3‑149 (MQ=255)                                                                                                                                        agagGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCttt           >  5:195013/1‑148 (MQ=255)                                                                                                                                        agagGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc          >  5:368410/1‑149 (MQ=255)                                                                                                                                          agGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc          >  5:218060/1‑147 (MQ=255)                                                                                                                                          agGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc          <  6:218060/147‑1 (MQ=255)                                                                                                                                          agGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCACCAGACTA‑GGAGCTAAACGATCGCTGCGTCTACTCCTCAt                                                    <  5:241080/105‑1 (MQ=255)                                                                                                                                          agGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCACCAGACTA‑GGAGCTAAACGATCGCTGCGTCTACTCCTCAt                                                    >  6:241080/1‑105 (MQ=255)                                                                                                                                              aGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGCTCGCCAGACTATGTCGAATAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAgg      <  1:292623/149‑1 (MQ=255)                                                                                                                                                cTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCaa  >  7:188652/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGG‑TCGCCAGACTA‑GGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAA  >  NC_002947/58293‑58584

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