Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,068,200 G→A 17.4% intergenic (+65/‑142) 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,2000GA17.4% 552.2 / 93.7 236intergenic (+65/‑142)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (94/97);  new base A (14/25);  total (108/128)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-02

Only 100 of 248 total aligned reads displayed.

AAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGA  >  NC_002947/2068052‑2068349
                                                                                                                                                      |                                                                                                                                                     
aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGa                                                                                                                                                      >  2:2645690/1‑150 (MQ=255)
aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAGTGCACAGGGTTCGGa                                                                                                                                                      >  2:5786568/1‑150 (MQ=255)
  gaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAATGCACCGGGTTCGGAcc                                                                                                                                                    >  1:4368085/1‑150 (MQ=255)
      cTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgt                                                                                                                                                <  1:2420528/150‑1 (MQ=255)
      cTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgt                                                                                                                                                <  1:7326748/150‑1 (MQ=255)
       tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCCCCGGCAATGCACCGGGGTCGGACCAtgtg                                                                                                                                               >  1:3562943/1‑150 (MQ=255)
       tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                                               <  1:180729/150‑1 (MQ=255)
          aaGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAg                                                                                                                                            <  1:4492636/150‑1 (MQ=255)
            gCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGG‑‑GGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatg                                                                                                                                          <  1:6132469/150‑4 (MQ=21)
            gCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGt                                                                                                                                          <  1:7183418/150‑1 (MQ=255)
             cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                                                         <  1:6502619/150‑5 (MQ=25)
               tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGTACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCt                                                                                                                                       <  1:1431464/150‑1 (MQ=255)
               tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                                                         <  1:1189781/150‑5 (MQ=18)
               tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgcct                                                                                                                                       <  2:1687258/150‑7 (MQ=21)
               tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgcct                                                                                                                                       <  2:237520/150‑7 (MQ=21)
               tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCt                                                                                                                                       <  1:6991651/150‑1 (MQ=255)
               tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCt                                                                                                                                       <  1:6355743/150‑1 (MQ=255)
                atatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCGTGTGTAGGTCCtg                                                                                                                                      >  2:2031176/1‑150 (MQ=255)
                 tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCTGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                                     <  1:2883506/150‑1 (MQ=255)
                 tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                                     <  1:2796887/150‑1 (MQ=255)
                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtgg                                                                                                                                   <  1:3544641/150‑2 (MQ=21)
                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtgg                                                                                                                                   <  2:7308875/150‑2 (MQ=21)
                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                   <  2:5430578/150‑1 (MQ=255)
                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                   >  2:229969/1‑150 (MQ=255)
                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                   <  2:2740195/150‑1 (MQ=255)
                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                   >  2:7513104/1‑150 (MQ=255)
                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                   <  1:5458152/150‑1 (MQ=255)
                    ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAGTGCACCGGGTTCGGACCATGTGCAGGTCCAGTGCg                                                                                                                                  >  2:2568704/1‑150 (MQ=255)
                      aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTTCACAGGGTTCGGACCATGTGCAGGTCCAGTGCGAg                                                                                                                                >  2:5999241/1‑150 (MQ=255)
                       gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGc                                                                                                                               <  2:7937718/150‑1 (MQ=255)
                         ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                             <  1:7225535/150‑1 (MQ=255)
                         ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                             <  1:5483516/150‑1 (MQ=255)
                         ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCGTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                             <  2:1145439/150‑1 (MQ=255)
                             gACCAAGAGCCTGCTGGATAACACTGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCaa                                                                                                                         <  2:7927787/150‑1 (MQ=255)
                             gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCGTGCCCCCGCCCACCCGCCGTGCCCAGGGGTCGGACCCTGTGCAGGTCCGGTGCGAGCGGGccc                                                                                                                         >  2:1927005/1‑148 (MQ=255)
                                   gagCCTGCTGGATACCACCGCCGGCGGGGTGCACGGCACACCGACGTAACCCGGGGCCGCGTTGC‑‑GGCCAATTCGCGGGCTTGCCCGCTCCCACCGGCTTTGCACACGGTCCGACCCCGGTGCAAGGCCTGTGGGAGCGGGCAAGCCcgc                                                                                                                   <  1:615815/150‑1 (MQ=21)
                                       ctgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGAAAGCCCGCGAAt                                                                                                               <  2:7189743/150‑1 (MQ=255)
                                           gggATACCACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCTGTTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGc                                                                                                           <  2:6202632/149‑1 (MQ=255)
                                            ggATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGGTTGCCCGCTCCCCAAGGCAATGCACAGGGTTCGGACCATTTGCAGGTCCTGTGCGAGCGGGCAAGCCCCCGAATGGGcc                                                                                                          >  2:3196477/1‑150 (MQ=255)
                                                      gccgGCGTGGGGCACGCCACCCCGACGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCGTGCCCGCCCCCCCCGGAATTGCACAGGGTCCGGACCTTGTGCAGGTCCGGTGCGAGGGGGCAAGCCCGCGAATGGGCCCCGCACAAtt                                                                                                <  2:3115245/150‑1 (MQ=25)
                                                         ggCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCATGTGCAGGTCCTGTGCGAGCGGGCCAGCCCGCGAATGGGCCCCGCACAATTTTc                                                                                             >  2:7652226/1‑150 (MQ=255)
                                                           cGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGTTTCGGACCTTGTGATGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTccc                                                                                           <  2:1221491/150‑1 (MQ=255)
                                                                  cACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCGTGCCCGCGCCCACAGGCATTTCACAGGGTTCGGACCATTTGCAGGTCCAGTGCGAGCGGGCCACCCCCCCAATGGGGCCCCCACAAATTTACCCGACGGa                                                                                    >  2:2994438/1‑150 (MQ=255)
                                                                   aCGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGaa                                                                                   <  2:8282695/150‑1 (MQ=255)
                                                                       cacaCCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCTTTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACAtt                                                                               <  2:6940620/150‑1 (MQ=255)
                                                                           ccGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATcc                                                                           <  2:8727300/150‑1 (MQ=255)
                                                                              aaGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTTCCCCCAACCCCAGGCAATGCACATGGTTGCGACCTTGTTCCGGGACAGTTCGAGCGGGGAATCCCGCCAATGGGGCCCGCAAAATTTTACCCGACAGAACATTATACCCg                                                                        >  2:3996457/1‑150 (MQ=255)
                                                                                    gaatgggCCGCAACGC‑‑GGCCCCTTCGCGGGCTTGCCCGCTCCCCCAGGCATTGCCCCGGGTTCGGACCTTGGGCAGGTCCGGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGCCGGAACAATATCccccccccca                                                                  >  2:4924767/5‑148 (MQ=25)
                                                                                         gcccGCGTTGC‑‑GCCCCATCCGCGGGCTTCCCCGCACCCCCCGCAATTCAACAGGGTCCGCACCTTTTCCAGTACCTGTGCCAGCGCCAAACCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                                             <  1:5266600/148‑1 (MQ=255)
                                                                                               gTTGG‑‑GGCCCTTTCGGGGGCTTGCCCTCCCCCACAGCCATTGCCCCGGGTTCGCAGCTGGTGCCGGCCGGGGGCGCGCGGGAAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTAtt                                                       <  1:7503721/150‑1 (MQ=255)
                                                                                                ttGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATGGCACTGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTc                                                      <  2:7469350/150‑1 (MQ=255)
                                                                                                 tGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                     <  2:4330367/150‑1 (MQ=255)
                                                                                                 ggC‑‑GGCCCTTCCGCGGGCTGGCCCGCTCCCCCCGGCTTGGCACCGGGTTCGGCCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                     <  1:5757024/149‑1 (MQ=255)
                                                                                                         ccATTCGCGGGCTTGCCCGCTCCCACGGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCCAGGTCAGTTATTCAGGCATc                                               >  2:1109487/1‑150 (MQ=255)
                                                                                                          cttgcgggggCTTGCCCGCTCCCCCCGGCATTCCAGTGGGTTCGGACCTTGTGCCGGCCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATca                                              <  1:4293497/143‑1 (MQ=255)
                                                                                                            ttCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATcacc                                            >  2:6504194/1‑149 (MQ=255)
                                                                                                              cgcgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATcacaca                                          <  2:4851525/150‑1 (MQ=255)
                                                                                                                cgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGGAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACAcc                                        <  2:5097483/150‑1 (MQ=255)
                                                                                                                 gggCTTGCCCGCTCCCACAGGCATTGCACAGGGTTGGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCt                                       <  1:5789446/150‑1 (MQ=255)
                                                                                                                    ggtGCCCGCGCCCACCGGCTTTGCACCGGGTTCGGACGTTGTGCAGGTCCTGGGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt                                    <  1:7911451/148‑1 (MQ=255)
                                                                                                                    cTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt                                    >  2:8208396/1‑150 (MQ=255)
                                                                                                                        cccGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                <  2:2464958/150‑1 (MQ=255)
                                                                                                                         ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                               >  1:364650/1‑150 (MQ=255)
                                                                                                                         ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                               >  2:5596462/1‑150 (MQ=255)
                                                                                                                         ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                               >  2:1898714/1‑150 (MQ=255)
                                                                                                                          cGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGCCTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt                              >  2:3968986/1‑150 (MQ=255)
                                                                                                                          cGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt                              >  1:2593744/1‑150 (MQ=255)
                                                                                                                           gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                             >  1:6803461/1‑150 (MQ=255)
                                                                                                                           gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                             >  2:4552056/1‑150 (MQ=255)
                                                                                                                           gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                             >  2:15042/1‑150 (MQ=255)
                                                                                                                               ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                         >  2:1772040/1‑150 (MQ=255)
                                                                                                                               ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                         >  2:6644897/1‑150 (MQ=255)
                                                                                                                               ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                         <  2:4602013/150‑1 (MQ=255)
                                                                                                                               ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                         <  2:3429726/150‑1 (MQ=255)
                                                                                                                                cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTg                        >  2:3487966/1‑150 (MQ=255)
                                                                                                                                cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTg                        >  2:2347747/1‑150 (MQ=255)
                                                                                                                                  caGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgag                      >  2:5274077/1‑150 (MQ=255)
                                                                                                                                  caGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTACGTTGACagag                      >  2:6601702/1‑150 (MQ=255)
                                                                                                                                   aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                     >  2:7736686/1‑150 (MQ=255)
                                                                                                                                    ggacctGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCGCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGagag                    >  1:5199875/6‑150 (MQ=21)
                                                                                                                                     gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg                   <  1:5817542/150‑1 (MQ=255)
                                                                                                                                      acctGCACAAGGTCCGAACCCTGTGCAAGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAggg                  >  1:4627717/4‑150 (MQ=37)
                                                                                                                                          gCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCg              >  2:8549447/1‑150 (MQ=255)
                                                                                                                                           cacaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGt             >  2:6980112/1‑150 (MQ=255)
                                                                                                                                            acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt            <  2:3287525/150‑1 (MQ=255)
                                                                                                                                              aaggTCCGAACCCTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGaa                                                                                                                <  1:3764964/46‑1 (MQ=11)
                                                                                                                                              aaggTCCGAACCCTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGaa                                                                                                                >  2:3764964/3‑48 (MQ=11)
                                                                                                                                              aaggTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCa          >  1:4975312/3‑150 (MQ=37)
                                                                                                                                                ggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt        >  2:2704477/1‑150 (MQ=255)
                                                                                                                                                ggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt        >  2:107580/1‑150 (MQ=255)
                                                                                                                                                ggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt        <  2:3612747/150‑1 (MQ=255)
                                                                                                                                                ggTTCGCACGTTGGGCTGGGCGGGGGCGAGCGGCCACCCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt        <  1:1518411/150‑1 (MQ=255)
                                                                                                                                                ggTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt        <  1:3586574/150‑1 (MQ=37)
                                                                                                                                                  ttCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGa      >  1:2273265/1‑150 (MQ=255)
                                                                                                                                                  ttCGGACATTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGATTGTCCGTTCAATGa      <  2:6724067/150‑1 (MQ=255)
                                                                                                                                                    cGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGaaa    <  2:465339/150‑1 (MQ=255)
                                                                                                                                                      gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga  >  1:170045/1‑150 (MQ=255)
                                                                                                                                                      gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga  >  2:1331380/1‑150 (MQ=255)
                                                                                                                                                      gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga  >  1:493003/1‑150 (MQ=255)
                                                                                                                                                      gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga  >  1:8014065/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGC‑‑GGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGA  >  NC_002947/2068052‑2068349

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: