Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,631,094 T→G 10.6% intergenic (+28/‑153) 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,0940TG10.6% 241.2 / 18.5 85intergenic (+28/‑153)PP_1430/lepADegP‑like serine endoprotease/elongation factor 4
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base G (4/5);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

CATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGTT  >  NC_002947/1630946‑1631241
                                                                                                                                                    |                                                                                                                                                    
cATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCtt                                                                                                                                                      >  7:88900/1‑149 (MQ=255)
     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                                                                                                                                         >  7:157301/1‑141 (MQ=255)
                    tCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttc                                                                                                                                         <  8:157301/142‑2 (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                                                                                                    >  6:221816/1‑130 (MQ=255)
                                                                    gcgTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    >  6:217069/1‑130 (MQ=255)
                                                                    gcgTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    <  5:221816/130‑1 (MQ=255)
                                                                    gcgTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    <  5:217069/130‑1 (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                                                          <  2:163921/149‑1 (MQ=255)
                                                                                           aGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAg                                                          <  4:295209/149‑1 (MQ=255)
                                                                                                   caCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCTGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATa                                                                         >  5:302499/1‑126 (MQ=255)
                                                                                                   caCTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCTGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATa                                                                         <  6:302499/126‑1 (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                                                                                                    >  6:62304/1‑75 (MQ=255)
                                                                                                                           aGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgtg                                                                                                    <  5:62304/75‑1 (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                    <  2:245601/149‑1 (MQ=255)
                                                                                                                                 gCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTc                    <  1:93798/149‑1 (MQ=255)
                                                                                                                                  cAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCg                   >  5:269185/1‑149 (MQ=255)
                                                                                                                                    aGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCa                 <  1:290502/149‑1 (MQ=255)
                                                                                                                                    aGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCa                 <  4:284749/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                 <  8:131536/145‑1 (MQ=255)
                                                                                                                                        gAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCa                 >  7:131536/1‑145 (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)
                                                                                                                                                    |                                                                                                                                                    
CATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAATAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGTT  >  NC_002947/1630946‑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: