Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,631,095 C→A 8.6% intergenic (+29/‑152) PP_1430 → / → lepA DegP‑like serine endoprotease/elongation factor 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,631,0950CA 8.6% 249.6 / 10.1 86intergenic (+29/‑152)PP_1430/lepADegP‑like serine endoprotease/elongation factor 4
Reads supporting (aligned to +/- strand):  ref base C (40/38);  new base A (3/5);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

GTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGTT  >  NC_002947/1630951‑1631241
                                                                                                                                                |                                                                                                                                                   
gTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCt                                                                                                                                                 >  5:71237/1‑149 (MQ=255)
gTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCt                                                                                                                                                 >  4:175356/1‑149 (MQ=255)
    ggCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGcc                                                                                                                                              <  2:36259/148‑1 (MQ=255)
       cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCttt                                                                                                                                          >  1:163921/1‑149 (MQ=255)
       cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCAAAAGGGCAGCCGAAGCTGCCCttt                                                                                                                                          <  4:89504/149‑1 (MQ=255)
        aaGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCtttt                                                                                                                                         >  8:283675/1‑149 (MQ=255)
        aaGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCtttt                                                                                                                                         >  6:237521/1‑149 (MQ=255)
        aaGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttc                                                                                                                                         <  7:198344/149‑2 (MQ=255)
         aGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCttttt                                                                                                                                        <  8:88900/149‑1 (MQ=255)
         aGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttcg                                                                                                                                        <  6:260533/149‑3 (MQ=255)
            aaTTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAATCATTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCAt                                                                                                                                     <  7:205153/149‑1 (MQ=255)
            aaTTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCAt                                                                                                                                     >  2:313064/1‑149 (MQ=255)
             aTTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttcgcctt                                                                                                                                    <  4:165352/149‑7 (MQ=255)
               tCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttc                                                                                                                                         <  8:157301/142‑2 (MQ=255)
               tCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttc                                                                                                                                         >  7:157301/1‑141 (MQ=255)
                  ccGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATGACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCt                                                                                                                               >  8:298352/1‑149 (MQ=255)
                   cGAAATCGCCAAGGAACTGCCGTAGAATCGTTCCGTGTGGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCATCTTCGGCTGCCCTTTTTCATGAATCtt                                                                                                                              <  4:157228/149‑1 (MQ=255)
                   cGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCtt                                                                                                                              <  6:311861/149‑1 (MQ=255)
                   cGAAATCGCCAAGGAACTCCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATGACTTTCAAGCTTGCTGAATAAGCACTTTCGCAAGGCGAAACGGCAGCTTCGGCTGACCTTTTTCATGAATCt                                                                                                                               >  3:157228/1‑148 (MQ=255)
                          gCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAAtt                                                                                                                      <  6:15368/149‑1 (MQ=255)
                           ccAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCACGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTGTTTCATGAATCTTCCCAAttt                                                                                                                     >  1:326008/1‑148 (MQ=255)
                            cAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCa                                                                                                                                      >  4:19351/1‑132 (MQ=255)
                            cAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCa                                                                                                                                      <  3:19351/132‑1 (MQ=255)
                               ggAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCgg                                                                                                                 <  4:200274/149‑1 (MQ=255)
                                      ccGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGAtt                                                                                                          >  6:266499/1‑149 (MQ=255)
                                       cGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTg                                                                                                         <  3:175356/149‑1 (MQ=255)
                                        gaagaaTCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGc                                                                                                        >  5:312687/1‑149 (MQ=255)
                                         aagaaTCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCt                                                                                                       <  5:237521/149‑1 (MQ=255)
                                                 ttCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAg                                                                                               >  6:142366/1‑149 (MQ=255)
                                                    cGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTc                                                                                             <  5:142366/148‑1 (MQ=255)
                                                      tgtCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGa                                                                                          >  3:295209/1‑149 (MQ=255)
                                                       gtCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGa                                                                                          >  2:290502/1‑148 (MQ=255)
                                                               gcgTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    <  5:217069/130‑1 (MQ=255)
                                                               gcgTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    <  5:221816/130‑1 (MQ=255)
                                                               gcgTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    >  6:221816/1‑130 (MQ=255)
                                                               gcgTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    >  6:217069/1‑130 (MQ=255)
                                                                     tGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACaga                                                                                      >  2:81081/1‑138 (MQ=255)
                                                                     tGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACaga                                                                                      <  1:81081/138‑1 (MQ=255)
                                                                     tGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCAGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCa                                                                           >  2:260869/1‑149 (MQ=255)
                                                                       cGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCAt                                                                          >  2:93798/1‑148 (MQ=255)
                                                                       cGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCAt                                                                          <  5:266499/148‑1 (MQ=255)
                                                                        gTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATaa                                                                        <  1:313064/149‑1 (MQ=255)
                                                                         tCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGTATGATTGCTGTGTACAGTCGGACAGAGGAATATCACATAAc                                                                       >  7:122158/1‑149 (MQ=255)
                                                                            gggCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAAcc                                                                      >  1:174197/1‑147 (MQ=255)
                                                                             ggCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcag                                                                   >  3:9070/1‑149 (MQ=255)
                                                                              gCCGTGCCAGCTTCATCTCTTTCAAGCTTGCTTAATAAGCAGTTTCTCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATTAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcagc                                                                  <  2:326008/149‑1 (MQ=255)
                                                                              gCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcag                                                                   >  1:201801/1‑148 (MQ=255)
                                                                               ccGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcagca                                                                 <  7:283675/149‑1 (MQ=255)
                                                                                 gTGCCATCTTCATCACTTTCAATCTTGCTGTATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTTCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCagcagc                                                               <  3:182247/149‑1 (MQ=255)
                                                                                      aGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAg                                                          <  4:295209/149‑1 (MQ=255)
                                                                                      aGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAg                                                          <  2:163921/149‑1 (MQ=255)
                                                                                              caCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCTGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATa                                                                         <  6:302499/126‑1 (MQ=255)
                                                                                              caCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCTGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATa                                                                         >  5:302499/1‑126 (MQ=255)
                                                                                               aCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTc                                                 >  1:245601/1‑149 (MQ=255)
                                                                                                 tttCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTccc                                               <  2:174197/149‑1 (MQ=255)
                                                                                                       gCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTa                                                       >  1:11603/1‑135 (MQ=255)
                                                                                                       gCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTa                                                       <  2:11603/135‑1 (MQ=255)
                                                                                                                 tAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTaa                                                           >  6:129764/1‑121 (MQ=255)
                                                                                                                 tAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTaa                                                           <  5:129764/121‑1 (MQ=255)
                                                                                                                 tAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcggcgg                               <  7:298352/149‑1 (MQ=255)
                                                                                                                   aGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCggcgggc                             >  6:139884/1‑149 (MQ=255)
                                                                                                                    gCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCggcgggcg                            >  6:166223/1‑149 (MQ=255)
                                                                                                                      aGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    <  5:62304/75‑1 (MQ=255)
                                                                                                                      aGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    >  6:62304/1‑75 (MQ=255)
                                                                                                                      aGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTAttttt                                     >  8:210661/1‑138 (MQ=255)
                                                                                                                      aGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTAttttt                                     <  7:210661/138‑1 (MQ=255)
                                                                                                                       gTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCGTTTTTCATGAATCTTACACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTc                          >  1:77613/1‑149 (MQ=255)
                                                                                                                       gTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTcc                         <  6:312687/149‑1 (MQ=255)
                                                                                                                        tttCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCg                        >  2:22265/1‑149 (MQ=255)
                                                                                                                            gCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTc                    <  1:93798/149‑1 (MQ=255)
                                                                                                                            gCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTc                    <  2:245601/149‑1 (MQ=255)
                                                                                                                             cAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCg                   >  5:269185/1‑149 (MQ=255)
                                                                                                                               aGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCa                 <  4:284749/149‑1 (MQ=255)
                                                                                                                               aGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCa                 <  1:290502/149‑1 (MQ=255)
                                                                                                                                tgaaaaaGGGCAGCCGCAGCTGCCCTTTTTCATGCATCTTCCCAATGGCCGGAATGATTGCTGTGGACAGGCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATGCCCGGCGAGTTTTCGGCGGGCGTCCGGCTCGccg                >  5:176795/5‑147 (MQ=255)
                                                                                                                                ggCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTACCGGCTATTTTTCGGCGGGCGTCCGGATCGCAg                >  7:153215/1‑149 (MQ=255)
                                                                                                                                 gaaaaaGGGCAGCCGAAGCTGGCCTTATTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCAAGCAGCTTAATGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTagcagc               >  7:178478/4‑149 (MQ=255)
                                                                                                                                 gaaaaaGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGc               >  3:108452/4‑149 (MQ=255)
                                                                                                                                   gAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCa                 >  7:131536/1‑145 (MQ=255)
                                                                                                                                   gAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCa                 <  8:131536/145‑1 (MQ=255)
                                                                                                                                       gggCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTc         <  6:269185/149‑1 (MQ=255)
                                                                                                                                         gCAGCTTCTGCTGCCCTTTTTCATGAATCTTCACAATTGCCGCAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGa       <  6:71237/149‑1 (MQ=255)
                                                                                                                                         gCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGa       >  7:215248/1‑149 (MQ=255)
                                                                                                                                           atcTTCGGCTGCCCTTTTTCATGAATCTTCACAATTCCCCTAATGATTGCTGTGTACAGTCTGACAGTGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCGGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAgt     <  7:304455/147‑1 (MQ=255)
                                                                                                                                           aGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAgt     <  2:77613/149‑1 (MQ=255)
                                                                                                                                              ttCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAAACCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGtt  >  7:122736/1‑149 (MQ=255)
                                                                                                                                                cGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGtt  <  8:215248/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGTT  >  NC_002947/1630951‑1631241

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