Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,068,197 T→C 18.7% intergenic (+62/‑145) PP_1838 → / → PP_1839 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,068,1970TC18.7% 472.0 / 90.7 203intergenic (+62/‑145)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (86/79);  new base C (21/17);  total (107/96)
Fisher's exact test for biased strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-02

Only 100 of 209 total aligned reads displayed.

GCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGA  >  NC_002947/2068048‑2068345
                                                                                                                                                     |                                                                                                                                                    
gCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGtt                                                                                                                                                      <  2:4616218/150‑1 (MQ=255)
 cTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTc                                                                                                                                                     <  2:1020010/150‑1 (MQ=255)
 cTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTc                                                                                                                                                     <  1:3613913/150‑1 (MQ=255)
  tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCg                                                                                                                                                    >  2:3555446/1‑150 (MQ=38)
  tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCg                                                                                                                                                    >  2:4075834/1‑150 (MQ=255)
    aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGa                                                                                                                                                  >  2:1803638/1‑150 (MQ=255)
    aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCCGTGCACAGGGGTCGGa                                                                                                                                                  >  2:2032892/1‑150 (MQ=255)
     agaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAAc                                                                                                                                                 <  1:4407321/150‑1 (MQ=37)
        agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAAccct                                                                                                                                              <  2:5147583/150‑3 (MQ=37)
        agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                              <  1:1876342/150‑1 (MQ=255)
        agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                              <  1:2459536/150‑1 (MQ=255)
           tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                                           <  1:1361407/150‑1 (MQ=255)
           tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                                           <  1:4203125/150‑1 (MQ=255)
            ggAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGc                                                                                                                                          >  2:4145352/1‑150 (MQ=255)
            ggAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATGCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGc                                                                                                                                          <  2:4809880/150‑1 (MQ=255)
             gAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCa                                                                                                                                         <  1:1087332/150‑1 (MQ=255)
               aGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAgg                                                                                                                                       <  1:1561734/150‑1 (MQ=255)
                gCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGt                                                                                                                                      <  2:4534988/150‑1 (MQ=255)
                  gTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGGCCGCTCCCACAGGCATTTGACAGGGTTGGGGCCTGGGGCAGGTcc                                                                                                                                    >  2:1442898/1‑150 (MQ=255)
                   tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCt                                                                                                                                   <  2:1140108/150‑1 (MQ=255)
                   tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgcct                                                                                                                                   <  1:802706/150‑7 (MQ=25)
                   tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAATCGCGGGCGTGCCCCCCCCCAAAGGCCGAGCACCGGGTTGGGGCCCAGTGCCGGTccg                                                                                                                                   >  1:4358870/1‑149 (MQ=255)
                    atatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtg                                                                                                                                  <  1:658272/150‑1 (MQ=255)
                    atatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtg                                                                                                                                  <  1:2375156/150‑1 (MQ=255)
                       tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                               <  1:3963686/150‑1 (MQ=255)
                        ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCg                                                                                                                              >  2:3747649/1‑150 (MQ=255)
                         cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGa                                                                                                                             <  1:1373567/150‑1 (MQ=18)
                         cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGa                                                                                                                             <  1:863466/150‑1 (MQ=255)
                         cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGa                                                                                                                             <  1:5463743/150‑1 (MQ=255)
                         cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGa                                                                                                                             <  1:3875871/150‑1 (MQ=21)
                         cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCGCCCACACGCCGTGCCCAGGGGTCGGACCATGTGCAGGTCCAGTGCGa                                                                                                                             >  2:4053315/1‑150 (MQ=255)
                         cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTAGCGGGCATGCCAGCACCCACAGGCCATGCACAGGGTTAGGACCATTTGCAGGTACAGTGCGa                                                                                                                             >  2:606261/1‑150 (MQ=255)
                          aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAg                                                                                                                            <  1:3555446/150‑1 (MQ=18)
                          aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACACGCAATGCACAGGGGTCGGGCCTTGTGGAGGTGCTGTGCggg                                                                                                                            >  2:2060841/1‑148 (MQ=39)
                           gcccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGc                                                                                                                           <  2:2124944/148‑1 (MQ=255)
                             ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                         <  2:846953/150‑1 (MQ=18)
                             ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                         <  2:1243392/150‑1 (MQ=18)
                              cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCggg                                                                                                                        <  1:5070503/150‑1 (MQ=255)
                                 gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCaa                                                                                                                     <  2:426228/150‑1 (MQ=255)
                                  aCCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGAcc                                                                                                                                                >  1:3398159/1‑122 (MQ=255)
                                          ccTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTAGCGGGCGTGCCCCCTCCCCCAGGCAATGCACAGGGGTCGGGCCATGTGGAGGTCCAGTGCGAGCGGGCAAGccccccaa                                                                                                            >  1:928231/1‑145 (MQ=255)
                                            ggctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATg                                                                                                          <  2:3380137/149‑1 (MQ=255)
                                            ggctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATGGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATg                                                                                                          <  2:888053/149‑1 (MQ=255)
                                                   tAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAg                                                                                                                                        <  1:4507113/113‑1 (MQ=255)
                                                        ccgccgGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCGCCCCCAGGCAATGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCCCGAAAGGGCCCCGCACaa                                                                                              >  2:4771490/1‑150 (MQ=255)
                                                        ccgccgGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCCCGCCCACAGGCCGTGCACAGGGTTCGGGCCGGTTGGAGGTCCAGTGGGAGCGGGCAAGCCCGCGCAGGGGCCCCCcacca                                                                                              >  1:1307126/1‑148 (MQ=255)
                                                         cgccgGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCCCTCCCACAGGCATTGCCCAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAACCCCGCGAATGGGCCCCCCACAAt                                                                                             >  2:3879757/1‑150 (MQ=255)
                                                          gccgGCGTGGGGCACGGCACACCGACGTAACCCGGGGCCGGTTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCTTTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCACGCCCGCGAATGGGCCCCGCACAAtt                                                                                            <  1:5580855/150‑1 (MQ=37)
                                                            cgGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTGGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAAtttt                                                                                          <  2:926151/150‑1 (MQ=255)
                                                                   gtgCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCCTTGCACCGGGTTCGGACCTTGTGCAGGGCCTGGGCGAGCGGGCAAGCCCGCGAAGGGGCCCCGCACAAATTTCCCgggc                                                                                   >  2:3155619/1‑148 (MQ=37)
                                                                           cacaCCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTAGCGGGCTTGCCCGCTCCCACACGCCTTGCACAGGGTTGGGGACGATTTCAGGTGCTGTGCGGGCGGGCAACCCCGCGAAAGGGCCCCCCACAAATTTACCGGACTGaaaaat                                                                           >  2:658449/1‑146 (MQ=255)
                                                                              acCTAAGAAAACCGGGGCCGCGGTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCCCCGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATc                                                                        <  1:2267355/150‑1 (MQ=38)
                                                                              acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGGTTGCCCGCTCCCACCGGCCGTGCCCACGGGTCGGTCCTTTTGCAGGTCCGGTGCGGGCGGGCCAGCCCGCGCATGGGCCCCGCACCAATTTCCCGGACTGAACATTaac                                                                        >  2:4295024/1‑148 (MQ=255)
                                                                              acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTTTCCGCTCCCCCAGGCCTTGCACAGGGTTAGGACCTTGTGCAGGTCCGGGGCGAGCGGGCCAGCCCGCGCAGGGGCCCCCCACAAATTTCCCGGACAGAGCAGTaac                                                                        >  2:5167436/1‑148 (MQ=255)
                                                                                cGAAGTAACCCGGGGCCGCGTTGGGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATccc                                                                      <  2:5361216/150‑1 (MQ=255)
                                                                                      acaCCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCGGCTCCCACCGGAATTGCACGGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGcccc                                                                <  1:4686815/148‑1 (MQ=255)
                                                                                            gggCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTc                                                          <  2:4591001/150‑1 (MQ=255)
                                                                                                cgcgTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGtt                                                      <  2:4353732/150‑1 (MQ=255)
                                                                                                   gTTGCGGCCCATGCGGGGGCTGGCCCGCTCCCCCCGCCATGGCCCCGGGTTCGGACCTTGTGCAGGTCCTGGGCGCGGGCCCAACCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTAtt                                                   <  1:1717412/150‑1 (MQ=255)
                                                                                                     tGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                 <  2:4677177/150‑1 (MQ=255)
                                                                                                      gCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAg                                                <  2:4569513/150‑1 (MQ=255)
                                                                                                                 gcgGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAAATTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCacacac                                     >  1:3278366/1‑150 (MQ=17)
                                                                                                                 gcgGGCTTGCCCGCTCCCAAAGGCATTGCACACGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAAATTTCCCGGAATGAACATTATCCCCGCCCCATGGTCAGTTATTCAGGCATCacacac                                     >  2:2343646/1‑150 (MQ=38)
                                                                                                                  cgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACAcc                                    >  2:354618/1‑150 (MQ=255)
                                                                                                                   gggCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCt                                   >  2:3207991/1‑150 (MQ=255)
                                                                                                                    ggCTTGCCCGCTCCCACAGGCTTTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTg                                  <  1:4162532/150‑1 (MQ=255)
                                                                                                                       ttGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtg                               >  1:3835216/1‑150 (MQ=255)
                                                                                                                          cccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                            >  2:4125581/1‑150 (MQ=255)
                                                                                                                           ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                           >  1:2114619/1‑150 (MQ=18)
                                                                                                                            cGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt                          >  2:4852639/1‑150 (MQ=18)
                                                                                                                            cGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt                          <  2:3835216/150‑1 (MQ=255)
                                                                                                                             gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                         >  2:2942625/1‑150 (MQ=18)
                                                                                                                             gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                         >  2:3527577/1‑150 (MQ=18)
                                                                                                                             gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                         >  2:5400426/1‑150 (MQ=255)
                                                                                                                                 ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCaa                                                              <  2:2481918/109‑1 (MQ=255)
                                                                                                                                 ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                     >  2:2202896/1‑150 (MQ=255)
                                                                                                                                  cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTg                    >  1:2736799/1‑150 (MQ=255)
                                                                                                                                   acaGGCATTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                   >  2:2422048/1‑150 (MQ=255)
                                                                                                                                   acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                   >  1:450877/1‑150 (MQ=255)
                                                                                                                                   acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                   >  1:5437051/1‑150 (MQ=255)
                                                                                                                                   acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                   >  1:3652281/1‑150 (MQ=255)
                                                                                                                                    caGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgag                  >  2:1647196/1‑150 (MQ=255)
                                                                                                                                    caGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgag                  >  1:1123772/1‑150 (MQ=255)
                                                                                                                                     cggCATTGCGCTGGGTTCGGACCTGGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                 <  2:2867437/149‑1 (MQ=255)
                                                                                                                                     aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                 >  2:170007/1‑150 (MQ=255)
                                                                                                                                         aTTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGt             >  1:3738080/1‑150 (MQ=255)
                                                                                                                                             cacaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGt         <  2:5404600/150‑1 (MQ=255)
                                                                                                                                             cacaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGt         >  1:1023664/1‑150 (MQ=255)
                                                                                                                                             cacaAGGTCCGAACCCTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                         >  1:1524091/1‑134 (MQ=255)
                                                                                                                                             cacaAGGTCCGAACCCTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                         <  2:1524091/134‑1 (MQ=255)
                                                                                                                                              acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt        >  2:4875756/1‑150 (MQ=255)
                                                                                                                                              acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATGTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt        <  2:4622551/150‑1 (MQ=255)
                                                                                                                                               caGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                         >  1:2535206/1‑132 (MQ=255)
                                                                                                                                               caGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                         <  2:2535206/132‑1 (MQ=255)
                                                                                                                                                 gggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCaa     >  1:2429986/1‑150 (MQ=255)
                                                                                                                                                  ggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACCCACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt    <  2:2614619/150‑1 (MQ=255)
                                                                                                                                                  ggTCCGAACCCTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCACGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt    >  2:4058547/1‑150 (MQ=255)
                                                                                                                                                  ggTCCGAACCCTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt    >  2:5318747/1‑150 (MQ=255)
                                                                                                                                                  ggTCCGAACCCTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt    >  2:4648478/1‑150 (MQ=255)
                                                                                                                                                   gTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATg   >  1:3134654/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGA  >  NC_002947/2068048‑2068345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: