Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,536,592 A→C 6.5% E77A (GAG→GCG)  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,5920AC 6.5% 186.3 / 4.7 62E77A (GAG→GCG) PP_3125Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base C (2/2);  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 = 4.75e-01

TACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGT  >  NC_002947/3536447‑3536738
                                                                                                                                                 |                                                                                                                                                  
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                                                                                                                    >  5:215316/1‑144 (MQ=255)
                                  ccTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGc                                                                                                                    <  6:215316/144‑1 (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:10067/149‑1 (MQ=255)
                                                      aaaCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCa                                                                                           <  2:130710/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                                                                                             >  2:9282/1‑136 (MQ=255)
                                                                 ggCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGt                                                                                             <  1:9282/136‑1 (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                                                             >  5:227735/1‑149 (MQ=255)
                                                                                    gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCt                                                             >  1:136072/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                                                                                                >  1:274655/1‑103 (MQ=255)
                                                                                               cGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCAcc                                                                                                <  2:274655/103‑1 (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                                <  5:288703/140‑1 (MQ=255)
                                                                                                                          cTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGAcc                                >  6:288703/1‑140 (MQ=255)
                                                                                                                           tACAGCCTGGGGCGCCGTCCCGTGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGGCTGCTGGCGTTTCTGATTCAGCCGCCGTCCGTGTTCAGCGACCCTACGTTCaa                      <  2:284746/149‑1 (MQ=255)
                                                                                                                               gCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGaa                  <  8:243652/149‑1 (MQ=255)
                                                                                                                                    gTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAacac               >  6:171278/1‑147 (MQ=255)
                                                                                                                                             aCCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGaa      >  5:53921/1‑147 (MQ=255)
                                                                                                                                               cGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGt  >  8:109021/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGT  >  NC_002947/3536447‑3536738

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