Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,065,150 C→G 13.1% intergenic (+42/‑165) 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,065,1500CG13.1% 271.2 / 27.6 92intergenic (+42/‑165)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/48);  new base G (6/6);  total (38/54)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CCGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTG  >  NC_002947/2065004‑2065296
                                                                                                                                                  |                                                                                                                                                  
ccGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCcac                                                                                                                                                 <  2:3906118/150‑1 (MQ=255)
  ggACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGAACGCGACCAAGAGCCAGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGAGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAg                                                                                                                                               >  2:4616413/1‑150 (MQ=255)
  ggACCAGCCGGAAATGCTAAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGAGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAg                                                                                                                                               >  2:4616367/1‑150 (MQ=255)
        gCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTg                                                                                                                                         >  2:297837/1‑150 (MQ=255)
        gCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAATg                                                                                                                                         >  2:638262/1‑150 (MQ=255)
        gCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCGCCCACAGGCAATg                                                                                                                                         >  2:638152/1‑150 (MQ=255)
           ggAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcac                                                                                                                                      <  2:3173688/150‑1 (MQ=255)
            gAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcaca                                                                                                                                     <  1:2069927/150‑1 (MQ=255)
            gAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcaca                                                                                                                                     <  1:2489629/150‑1 (MQ=255)
            gAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcaca                                                                                                                                     >  2:3004548/1‑150 (MQ=255)
             aaaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAg                                                                                                                                    <  2:4372110/150‑1 (MQ=255)
              aaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAgg                                                                                                                                   <  1:1027978/150‑1 (MQ=255)
              aaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAgg                                                                                                                                   <  2:4652932/150‑1 (MQ=255)
              aaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAgg                                                                                                                                   <  2:4652953/150‑1 (MQ=255)
              aaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAgg                                                                                                                                   <  1:4616271/150‑1 (MQ=255)
                 gCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGtt                                                                                                                                <  1:1499523/150‑1 (MQ=255)
                 gCTCAAGAAGCTGGAAGCGTATATCCAGACCGCCACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGtt                                                                                                                                <  2:3232680/150‑1 (MQ=255)
                  cTCAATAAGCTGGAAGCGTCTATCCAGACCGCGACCAAGAGCCTGCTGTATAACACCTCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCTCCCACAGGCATTGCACAGGGTTc                                                                                                                               <  1:1095389/150‑1 (MQ=255)
                  cTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTc                                                                                                                               >  1:3549583/1‑150 (MQ=255)
                  cTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTc                                                                                                                               <  1:4599475/150‑1 (MQ=255)
                     aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCCATGCACACGGTTCGGa                                                                                                                            >  1:5039410/1‑150 (MQ=255)
                     aagaagCTGGAAGCGTATACCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGa                                                                                                                            <  2:2384943/150‑1 (MQ=255)
                      agaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACAGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGAc                                                                                                                           <  2:2581469/150‑1 (MQ=255)
                         agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                        <  1:40894/150‑1 (MQ=255)
                         agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                        <  2:1321971/150‑1 (MQ=255)
                         agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAAGGTCCGAAccct                                                                                                                        <  2:3804627/150‑3 (MQ=255)
                           cTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCCCTCCCACACGCATTGCACAGGGTTCGGACCTtgt                                                                                                                      >  2:4870497/1‑150 (MQ=255)
                           cTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCAATGCACAGGGTTCGGACCgtgt                                                                                                                      >  2:1467823/1‑150 (MQ=255)
                            tGGAAGCGTATATCCAGACCGCGACCACGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTATACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                     <  2:690970/150‑1 (MQ=255)
                            tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCTGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                     <  2:82798/150‑1 (MQ=255)
                            tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGGGGGCTTGGCCGCTGCCACAGGGATTGGGCGGGGTTCGGACCTtgtg                                                                                                                     >  2:2589003/1‑150 (MQ=255)
                               aaGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAg                                                                                                                  <  1:4030451/150‑1 (MQ=255)
                                  cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                               <  1:3271386/150‑5 (MQ=25)
                                  cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTc                                                                                                               <  1:4115180/150‑1 (MQ=255)
                                  cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCACTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                               <  2:652466/150‑5 (MQ=37)
                                  cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAACCCTGTGCaaggc                                                                                                               <  2:2868263/150‑5 (MQ=25)
                                    tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgcct                                                                                                             <  2:3711486/150‑7 (MQ=21)
                                    tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgcct                                                                                                             <  2:3711628/150‑7 (MQ=21)
                                    tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCt                                                                                                             <  2:3207984/150‑1 (MQ=255)
                                    tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTCGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAAGGTCCGAACCCTGTGCaatgcct                                                                                                             <  2:3668115/150‑7 (MQ=37)
                                      tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                           <  2:1644682/150‑1 (MQ=255)
                                      tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                           <  2:1961903/150‑1 (MQ=255)
                                       atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                          >  2:4557781/1‑150 (MQ=255)
                                       atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                          <  2:4740604/150‑1 (MQ=255)
                                       atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCcgtg                                                                                                          >  2:2067952/1‑150 (MQ=255)
                                         ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtggg                                                                                                        <  1:1165064/150‑3 (MQ=18)
                                          cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGa                                                                                                       <  2:4105305/150‑1 (MQ=255)
                                          cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGa                                                                                                       <  1:981205/150‑1 (MQ=255)
                                          cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtg                                                                                                          <  2:2477993/147‑1 (MQ=21)
                                          cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtg                                                                                                          >  1:2477993/1‑147 (MQ=21)
                                           aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAGTCGCGGGCTTGCCCCCGCCCACAGGCCGTGGACCCGGGTCGGGCCGTGTGCCGGGCCGGGGCgtg                                                                                                      >  2:1013293/1‑148 (MQ=255)
                                             aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCACTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCg                                                                                                    <  1:2680909/150‑1 (MQ=25)
                                              ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                   <  1:1425261/150‑1 (MQ=18)
                                              ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                   <  1:368006/150‑1 (MQ=255)
                                              ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                   <  2:2071133/150‑1 (MQ=255)
                                              ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                   <  1:368131/150‑1 (MQ=255)
                                              ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                   <  2:3549583/150‑1 (MQ=255)
                                                gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCCCTCCCACAGGCATTGCACAGGGTTAGGACCTTGTGCAGGGCCAGTGCGTGCGGGc                                                                                                 >  1:5116997/1‑150 (MQ=255)
                                                   aCCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCCCGCCCCCAGGGAGTGCGCCGGGGTCGGACCGTGTGCAGGTCCGGTGCGGGCGGGCCAg                                                                                              >  2:1381560/1‑150 (MQ=255)
                                                    ccAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGc                                                                                             <  2:4507190/150‑1 (MQ=255)
                                                             tgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCAGTGCGAGCGGGCAAGCCCGCGAATg                                                                                    >  2:3349119/1‑150 (MQ=255)
                                                                  gATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACACGCAATGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAAAGGGccc                                                                               >  2:2016099/1‑150 (MQ=255)
                                                                  gATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAGTCGCGGGCTTGCCCGCTGCCACAGGCCGTGCACAGGGGTCGGGCCGTGTGCAGGGCCTGTGGGTGCGGGCAAGCCCGCGCCTGGGccc                                                                               >  2:3106528/1‑150 (MQ=255)
                                                                       cacCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGGTCGGACCGTGTGGAGGTGCTGTGCGAGCGGGGAAGCCCGCGCCAGGGCCCCGcac                                                                          >  1:2527583/1‑150 (MQ=38)
                                                                       cacCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTAGCGGGCTTGCCCGCTCCCACACGCATTGCACCGGGTTCGGACCTTGTGCAGGTCCAGTGCGAGCGGGCCAGCCCGCGAATGGGCCCCGcac                                                                          >  2:4325644/1‑150 (MQ=39)
                                                                        acCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGcaca                                                                         >  2:2085870/1‑150 (MQ=255)
                                                                              ggCGTGGTGCACGGCACACCGACGTAACCCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACTGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTc                                                                   <  2:1123662/150‑1 (MQ=255)
                                                                                             acacCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGGTGCCACAGGCATTGCACAGGGTTCGGTCCTTGTGGAGGTGCTGTGCGAGCGGGGAAGGCCGGGTATTGGGCCCGCACAATTTTCCCGGAGTGTAAATTa                                                    >  2:2490975/1‑150 (MQ=255)
                                                                                                ccGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAGTGCACAGGGGTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATcc                                                 <  2:3702044/150‑1 (MQ=255)
                                                                                                  gAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCACGCCCGTGAATGGGCCCCGCCCAATTTTCCCGGACTGAACCTTATcccc                                               >  2:3144442/1‑150 (MQ=255)
                                                                                                      tAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcac                                                                                                                                      >  2:4556195/1‑59 (MQ=255)
                                                                                                      tAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcac                                                                                                                                      <  1:4556195/59‑1 (MQ=255)
                                                                                                          ccGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACAGTATCCCCGCCCCAAg                                       >  2:3817847/1‑150 (MQ=255)
                                                                                                            tgggCCGCGTTGCGGCCCATTCGCGGGCGTGCCCCCTCCCACAGGCCATGCACAGGGTTCGGGCCGTTTTCAGGTGCTTTGCGAGCGGGCCACCCCGCGCAAGGGCCCCCCCCAAATTTACCCGACAGAACATTATACCCCCCCCAAggg                                     >  2:4331737/2‑149 (MQ=255)
                                                                                                                ccGCGTTGCGGCCCATTCGCGGGCTTGCCCGCCCCCACAGGCATTGCACCGGGTTCGGACCTTGTGCAGGTCCTGTGCGCGCGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGt                                 <  1:834362/150‑1 (MQ=255)
                                                                                                                         ggCCCAGCCGCGGGCTGGCCCGCCCCCCCCGGCTTTGCCCGGGGTGGGGCCGTGGTGCCGGTCGTGGGGGCGGGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGc                        <  1:3656066/150‑1 (MQ=255)
                                                                                                                          gCCCATTCGCGGGCTTGCCCGCTCGCACAGGCACGGCACAAGGTGCGCACCCAGTGCAATGCCGGTGGGGGCGGGCAAGCCCGCGAAAGGGCCCCCCACAAATTTACCGGACTGAAAAATAACCCCGCCCCACGGTCAGTTATTCAGGCa                       >  2:891565/1‑150 (MQ=9)
                                                                                                                             cATTGGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGCGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATca                    <  2:897741/150‑1 (MQ=255)
                                                                                                                                  gcgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCATGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTACCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCacacac               >  1:2667879/1‑150 (MQ=255)
                                                                                                                                     ggCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTg            <  2:4527805/150‑1 (MQ=255)
                                                                                                                                     ggCTTGCCCGCACCCACTGGCTTTGCACGGGGTTGGGCCGTTGTGCCGGGCCTGTGCGCGCGCGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTg            <  1:3004548/150‑1 (MQ=255)
                                                                                                                                        ttGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtg         >  2:3498281/1‑150 (MQ=255)
                                                                                                                                           cccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc      <  2:4035398/150‑1 (MQ=255)
                                                                                                                                            ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                      >  2:4722116/1‑133 (MQ=17)
                                                                                                                                            ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                      <  1:4722116/133‑1 (MQ=17)
                                                                                                                                            ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg     >  2:458275/1‑150 (MQ=18)
                                                                                                                                             cGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTCTTCAGGCATCACACACCTGGTGTTCGt    >  1:3545844/1‑150 (MQ=18)
                                                                                                                                             cGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt    >  1:827897/1‑150 (MQ=18)
                                                                                                                                               cTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCACTGCCGGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTg  >  2:4687063/1‑150 (MQ=18)
                                                                                                                                               cTCCCCCAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTg  <  2:446245/150‑1 (MQ=255)
                                                                                                                                               cTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTg  >  2:4489287/1‑150 (MQ=255)
                                                                                                                                               cTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTg  >  2:1507267/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTG  >  NC_002947/2065004‑2065296

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: