Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,886,343 C→G 6.8% G505A (GGC→GCC)  PP_3430 ← two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,886,3430CG 6.8% 185.0 / 4.4 59G505A (GGC→GCC) PP_3430two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base G (2/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

CGGCGTCGACTCGATCCAGGCATACCCCCCCGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAA  >  NC_002947/3886201‑3886491
                                                                                                                                               |                                                                                                                                                    
cGGCGTCGACTCGATCCAGGCATACCCCCCCGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATgg                                                                                                                                                 <  3:34558/148‑1 (MQ=255)
  gCGTCGACTCGATCCAGGCATACCCCCCCGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGccatgccatg                                                                                                                                                  >  2:7291/1‑145 (MQ=255)
  gCGTCGACTCGATCCAGGCATACCCCCCCGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGccatgccatg                                                                                                                                                  <  1:7291/145‑1 (MQ=255)
    gTCGACTCGATCCAGGCATACCCCCCCGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCcatc                                                                                                                                            <  5:78284/149‑1 (MQ=255)
       gACTCGATCCAGGCATACCCCCCCGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGGCATGGccgcggc                                                                                                                                          <  8:79481/148‑6 (MQ=255)
                cAGGCATACCCCCCCGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATg                                                                                                                                <  6:272522/149‑1 (MQ=255)
                        cccccccGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGcc                                                                                                                         <  2:310427/148‑1 (MQ=255)
                           ccccGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCcaca                                                                                                                      >  4:222013/1‑148 (MQ=255)
                               gACTGTCGGACGAAGCCATACACCCTCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTc                                                                                                                 <  6:227985/149‑1 (MQ=255)
                                    tCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGccc                                                                                                             <  3:294849/148‑1 (MQ=255)
                                       gACGAAGCCATACACCATCGTCAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGc                                                                                                                               >  7:100239/1‑127 (MQ=255)
                                       gACGAAGCCATACACCATCGTCAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGc                                                                                                                               <  8:100239/127‑1 (MQ=255)
                                       gACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATTCAAGGCCACATAGTCGCCccgcg                                                                                                         >  8:268603/1‑149 (MQ=255)
                                            aGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAg                                                                                                    <  6:73444/149‑1 (MQ=255)
                                              ccATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAACTCGGTTCGAAAGCCCTATTGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGc                                                                                                   <  1:265452/149‑1 (MQ=255)
                                                    acCATCGACTGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGccatgccatg                                                                                                                                                  <  4:193638/95‑1 (MQ=255)
                                                    acCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGccatgccatg                                                                                                                                                  >  3:193638/1‑95 (MQ=255)
                                                    acCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCtt                                                                                             >  8:55761/1‑148 (MQ=255)
                                                    acCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTc                                                                                            <  3:222013/149‑1 (MQ=255)
                                                      cATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCAt                                                                                          >  2:266752/1‑149 (MQ=255)
                                                       aTCGACAGGCCAGGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATc                                                                                         <  7:32666/149‑1 (MQ=255)
                                                       aTCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATc                                                                                         <  5:120177/149‑1 (MQ=255)
                                                       aTCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATc                                                                                         >  8:22020/1‑149 (MQ=255)
                                                       aTCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATc                                                                                         >  1:95578/1‑149 (MQ=255)
                                                               gCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGaaa                                                                                  >  1:159530/1‑148 (MQ=255)
                                                                 cAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAAc                                                                               >  5:123150/1‑149 (MQ=255)
                                                                    gCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCTCCCGGCTGCAAGCCACCTTCATCCGGGAAAAAcgcg                                                                            <  7:268603/149‑1 (MQ=255)
                                                                    gCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAAcgcg                                                                            >  7:286381/1‑149 (MQ=255)
                                                                            ccACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCcg                                                                    >  1:191233/1‑149 (MQ=255)
                                                                             cACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATACTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCcgt                                                                   >  7:224984/1‑149 (MQ=255)
                                                                                gCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCcgttcg                                                                >  7:306446/1‑149 (MQ=255)
                                                                                  ccGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCcgttcgtt                                                              <  2:95578/149‑1 (MQ=255)
                                                                                        ggTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACgct                                                        >  4:267132/1‑149 (MQ=255)
                                                                                          tttGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACgctgc                                                      >  3:41993/1‑149 (MQ=255)
                                                                                              gTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGAt                                                  >  8:10150/1‑149 (MQ=255)
                                                                                               ttGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAAcgcg                                                                            >  1:216631/1‑122 (MQ=255)
                                                                                               ttGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAAcgcg                                                                            <  2:216631/122‑1 (MQ=255)
                                                                                               ttGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATg                                                 >  8:55262/1‑149 (MQ=255)
                                                                                               ttGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATg                                                 >  3:260453/1‑149 (MQ=255)
                                                                                                  tgtAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGt                                               <  8:306446/148‑1 (MQ=255)
                                                                                                   gtAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGa                                             >  7:292814/1‑149 (MQ=255)
                                                                                                    tAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGa                                             >  1:283140/1‑148 (MQ=255)
                                                                                                      aa‑aCGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGGTGCTGATGGTGACAt                                          >  7:90112/1‑149 (MQ=255)
                                                                                                       a‑aCGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACAtt                                         >  3:32923/1‑149 (MQ=255)
                                                                                                       a‑aCGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACAtt                                         <  4:41993/149‑1 (MQ=255)
                                                                                                           ggTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGcc                                      <  2:191233/149‑1 (MQ=255)
                                                                                                             ttCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCt                                    >  2:316020/1‑149 (MQ=255)
                                                                                                                 aaaGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTc                                >  4:225276/1‑149 (MQ=255)
                                                                                                                        tGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCAc                         >  6:30058/1‑149 (MQ=255)
                                                                                                                           cGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCt                      <  8:90112/149‑1 (MQ=255)
                                                                                                                             atatCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTc                    <  6:123150/149‑1 (MQ=255)
                                                                                                                              tatCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCa                   >  5:59265/1‑149 (MQ=255)
                                                                                                                               atCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCAc                  >  5:265978/1‑149 (MQ=255)
                                                                                                                                 ccgccgCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCAcgc                >  5:73124/1‑149 (MQ=255)
                                                                                                                                  cgccgcGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCAcgc                <  4:32923/148‑1 (MQ=255)
                                                                                                                                     gatggCCATGGCATGGCCGTCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCAt            >  4:132913/4‑149 (MQ=255)
                                                                                                                                     gatggCCATGGCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGa                                                            <  3:24984/98‑1 (MQ=255)
                                                                                                                                     gatggCCATGGCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGa                                                            >  4:24984/4‑101 (MQ=255)
                                                                                                                                      gcgGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCAtt           <  3:225276/149‑1 (MQ=255)
                                                                                                                                        ggCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTgg         <  6:73124/149‑1 (MQ=255)
                                                                                                                                        ggCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTgg         <  6:59265/149‑1 (MQ=255)
                                                                                                                                        ggCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTgg         <  8:286381/149‑1 (MQ=255)
                                                                                                                                        ggCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTgg         <  8:292814/149‑1 (MQ=255)
                                                                                                                                           catggcatggCCAACATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGca       >  7:196410/1‑148 (MQ=255)
                                                                                                                                               ccatgGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACaa  <  7:55761/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CGGCGTCGACTCGATCCAGGCATACCCCCCCGACTGTCGGACGAAGCCATACACCATCGACAGGCCAAGCCCGGCCCCACGCCCGAGAGGTTTGGTTGTGTAAA‑ACGGTTCGAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAA  >  NC_002947/3886201‑3886491

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