Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,143,926 T→A 11.1% intergenic (+35/‑88) 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,143,9260TA11.1% 177.7 / 14.6 63intergenic (+35/‑88)PP_3647/PP_3648oxidoreductase/carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base A (3/4);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

GGTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGG  >  NC_002947/4143778‑4144066
                                                                                                                                                    |                                                                                                                                            
ggTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCTCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCtt                                                                                                                                             <  1:1922296/150‑1 (MQ=255)
ggTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCtt                                                                                                                                             <  2:77940/150‑1 (MQ=255)
      cACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgc                                                                                                                                       >  1:780452/1‑150 (MQ=255)
         gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGaa                                                                                                                                    <  2:410668/150‑1 (MQ=255)
              cATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGc                                                                                                                               <  1:70426/150‑1 (MQ=255)
              cATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGc                                                                                                                               <  2:807599/150‑1 (MQ=255)
                ttGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAATCGATCTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGccc                                                                                                                             >  1:1876013/1‑150 (MQ=255)
                         cgGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCa                                                                                                                    <  2:1889087/150‑1 (MQ=255)
                          ggTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTACTCGTGCTTGCCCGCGAAGATGCCCCTGTTTTCAg                                                                                                                   <  1:2301624/150‑1 (MQ=255)
                           gTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAgg                                                                                                                  >  2:362434/1‑150 (MQ=255)
                                   ccTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGa                                                                                                                                  <  1:769279/126‑1 (MQ=255)
                                   ccTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGa                                                                                                                                  >  2:769279/1‑126 (MQ=255)
                                    cTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCt                                                                                                         >  2:859/1‑150 (MQ=255)
                                        tgtTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATc                                                                                                     >  2:777333/1‑150 (MQ=255)
                                            gCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtct                                                                                                 >  1:73640/1‑150 (MQ=255)
                                                   cagcCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTc                                                                                          <  2:73640/150‑1 (MQ=255)
                                                        gATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCGCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTg                                                                                     >  1:1739311/1‑150 (MQ=255)
                                                         aTGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGc                                                                                    <  1:777333/150‑1 (MQ=255)
                                                              cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                               >  2:1182965/1‑150 (MQ=255)
                                                              cacgttGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAgg                                                                                                                  >  2:1013340/6‑115 (MQ=255)
                                                              cacgttGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAgg                                                                                                                  <  1:1013340/110‑1 (MQ=255)
                                                                 gATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGc                                                                                                                               <  2:1481827/99‑1 (MQ=255)
                                                                 gATGGTGTACCACAAGGGCAAGCTAGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGc                                                                                                                               >  1:1481827/1‑99 (MQ=255)
                                                                 gATGGTGTACCACAAGGGCAAGCTAGGCGCCCTGACGTAACATGTATCATACAGGTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCGAATCGTCTCACCCTCGGTTGCTGCCCCtt                                                                            >  1:1572613/1‑150 (MQ=255)
                                                                  aTGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttt                                                                           <  1:1467290/150‑1 (MQ=255)
                                                                        tACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCAGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATg                                                                     >  2:597204/1‑150 (MQ=255)
                                                                           cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTcgcg                                                                                                           >  1:951939/1‑109 (MQ=255)
                                                                           cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTcgcg                                                                                                           <  2:951939/109‑1 (MQ=255)
                                                                           cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCg                                                                  >  1:230192/1‑150 (MQ=255)
                                                                            acaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATTCCCCTGTTTTCAGGGGTCGCGCTAAt                                                                                                      >  1:1032058/1‑113 (MQ=255)
                                                                            acaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATTCCCCTGTTTTCAGGGGTCGCGCTAAt                                                                                                      <  2:1032058/113‑1 (MQ=255)
                                                                            acaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCgg                                                                 >  2:58062/1‑150 (MQ=255)
                                                                             caAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                <  1:58062/150‑1 (MQ=255)
                                                                                   cAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGTAGTACTg                                                          <  2:1739311/150‑1 (MQ=255)
                                                                                    aaGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGc                                                         >  1:172923/1‑150 (MQ=255)
                                                                                      gCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGcac                                                       >  2:1325570/1‑150 (MQ=255)
                                                                                        tCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACAt                                                     <  2:1572613/150‑1 (MQ=255)
                                                                                               ccTGCCGCCACTTGGCTGATACATCTTCGGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAAt                                              <  1:859/150‑1 (MQ=255)
                                                                                                cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtct                                                                                                 >  1:2142093/1‑98 (MQ=255)
                                                                                                cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtct                                                                                                 <  2:2142093/98‑1 (MQ=255)
                                                                                                cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATg                                             <  1:1182965/150‑1 (MQ=255)
                                                                                                      ccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGt                                       <  2:2337981/150‑1 (MQ=255)
                                                                                                        aCTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACAt                                                     >  2:900316/1‑134 (MQ=255)
                                                                                                        aCTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACAt                                                     <  1:900316/134‑1 (MQ=255)
                                                                                                               tGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGc                                                                    <  2:1401117/112‑1 (MQ=255)
                                                                                                               tGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGc                                                                    >  1:1401117/1‑112 (MQ=255)
                                                                                                                 aTACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGc                            >  1:251280/1‑150 (MQ=255)
                                                                                                                     aTCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGtt                                                                                      <  1:422419/88‑1 (MQ=255)
                                                                                                                     aTCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGtt                                                                                      >  2:422419/1‑88 (MQ=255)
                                                                                                                        ttGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGaaaa                     >  1:428386/1‑150 (MQ=255)
                                                                                                                          gAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATc                   <  2:251280/150‑1 (MQ=255)
                                                                                                                            gTCGATGTTGGCATGCTCGCGGGCTTGCCAGCGAAGGTGCCACTGGTTTCCGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTGTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCg                 <  2:1876013/150‑1 (MQ=255)
                                                                                                                            gTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGCTTTCAGGGGCCGCGCTAATCGTCCCACCCTCGGCTGCGGCCCCTTTTTTATCCCCGAGAACTGCACATGAGCCCTGACGAGTACGAAAAGAACCTGAAAATCCg                 >  2:1572700/1‑150 (MQ=255)
                                                                                                                             tCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGc                >  1:1298713/1‑150 (MQ=255)
                                                                                                                               gATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGa                                                   <  1:1990771/113‑1 (MQ=255)
                                                                                                                               gATATTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGa                                                   >  2:1990771/1‑113 (MQ=255)
                                                                                                                                 tGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGATTTCAGGGGTAGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCAccc            >  1:2222714/1‑150 (MQ=255)
                                                                                                                                    gggCATCTTCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtc                                                                                                  >  2:2087504/2‑61 (MQ=21)
                                                                                                                                    gggCATCTTCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtc                                                                                                  <  1:2087504/60‑1 (MQ=21)
                                                                                                                                          ccgcgcgGGCTTGTCCGCTAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTgg   <  2:1298713/147‑1 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg  >  1:2113783/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg  >  2:1217850/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg  <  1:1325570/150‑1 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                         >  1:1837689/1‑79 (MQ=38)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                         <  2:1837689/79‑1 (MQ=38)
                                                                                                                                                    |                                                                                                                                            
GGTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGG  >  NC_002947/4143778‑4144066

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