Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,148,876 T→A 6.8% intergenic (+36/‑87) PP_3647 → / → PP_3648 oxidoreductase/carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,148,8760TA 6.8% 690.4 / 29.8 220intergenic (+36/‑87)PP_3647/PP_3648oxidoreductase/carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base T (105/97);  new base A (7/8);  total (112/108)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 226 total aligned reads displayed.

GGTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTAC  >  NC_002947/4148727‑4149025
                                                                                                                                                     |                                                                                                                                                     
      cACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgc                                                                                                                                                 >  1:1173169/1‑150 (MQ=255)
      cACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgc                                                                                                                                                 >  1:526006/1‑150 (MQ=255)
        cgcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGa                                                                                                                                               >  1:5283113/1‑150 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              <  1:3522895/150‑1 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              <  2:789129/150‑1 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              <  2:4423078/150‑1 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGGTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              <  2:6534749/150‑1 (MQ=255)
          cgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAg                                                                                                                                             <  1:2119245/150‑1 (MQ=255)
           gATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGa                                                                                                                                            >  2:583000/1‑150 (MQ=255)
           gATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGa                                                                                                                                            <  1:7085560/150‑1 (MQ=255)
            aTCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGCGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGAt                                                                                                                                           <  1:3224949/150‑1 (MQ=255)
            aTCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGAt                                                                                                                                           <  2:1972043/150‑1 (MQ=255)
                  gTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGAGGGTGTACCACAAGGGCACGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTGGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCt                                                                                                                                     <  2:2523543/150‑1 (MQ=255)
                  gTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCAGATGACCCCGATGTTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATAAATCTTGAGTCGAGGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCt                                                                                                                                     <  2:1409548/150‑1 (MQ=255)
                    cGGCGCGGTGAAGAACCTGGTGGTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGt                                                                                                                                   <  2:5400188/150‑1 (MQ=255)
                     ggCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGtt                                                                                                                                  >  2:3358558/1‑150 (MQ=255)
                         cgGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGTTGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTTATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCa                                                                                                                              <  1:639960/150‑1 (MQ=255)
                         cgGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCa                                                                                                                              <  2:3392226/150‑1 (MQ=255)
                         cgGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCa                                                                                                                              <  1:6801045/150‑1 (MQ=255)
                              aagaaCCTGGTGTAGCCTGAGCAGCAGATGACCCAGATCGTGTACCACAAGTGAAAGCTCGGCGCCCTGACGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCACGCGAAGATGCCCCTGTTTTCAGGGGt                                                                                                                         >  2:2027222/1‑150 (MQ=255)
                               agaaCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTc                                                                                                                        <  1:4520229/150‑1 (MQ=255)
                                  aCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTcgcg                                                                                                                     >  2:6168140/1‑150 (MQ=255)
                                   ccTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCgcgc                                                                                                                    <  2:1173169/150‑1 (MQ=255)
                                   ccTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCgcgc                                                                                                                    >  1:1077485/1‑150 (MQ=255)
                                   ccTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCgcgc                                                                                                                    <  1:2226319/150‑1 (MQ=255)
                                   ccTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCgcgc                                                                                                                    <  2:4755725/150‑1 (MQ=255)
                                    cTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCt                                                                                                                   <  1:5666476/150‑1 (MQ=255)
                                     tGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTa                                                                                                                  >  2:2990400/1‑150 (MQ=255)
                                       gtgtTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAAt                                                                                                                >  2:522039/1‑150 (MQ=255)
                                        tgtTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCACGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGTTGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATc                                                                                                               <  1:5720909/150‑1 (MQ=255)
                                          ttGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGt                                                                                                             <  2:628233/150‑1 (MQ=255)
                                          ttGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGt                                                                                                             >  1:770533/1‑150 (MQ=255)
                                           tGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtc                                                                                                            <  1:4542748/150‑1 (MQ=255)
                                            gCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtca                                                                                                           >  2:6003592/1‑149 (MQ=255)
                                             ccTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtctc                                                                                                          >  2:1636088/1‑150 (MQ=255)
                                             ccTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtctc                                                                                                          >  1:257694/1‑150 (MQ=255)
                                                gAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCAcc                                                                                                       <  2:2730851/150‑1 (MQ=255)
                                                  gcagcCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCt                                                                                                     >  2:3868494/1‑150 (MQ=255)
                                                      ccGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGt                                                                                                 >  1:3811702/1‑150 (MQ=255)
                                                          tGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCg                                                                                             <  2:2642667/150‑1 (MQ=255)
                                                          tGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTATCACCCACGGTTGCg                                                                                             >  1:2697155/1‑150 (MQ=255)
                                                           gACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTGGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCgg                                                                                            <  1:5016481/150‑1 (MQ=255)
                                                              cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                         >  2:4854791/1‑150 (MQ=255)
                                                              cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                         >  1:7155979/1‑150 (MQ=255)
                                                              cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                         >  1:1456892/1‑150 (MQ=255)
                                                              cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                         >  1:2903229/1‑150 (MQ=255)
                                                              cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTACCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                         >  2:6384345/1‑150 (MQ=255)
                                                                cGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCt                                                                                       >  1:852770/1‑150 (MQ=255)
                                                                cGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCt                                                                                       >  2:1326918/1‑150 (MQ=255)
                                                                 gATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCtt                                                                                      >  2:2453698/1‑150 (MQ=255)
                                                                       gtACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTAt                                                                                >  1:5729273/1‑150 (MQ=255)
                                                                        tACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATg                                                                               <  2:6908642/150‑1 (MQ=255)
                                                                            acaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCgg                                                                           >  1:1228670/1‑150 (MQ=255)
                                                                            acaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCgg                                                                           >  1:3880274/1‑150 (MQ=255)
                                                                             caAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                          <  2:2230000/150‑1 (MQ=255)
                                                                             caAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                          <  1:2306092/150‑1 (MQ=255)
                                                                                 ggCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTAc                                                                      >  2:608328/1‑150 (MQ=255)
                                                                                  gCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACt                                                                     <  2:7443191/150‑1 (MQ=255)
                                                                                  gCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACt                                                                     >  1:4347488/1‑150 (MQ=255)
                                                                                  gCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACt                                                                     <  1:4802272/150‑1 (MQ=255)
                                                                                      gCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGcac                                                                 <  1:2990400/150‑1 (MQ=255)
                                                                                        tCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACAt                                                               >  1:5089065/1‑150 (MQ=255)
                                                                                        tCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACAt                                                               <  2:5520781/150‑1 (MQ=255)
                                                                                          ggCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGa                                                             <  1:4184477/150‑1 (MQ=255)
                                                                                              cccTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCaa                                                         <  1:7103642/150‑1 (MQ=255)
                                                                                                cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                          >  2:3436681/1‑131 (MQ=255)
                                                                                                cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                          <  1:3436681/131‑1 (MQ=255)
                                                                                                    cgccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaaaa                                                   <  2:5614883/150‑1 (MQ=255)
                                                                                                    cgccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaaaa                                                   <  1:1010720/150‑1 (MQ=255)
                                                                                                     gccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAg                                                  >  1:1548303/1‑150 (MQ=255)
                                                                                                               tGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAgg                                        >  2:2102080/1‑150 (MQ=255)
                                                                                                               tGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAgg                                        >  2:7171069/1‑150 (MQ=255)
                                                                                                                gATACCTCTTGAGTCGTTGTTTGCATCTTCGCGGGCTTGCCCGCGAAGATGCCGCTGTTTTCAGGGGTCGCGATAATCGTATCACCCTCGGTTGAGGCCACTTTTTTATGCCGTATTACTGCTCTTGAGCAATGTAAAGTACGAAAAggg                                       <  2:4990042/150‑1 (MQ=255)
                                                                                                                   aaaTCTTGAGTCGATGTTGGCATCTTCGCGGGCTGGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCt                                    <  1:133102/148‑1 (MQ=255)
                                                                                                                   aCATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCt                                    <  1:1784590/150‑1 (MQ=255)
                                                                                                                      tCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGaa                                 >  2:2872237/1‑150 (MQ=255)
                                                                                                                      tCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGaa                                 >  1:6995482/1‑150 (MQ=255)
                                                                                                                      tCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGaa                                 >  1:1873023/1‑150 (MQ=255)
                                                                                                                       cTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGaaa                                >  1:1771213/1‑150 (MQ=255)
                                                                                                                        ttGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGaaaa                               >  1:5760199/1‑150 (MQ=255)
                                                                                                                          gAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATc                             <  2:852770/150‑1 (MQ=255)
                                                                                                                           aGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATcc                            <  2:106511/150‑1 (MQ=255)
                                                                                                                             tCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGc                          >  2:3623367/1‑150 (MQ=255)
                                                                                                                                aTGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAATGTACGACAAGGGCCTGAAAATCCGCAcc                       >  2:7052532/1‑150 (MQ=255)
                                                                                                                                    tGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAgg                   >  2:3996396/1‑150 (MQ=255)
                                                                                                                                       cATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTAc                >  1:3305164/1‑150 (MQ=255)
                                                                                                                                          cTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTgg             <  1:3987892/150‑1 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTCCGAAAAGGGCCTGAAAATCCGCACCCAGGTCCgggg            >  2:4645914/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg            >  2:1355930/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg            >  1:4274030/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg            >  1:4641541/1‑150 (MQ=255)
                                                                                                                                            tCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGc           >  1:1414145/1‑150 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGa         <  2:1795922/150‑1 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGa         <  1:6648284/150‑1 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAg        <  1:2073834/150‑1 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAg        <  2:5839464/150‑1 (MQ=255)
                                                                                                                                                gggCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAgg       <  1:4385113/150‑1 (MQ=255)
                                                                                                                                                gggCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAgg       <  1:5878852/150‑1 (MQ=255)
                                                                                                                                                   cTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACt    <  2:7125556/150‑1 (MQ=255)
                                                                                                                                                   cTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACt    >  1:1675036/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTAC  >  NC_002947/4148727‑4149025

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