Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,536,585 C→G 8.9% R75G (CGT→GGT)  PP_3125 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,536,5850CG 8.9% 185.5 / 5.5 62R75G (CGT→GGT) PP_3125Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base G (2/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

ACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACAC  >  NC_002947/3536438‑3536725
                                                                                                                                                   |                                                                                                                                            
aCCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGAGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCg                                                                                                                                             >  5:97320/1‑149 (MQ=255)
 ccAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCg                                                                                                                                             >  4:79211/1‑148 (MQ=255)
         tACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAgcg                                                                                                                                    >  2:239549/1‑149 (MQ=255)
            cTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGCGCTCGGTACGGCGCACCa                                                                                                                                 <  6:171189/149‑1 (MQ=255)
                 cAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGccc                                                                                                                            >  7:72305/1‑149 (MQ=255)
                      tCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGtt                                                                                                                       >  1:130710/1‑149 (MQ=255)
                           cgcgGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCTGTACCGAGCGCCAGGCCCTGGTTGGCGa                                                                                                                  >  6:44058/1‑149 (MQ=255)
                              ggCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATgg                                                                                                               <  7:298266/149‑1 (MQ=255)
                                     acacGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTgg                                                                                                                         <  1:318174/132‑1 (MQ=255)
                                     acacGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTgg                                                                                                                         >  2:318174/1‑132 (MQ=255)
                                     acacGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAggg                                                                                                        >  4:57811/1‑149 (MQ=255)
                                         gCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCgg                                                                                                    <  1:38608/149‑1 (MQ=255)
                                          cccTCGAACGCTGCGGGGCTGAAACTGGCTAGGGCGCTGAACGTCAATGTCGAATAGCTCTTCGCCGACTAGATGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCggg                                                                                                   <  1:239549/149‑1 (MQ=255)
                                           ccTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGc                                                                                                       <  6:215316/144‑1 (MQ=255)
                                           ccTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGc                                                                                                       >  5:215316/1‑144 (MQ=255)
                                                aTCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAgg                                                                                                         >  5:342543/1‑137 (MQ=255)
                                                aTCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAgg                                                                                                         <  6:342543/137‑1 (MQ=255)
                                                   gCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGc                                                                                          >  8:86031/1‑149 (MQ=255)
                                                    cTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGcc                                                                                         <  8:72305/149‑1 (MQ=255)
                                                       cggcgCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTc                                                                                      <  4:289097/149‑1 (MQ=255)
                                                             tGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGt                                                                                >  4:305930/1‑149 (MQ=255)
                                                               aaaCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCa                                                                              <  2:130710/149‑1 (MQ=255)
                                                               aaaCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCa                                                                              <  2:10067/149‑1 (MQ=255)
                                                                 aCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGt                                                                                >  1:149492/1‑145 (MQ=255)
                                                                 aCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCaca                                                                            >  7:243652/1‑147 (MQ=255)
                                                                 aCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCATGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTgg                                                                                                                     >  5:265737/1‑108 (MQ=255)
                                                                 aCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCATGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTgg                                                                                                                     <  6:265737/108‑1 (MQ=255)
                                                                       gAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCg                                                                                                                   <  1:238726/104‑1 (MQ=255)
                                                                       gAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCg                                                                                                                   >  2:238726/1‑104 (MQ=255)
                                                                          ggCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGt                                                                                <  1:9282/136‑1 (MQ=255)
                                                                          ggCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGt                                                                                >  2:9282/1‑136 (MQ=255)
                                                                             gcTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGAtt                                                                          <  6:32880/139‑1 (MQ=255)
                                                                             gcTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGAtt                                                                          >  5:32880/1‑139 (MQ=255)
                                                                                 aaCGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGc                                                              <  3:79211/147‑1 (MQ=255)
                                                                                 aaCGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCct                                                            >  4:131149/1‑149 (MQ=255)
                                                                                       aaTGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCtgctgc                                                       <  3:305930/148‑1 (MQ=255)
                                                                                        aTGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCg                                                     <  3:57811/149‑1 (MQ=255)
                                                                                             gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTc                                                 <  3:74338/148‑1 (MQ=255)
                                                                                             gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCt                                                >  1:136072/1‑149 (MQ=255)
                                                                                             gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCt                                                >  5:227735/1‑149 (MQ=255)
                                                                                             gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGCGCTCGGTACGGCGCACCAg                                                                                                                                >  4:9668/1‑69 (MQ=25)
                                                                                             gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGCGCTCGGTACGGCGCACCAg                                                                                                                                <  3:9668/69‑1 (MQ=25)
                                                                                               agagCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCGGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGa                                              <  5:44058/149‑1 (MQ=255)
                                                                                                 agCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAg                                                                                                          <  8:87308/87‑1 (MQ=255)
                                                                                                 agCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAg                                                                                                          >  7:87308/1‑87 (MQ=255)
                                                                                                 agCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGAtt                                            <  3:131149/149‑1 (MQ=255)
                                                                                                   ctctTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGATCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCa                                          <  6:97320/149‑1 (MQ=255)
                                                                                                        cGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCAcc                                                                                   <  2:274655/103‑1 (MQ=255)
                                                                                                        cGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCAcc                                                                                   >  1:274655/1‑103 (MQ=255)
                                                                                                                 gCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAg                             <  2:149492/148‑1 (MQ=255)
                                                                                                                    ggCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGtt                           >  6:103238/1‑147 (MQ=255)
                                                                                                                     gCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCa                                                                <  6:254735/109‑1 (MQ=255)
                                                                                                                     gCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGtt                           <  7:136664/146‑1 (MQ=255)
                                                                                                                     gCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGtt                           >  8:136664/1‑146 (MQ=255)
                                                                                                                               gCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCCGCCTGCTGCCGTTTGTGATTCAGCCGCCGTCCGAGTTCCGCGACCCTACg              >  3:338574/1‑149 (MQ=255)
                                                                                                                                  gCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGtt            >  4:100406/1‑148 (MQ=255)
                                                                                                                                   ccagggCCTGGCGCTCGGTACGGCGCACCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTc           >  2:108149/6‑148 (MQ=255)
                                                                                                                                   cTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGAcc                   >  6:288703/1‑140 (MQ=255)
                                                                                                                                   cTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGAcc                   <  5:288703/140‑1 (MQ=255)
                                                                                                                                    tACAGCCTGGGGCGCCGTCCCGTGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGGCTGCTGGCGTTTCTGATTCAGCCGCCGTCCGTGTTCAGCGACCCTACGTTCaa         <  2:284746/149‑1 (MQ=255)
                                                                                                                                        gCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGaa     <  8:243652/149‑1 (MQ=255)
                                                                                                                                             gTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAacac  >  6:171278/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
ACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACAC  >  NC_002947/3536438‑3536725

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