Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,609,691 G→A 6.3% A275V (GCC→GTC)  PP_4929 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,609,6910GA 6.3% 190.0 / 2.0 63A275V (GCC→GTC) PP_4929LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/2);  total (30/33)
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.98e-01

GCGCAACCAGCGGCAGGCCCAGGTGAAACAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGT  >  NC_002947/5609543‑5609838
                                                                                                                                                    |                                                                                                                                                   
gcgcaACCAGCGGCAGGCCCAGGTGAAACAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCgccagg                                                                                                                                                     >  7:262217/1‑149 (MQ=255)
     aCCAGCGGCAGGCCCAGGTGAAACAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGc                                                                                                                                                <  2:295972/149‑1 (MQ=255)
        aGCGGCAGGCCCAGGTGAAACAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGctct                                                                                                                                             <  3:163838/149‑1 (MQ=255)
         gCGGCAGGCCCAGGTGAAACAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTg                                                                                                                                            <  8:262217/149‑1 (MQ=255)
                  ccAGGTGAAACAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAgg                                                                                                                                   >  4:158423/1‑149 (MQ=255)
                       tGAAACAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCcggc                                                                                                                              >  4:241479/1‑149 (MQ=255)
                        gAAACAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCcggcg                                                                                                                             >  7:147035/1‑149 (MQ=255)
                         aaaCAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCcggcgg                                                                                                                            <  3:53050/149‑1 (MQ=255)
                               cAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCa                                                                                                                      <  8:72368/149‑1 (MQ=255)
                                  ggTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTc                                                                                                                   >  8:129042/1‑149 (MQ=255)
                                     ataGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCACt                                                                                                                >  3:265769/1‑149 (MQ=255)
                                      taGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCt                                                                                                                <  8:100517/148‑1 (MQ=255)
                                        gagagTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTc                                                                                                             <  7:129042/149‑1 (MQ=255)
                                        gagagGTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTc                                                                                                             <  3:59203/149‑1 (MQ=255)
                                         agagTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTc                                                                                                             >  1:207486/1‑148 (MQ=255)
                                                        gcggcTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCAtgtg                                                                                              <  8:147035/148‑1 (MQ=255)
                                                        gcggcTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCAAGTTCGACCATGTgc                                                                                             <  2:156795/149‑1 (MQ=255)
                                                           gcTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGcc                                                                                          >  5:229000/1‑149 (MQ=255)
                                                            cTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCt                                                                                         <  6:221317/149‑1 (MQ=255)
                                                               aGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGt                                                                                      <  3:158423/149‑1 (MQ=255)
                                                                     cccGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGt                                                                                      <  5:289555/143‑1 (MQ=255)
                                                                     cccGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGt                                                                                      >  6:289555/1‑143 (MQ=255)
                                                                      ccGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTc                                                                                 <  2:172384/147‑1 (MQ=255)
                                                                      ccGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTc                                                                                 >  1:172384/1‑147 (MQ=255)
                                                                       cGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCggg                                                                              <  2:100106/149‑1 (MQ=255)
                                                                          tCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTAc                                                                           >  3:42923/1‑149 (MQ=255)
                                                                             cgcgAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGc                                                                        <  4:123748/149‑1 (MQ=255)
                                                                              gcgAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCt                                                                                                                >  2:159801/1‑108 (MQ=255)
                                                                              gcgAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCt                                                                                                                <  1:159801/108‑1 (MQ=255)
                                                                                gAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGcac                                                                      >  5:312378/1‑148 (MQ=255)
                                                                                   ggTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACcaca                                                                  >  5:277577/1‑149 (MQ=255)
                                                                                      gCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCt                                                                                         <  4:285822/123‑1 (MQ=255)
                                                                                      gCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCt                                                                                         >  3:285822/1‑123 (MQ=255)
                                                                                      gCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGACCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGc                                                               >  2:109691/1‑149 (MQ=255)
                                                                                      gCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGACCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGc                                                               >  2:113327/1‑149 (MQ=255)
                                                                                         cTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCgg                                                            <  2:207486/149‑1 (MQ=255)
                                                                                              aaGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTaa                                                       >  7:11298/1‑149 (MQ=255)
                                                                                              aaGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTaa                                                       >  7:16030/1‑149 (MQ=255)
                                                                                                  gcTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGACCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAg                                                   <  1:109691/149‑1 (MQ=255)
                                                                                                  gcTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGACCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAg                                                   <  1:113327/149‑1 (MQ=255)
                                                                                                   cTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAgg                                                  >  4:127101/1‑149 (MQ=255)
                                                                                                   cTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAgg                                                  >  3:203821/1‑149 (MQ=255)
                                                                                                   cTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAgg                                                  >  2:25867/1‑149 (MQ=255)
                                                                                                    tCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGcac                                                                      <  7:130783/128‑1 (MQ=255)
                                                                                                    tCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGcac                                                                      >  8:130783/1‑128 (MQ=255)
                                                                                                      ggCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGccc                                               <  3:241479/149‑1 (MQ=255)
                                                                                                        ccagccaCTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCggg                                                           <  4:26556/135‑1 (MQ=255)
                                                                                                        ccagccaCTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCggg                                                           >  3:26556/1‑135 (MQ=255)
                                                                                                         cagccaCTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGc                                            <  4:203821/149‑1 (MQ=255)
                                                                                                                  cgcgGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCcgc                                   >  8:74958/1‑149 (MQ=255)
                                                                                                                  cgcgGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGGCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCcg                                    <  6:229000/148‑1 (MQ=255)
                                                                                                                      gCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCAt                               <  3:127101/149‑1 (MQ=255)
                                                                                                                          ggCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATcagc                           <  6:312378/149‑1 (MQ=255)
                                                                                                                             ccAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGtt                                                                                                        >  6:327982/1‑69 (MQ=255)
                                                                                                                             ccAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTGGGTCGGGTACTGCACCACATGCCgg                                                            >  6:215252/1‑113 (MQ=255)
                                                                                                                             ccAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTGGGTCGGGTACTGCACCACATGCCgg                                                            <  5:215252/113‑1 (MQ=255)
                                                                                                                             ccAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGCTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGt                        >  2:110172/1‑149 (MQ=255)
                                                                                                                               agagGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCcgcg                                  >  4:239101/1‑137 (MQ=255)
                                                                                                                               agagGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCcgcg                                  <  3:239101/137‑1 (MQ=255)
                                                                                                                                   gTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCgg                                                                               <  1:26333/88‑1 (MQ=255)
                                                                                                                                   gTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCgg                                                                               >  2:26333/1‑88 (MQ=255)
                                                                                                                                       aTCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATg               <  6:277577/148‑1 (MQ=255)
                                                                                                                                       aTCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGCTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATg               <  1:110172/148‑1 (MQ=255)
                                                                                                                                              cgcCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGa        >  3:49819/1‑148 (MQ=255)
                                                                                                                                                 caggccagCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGaaa      <  5:302839/147‑1 (MQ=255)
                                                                                                                                                 caggccagCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACt    >  3:253815/1‑149 (MQ=255)
                                                                                                                                                   ggccagCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGt  <  6:281959/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGCAACCAGCGGCAGGCCCAGGTGAAACAGCAAGGTATAGAGAGTTCTGTTCATGGCGGCTAAGTCTACCCGATCGCGCGAAGGTGCACTGCAAAGCGCTCGGCCAGCCACTGCGCGGCCGGGCCCAGAGGTTCATCGCGGCGCCAGGCCAGCTCTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGT  >  NC_002947/5609543‑5609838

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑