Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,536,595 G→A 6.5% R78H (CGC→CAC)  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,5950GA 6.5% 183.4 / 2.5 61R78H (CGC→CAC) PP_3125Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (2/2);  total (30/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.63e-01

TACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCT  >  NC_002947/3536447‑3536742
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                          <  6:342543/137‑1 (MQ=255)
                                       aTCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAgg                                                                                                                          >  5:342543/1‑137 (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                                                                                                                                      <  6:265737/108‑1 (MQ=255)
                                                        aCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCATGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTgg                                                                                                                                      >  5:265737/1‑108 (MQ=255)
                                                              gAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCg                                                                                                                                    >  2:238726/1‑104 (MQ=255)
                                                              gAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCg                                                                                                                                    <  1:238726/104‑1 (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                                                                 >  1:136072/1‑149 (MQ=255)
                                                                                    gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCt                                                                 >  5:227735/1‑149 (MQ=255)
                                                                                    gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGCGCTCGGTACGGCGCACCAg                                                                                                                                                 <  3:9668/69‑1 (MQ=25)
                                                                                    gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGCGCTCGGTACGGCGCACCAg                                                                                                                                                 >  4:9668/1‑69 (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)
                                                                                                                                                   cgcCAGGCCCTGGTTGGCGATTGGGAAGGGCCGTGGTACGCGTCCCTCACCAGTCAGATTGGACAGCGGAGCCTGCTGCGGTTTGTGATTCAGCCGACGTCCGAGTTGAGCGACCCTACGTTCAAGGAACACCAGGGCGAATATTTCCt  <  4:338574/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCT  >  NC_002947/3536447‑3536742

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