Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,120 A→C 12.7% intergenic (‑130/+24) PP_4560 ← / ← csbD ribonuclease BN/stress response protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,1200AC12.7% 156.9 / 16.2 71intergenic (‑130/+24)PP_4560/csbDribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base A (34/28);  new base C (2/7);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 8.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

AGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGT  >  NC_002947/5181974‑5182264
                                                                                                                                                  |                                                                                                                                                  
aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                                 <  1:2107480/150‑1 (MQ=255)
  gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTTCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaaga                                                                                                                                               <  1:344199/150‑1 (MQ=255)
  gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaaga                                                                                                                                               <  2:1893474/150‑1 (MQ=255)
   tCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGAAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                              <  2:1639136/150‑1 (MQ=255)
    cGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAgg                                                                                                                                             <  2:2078491/150‑1 (MQ=255)
       ggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGGCCg                                                                                                                                          <  2:553558/150‑1 (MQ=39)
           aaTCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCAc                                                                                                                                      <  1:1063615/150‑1 (MQ=255)
              cATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATATGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGc                                                                                                                                   <  1:1870475/150‑1 (MQ=255)
              cATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGc                                                                                                                                   >  2:148630/1‑150 (MQ=255)
                          tttCGCCGCTTGATCGTTACACAGGCAAGTCCAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAgg                                                                                                                       >  2:1612670/1‑150 (MQ=255)
                                      aTCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGc                                                                                                           <  2:1099087/150‑1 (MQ=255)
                                       tCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAGGTGCCGGCCTCTTCTCGGGTGAACCCGcgaa                                                                                                                                                 >  1:371942/1‑111 (MQ=38)
                                       tCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAGGTGCCGGCCTCTTCTCGGGTGAACCCGcgaa                                                                                                                                                 <  2:371942/111‑1 (MQ=38)
                                         gTTACACAGGCAAGTCTAATAAACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGCATAAGTGCCGGCCTCTTCGCGGGGTTACCCGAGAAGAGGCGGGGCCCGCCTCAAGTCAAGGCTTGTTCACCGCTTc                                                                                                        >  2:568380/1‑150 (MQ=255)
                                          ttACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGc                                                                                                                      <  1:1439067/135‑1 (MQ=255)
                                          ttACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGc                                                                                                                      >  2:1439067/1‑135 (MQ=255)
                                          ttACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATcc                                                                                                       >  2:646676/1‑150 (MQ=255)
                                              acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTca                                                                                                   <  1:1443287/150‑1 (MQ=255)
                                                aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTcaca                                                                                                 >  1:2271483/1‑150 (MQ=255)
                                                aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTcaca                                                                                                 >  1:2271904/1‑150 (MQ=255)
                                                  gCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAGGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                              >  2:1380386/1‑103 (MQ=37)
                                                  gCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAATTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                              <  1:1380386/103‑1 (MQ=37)
                                                  gCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACAtt                                                                                               >  1:2223312/1‑150 (MQ=255)
                                                       tCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCt                                                                                          >  1:1836884/1‑150 (MQ=255)
                                                        cTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCGGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCtt                                                                                         <  2:1687991/150‑1 (MQ=255)
                                                           aTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGAc                                                                                      <  2:1151574/150‑1 (MQ=255)
                                                            tAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGAcc                                                                                     >  2:599262/1‑150 (MQ=255)
                                                            tAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGAcc                                                                                     >  2:1811760/1‑150 (MQ=255)
                                                              gACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaa                                                                                                                                                 <  1:1845888/88‑1 (MQ=37)
                                                                 cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGc                                                                                                                      >  2:118048/1‑112 (MQ=255)
                                                                 cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGc                                                                                                                      <  1:118048/112‑1 (MQ=255)
                                                                      ccGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGcc                                                                                                                                           >  1:1567662/1‑86 (MQ=255)
                                                                      ccGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGcc                                                                                                                                           <  2:1567662/86‑1 (MQ=255)
                                                                      ccGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcct                                                                           <  2:1836884/150‑1 (MQ=255)
                                                                       cGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGCCCTCTTCGCGGGATCACCCGAGAAGAGGCCGCCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACCTTTCCCTTGACCTGCTCcgcctc                                                                          >  2:763742/1‑150 (MQ=255)
                                                                             gtgCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcctc                                                                          >  2:1017481/1‑144 (MQ=255)
                                                                             gtgCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcctc                                                                          <  1:1017481/144‑1 (MQ=255)
                                                                             gtgCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCttt                                                                    <  1:1612670/150‑1 (MQ=255)
                                                                              tgctgcTCCCGGCCGTGGTTGCAGGATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCAGAGCAGGGGACGGCAACGCCTAATGTGGAGGATTCTTCACTTTATCGTTCAGAGTGCCCTGGACCTGCTGCGCCTCGCCTTTg                                                                   <  2:1206630/150‑1 (MQ=255)
                                                                              tgctgcTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTg                                                                   <  2:2271483/150‑1 (MQ=255)
                                                                              tgctgcTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTg                                                                   <  2:2271904/150‑1 (MQ=255)
                                                                               gctgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGa                                                                  <  2:2438833/150‑1 (MQ=255)
                                                                               gctgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGa                                                                  <  1:1782722/150‑1 (MQ=255)
                                                                                    tCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaagtcaag                                                                                                                        >  1:773330/1‑91 (MQ=255)
                                                                                    tCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaagtcaag                                                                                                                        <  2:773330/91‑1 (MQ=255)
                                                                                    tCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTcc                                                           >  1:1549986/1‑150 (MQ=255)
                                                                                                  gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTc                                             >  2:1146945/1‑150 (MQ=255)
                                                                                                     gaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGc                                          >  2:1339240/1‑150 (MQ=255)
                                                                                                            cTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCACGCACCGCCTCAAGTCAAGGCTTCTTCACCGCAGCCTTCACATTTCCCTTGACCGGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCcgccgc                                   >  1:958542/1‑150 (MQ=255)
                                                                                                                 gAATACGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGCCCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                              <  1:568380/150‑1 (MQ=255)
                                                                                                                        tGCCGTCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                       <  2:958542/150‑1 (MQ=255)
                                                                                                                        tGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                       <  1:1146945/150‑1 (MQ=255)
                                                                                                                        tGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                       >  2:1556768/1‑150 (MQ=255)
                                                                                                                        tGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                       >  1:14613/1‑150 (MQ=255)
                                                                                                                           cGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTc                    >  1:548846/1‑150 (MQ=255)
                                                                                                                            ggCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGAGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGc                                   <  2:1330134/134‑1 (MQ=38)
                                                                                                                            ggCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGAGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGc                                   >  1:1330134/1‑134 (MQ=38)
                                                                                                                             gCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                              >  1:197753/1‑138 (MQ=255)
                                                                                                                             gCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                              <  2:197753/138‑1 (MQ=255)
                                                                                                                               ctctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTcc                                                           <  1:1281882/107‑1 (MQ=255)
                                                                                                                               ctctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCCCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTcc                                                           >  2:1281882/1‑107 (MQ=255)
                                                                                                                                tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCGTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCtt               >  1:2043966/1‑150 (MQ=255)
                                                                                                                                tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt                         >  1:927277/1‑140 (MQ=255)
                                                                                                                                tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt                         <  2:927277/140‑1 (MQ=255)
                                                                                                                                tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCtt               <  1:599464/150‑1 (MQ=255)
                                                                                                                                tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCtt               >  1:1211522/1‑150 (MQ=255)
                                                                                                                                  ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc             >  1:200593/1‑150 (MQ=255)
                                                                                                                                   tCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGcc            <  1:599262/150‑1 (MQ=255)
                                                                                                                                   tCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGcc            <  1:1811760/150‑1 (MQ=255)
                                                                                                                                     gcgGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCACAGCCTTTCGTTGTCCGTCACCTTGCCaa          >  2:1101704/1‑150 (MQ=255)
                                                                                                                                           gaaCCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGACACTTACTTTGGACTGAACAAGGGCGTT‑‑CTTGTACCCGCACCTGTTCGTTGTCGGTCACCTTGCCAACACCCt    <  1:44468/148‑1 (MQ=255)
                                                                                                                                             aCCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGt  >  1:28820/1‑150 (MQ=255)
                                                                                                                                             aCCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGt  <  1:1716048/150‑1 (MQ=255)
                                                                                                                                             aCCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGt  >  1:1707875/1‑150 (MQ=255)
                                                                                                                                             aCCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTCCTGCGCCTCGCACCTCT‑‑ACTCACCTGTCTTTCCGTCGTCTCGCAGCTTTTCGTTGTCGACCACCTTGCCAGCACGCttg  >  1:2112513/1‑148 (MQ=255)
                                                                                                                                              cccGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt                         <  1:79280/126‑1 (MQ=255)
                                                                                                                                              cccGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt                         >  2:79280/1‑126 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGA‑‑GCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGT  >  NC_002947/5181974‑5182264

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