Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,797,370 C→G 8.8% I141M (ATC→ATG PP_2449 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,797,3700CG 8.8% 166.4 / 7.3 57I141M (ATC→ATGPP_2449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/20);  new base G (2/3);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

AGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACGGCGCAGCACGCTA  >  NC_002947/2797225‑2797516
                                                                                                                                                 |                                                                                                                                                  
aGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATccc                                                                                                                                                  <  7:342020/148‑1 (MQ=255)
       cGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATTCCCTCTTc                                                                                                                                            <  2:294467/147‑1 (MQ=255)
         gtgtAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAGAGGGGATGCCCTTTTCgatc                                                                                                                                        <  4:11808/149‑4 (MQ=255)
            gAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAg                                                                                                                                     >  2:344061/1‑149 (MQ=255)
                 gACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTAcc                                                                                                                                 <  6:109697/148‑1 (MQ=255)
                 gACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTAcc                                                                                                                                 >  5:109697/1‑148 (MQ=255)
                 gACCTCCAGAGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGACCAAGGGCTGGGAAGTGGAAGTGGTACATGGCGCCAACATGCGCGCACTCAGCGATATGACGGAGCGCCTGATCGAAAAGGGCATCACATCATCCGTGTAGTACCg                                                                                                                                >  6:126078/1‑149 (MQ=255)
                  aCCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGt                                                                                                                               >  4:239421/1‑149 (MQ=255)
                    cTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAGAGGGGATGCCCtttt                                                                                                                                             <  4:231929/133‑1 (MQ=255)
                    cTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAGAGGGGATGCCCtttt                                                                                                                                             >  3:231929/1‑133 (MQ=255)
                      gCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACg                                                                                                                           >  8:53479/1‑149 (MQ=255)
                                gCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGCAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCggt                                                                                                                 >  6:385462/1‑149 (MQ=255)
                                    aGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGc                                                                                                                     <  1:234990/141‑1 (MQ=255)
                                    aGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGc                                                                                                                     >  2:234990/1‑141 (MQ=255)
                                        gcTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTATAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGtgctgc                                                                                                         >  3:83798/1‑149 (MQ=255)
                                                      cGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTcccc                                                                                           <  5:365699/149‑1 (MQ=255)
                                                        gcgcCCCAGGGGTGGGATGTGGAAGTGGTATAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCAc                                                                                         <  4:83798/149‑1 (MQ=255)
                                                          gcCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCAccc                                                                                       >  8:190451/1‑149 (MQ=255)
                                                              cAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTg                                                                                    >  4:141011/1‑148 (MQ=255)
                                                                ggggTTGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTgaag                                                                                 <  7:53479/149‑1 (MQ=255)
                                                                 gggTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTgaaga                                                                                <  3:239421/149‑1 (MQ=255)
                                                                 gggTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTgaag                                                                                 >  3:110379/1‑148 (MQ=255)
                                                                        aTGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTATTCGGTGGTGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAg                                                                         >  5:294634/1‑149 (MQ=255)
                                                                         tgtgGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAgg                                                                        <  8:248280/149‑1 (MQ=255)
                                                                           tgGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGgcg                                                                      >  5:424565/1‑149 (MQ=255)
                                                                            ggAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGgcgc                                                                     >  6:304868/1‑149 (MQ=255)
                                                                            ggAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGgcgc                                                                     >  1:289427/1‑149 (MQ=255)
                                                                              aaGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGgcgcgc                                                                   <  7:190451/149‑1 (MQ=255)
                                                                               aGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCa                                                                  <  3:141011/149‑1 (MQ=255)
                                                                                gTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAg                                                                 <  5:304868/149‑1 (MQ=255)
                                                                                   gTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGa                                                              >  3:86560/1‑149 (MQ=255)
                                                                                   gTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGa                                                              >  4:278418/1‑149 (MQ=255)
                                                                                       aGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGaa                                                          <  3:363659/149‑1 (MQ=255)
                                                                                            gcCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGAt                                                     >  1:82533/1‑149 (MQ=255)
                                                                                            gcCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGAt                                                     <  4:86560/149‑1 (MQ=255)
                                                                                              cAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATgg                                                   >  6:149494/1‑149 (MQ=255)
                                                                                                       cgcACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTgg                                          <  6:424565/149‑1 (MQ=255)
                                                                                                              aGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCTCAGGGACGAAGTGATGGCCAAGCTGGGTACCg                                    <  5:385462/148‑1 (MQ=255)
                                                                                                              aGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGa                                   >  1:87309/1‑149 (MQ=255)
                                                                                                                   tatGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCAGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCgg                              >  6:385889/1‑149 (MQ=255)
                                                                                                                     tGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCggt                             <  1:344061/148‑1 (MQ=255)
                                                                                                                      gACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCggtgg                           >  3:340128/1‑149 (MQ=255)
                                                                                                                      gACGGCGCGCCTGATCGAAAAGGGCATCCCCCCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGaa                                                          >  4:405664/1‑118 (MQ=255)
                                                                                                                      gACGGCGCGCCTGATCGAAAAGGGCATCCCCCCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGaa                                                          <  3:405664/118‑1 (MQ=255)
                                                                                                                       aCGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCggtggt                          >  8:19557/1‑149 (MQ=255)
                                                                                                                       aCGGCGCGCCTGATCGAAAAGGGCATCCACTCTTCGGTGGAGTAACGTGACGGCAAGCCGGTGGTGCTGCTGTGGCCGTTCCCCACCCGTTAAGAAGAGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCggtggt                          >  8:13606/1‑149 (MQ=255)
                                                                                                                           cgcgcCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATc                      <  3:278418/149‑1 (MQ=255)
                                                                                                                             cgcCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTTGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCa                    <  4:110379/149‑1 (MQ=255)
                                                                                                                                    ccGAAGAGGGGATGCCCTTTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCgg                              >  4:309926/2‑132 (MQ=255)
                                                                                                                                    ccGAAGAGGGGATGCCCTTTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCgg                              <  3:309926/131‑1 (MQ=255)
                                                                                                                                     cGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACGgcg            >  8:198075/1‑149 (MQ=255)
                                                                                                                                      gAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATcc                     >  7:370083/1‑139 (MQ=255)
                                                                                                                                      gAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATcc                     <  8:370083/139‑1 (MQ=255)
                                                                                                                                      gAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACGgcg            >  8:119964/1‑148 (MQ=255)
                                                                                                                                       aaaaGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACGgcg            >  4:238941/1‑147 (MQ=255)
                                                                                                                                             gcgtCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACGGCGCAGCACGc    >  4:103016/4‑149 (MQ=255)
                                                                                                                                             gCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGCGGTCATCCAGCACGGCGCAGCACGc    <  5:149494/149‑1 (MQ=255)
                                                                                                                                              cATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACGGCGCAGCACGc    >  7:137078/1‑148 (MQ=255)
                                                                                                                                                 ccccTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACGGCGCAGCACGCTa  >  2:130767/1‑147 (MQ=255)
                                                                                                                                                 ccccTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACGGCGCAGCACGCTa  <  1:130767/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACGGCGCAGCACGCTA  >  NC_002947/2797225‑2797516

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: