Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,846,028 T→C 10.3% intergenic (+44/+47) yfhL → / ← mutM 4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,846,0280TC10.3% 271.2 / 16.0 97intergenic (+44/+47)yfhL/mutM4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (45/42);  new base C (5/5);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTA  >  NC_002947/5845883‑5846176
                                                                                                                                                 |                                                                                                                                                    
cagcagGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGTAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTc                                                                                                                                                  <  1:1203957/150‑1 (MQ=255)
cagcagGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTc                                                                                                                                                  >  1:557015/1‑150 (MQ=255)
    agGCCTGCCCGGTGGACTGCATCCCGCTGGATGAAGCGTAGCCGGAAACCTAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCATGCACCGTCCCTGTCGGAGCAGGCATGCCCGCTCTcac                                                                                                                                              <  2:1954596/150‑1 (MQ=255)
        cTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGtt                                                                                                                                          <  2:1404166/150‑1 (MQ=255)
          gCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTa                                                                                                                                        <  2:1631048/150‑1 (MQ=255)
           cccGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGACACACTTCCGATGCCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAg                                                                                                                                       <  2:64802/150‑1 (MQ=255)
           cccGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAg                                                                                                                                       <  1:2367166/150‑1 (MQ=255)
           cccGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAg                                                                                                                                       <  1:13940/150‑1 (MQ=255)
                 acaCTGCATCCCGCTGGATGTAGCGCATTCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTTTCAAGGGCTGATCATTTCTCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCg                                                                                                                                 <  2:1069903/148‑1 (MQ=255)
                    cTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCtt                                                                                                                              <  1:57821/150‑1 (MQ=255)
                    cTGCATCACGCTGGATGAAGCGCAGACGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCtt                                                                                                                              <  2:557015/150‑1 (MQ=255)
                     tGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCcac                                                                                                                                              >  1:2056434/1‑133 (MQ=255)
                     tGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCcac                                                                                                                                              <  2:2056434/133‑1 (MQ=255)
                     tGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTc                                                                                                                             >  1:849560/1‑150 (MQ=255)
                       cATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCaa                                                                                                                           >  2:1688493/1‑150 (MQ=255)
                        aTCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAAc                                                                                                                          >  2:1150832/1‑150 (MQ=255)
                             gCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGtt                                                                                                                     >  1:1848451/1‑150 (MQ=255)
                                 gATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCa                                                                                                                 >  2:2209548/1‑150 (MQ=255)
                                  aTGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTCTGAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAg                                                                                                                <  1:25536/150‑1 (MQ=255)
                                       gcgcATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAAc                                                                                                           >  1:133870/1‑150 (MQ=255)
                                       gcgcATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAAc                                                                                                           >  1:2083994/1‑150 (MQ=255)
                                         gcATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGt                                                                                                         <  2:418521/150‑1 (MQ=255)
                                           aTCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGttt                                                                                                       >  2:1324795/1‑150 (MQ=255)
                                           aTCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGttt                                                                                                       <  2:1317922/150‑1 (MQ=255)
                                           aTCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGttt                                                                                                       >  1:139635/1‑150 (MQ=255)
                                               ggAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATc                                                                                                   >  1:324961/1‑150 (MQ=255)
                                                  aaCCGAAGAACAGCTGATGTCCAAGTATCGCCGGATCACTTGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCTCTGTGATAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACATGTTCGCAGCGAACGTTTCATCGtt                                                                                                <  2:467052/150‑1 (MQ=255)
                                                  aaCCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGtt                                                                                                >  2:154139/1‑150 (MQ=255)
                                                      gaagaaCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCa                                                                                            <  1:1168043/150‑1 (MQ=255)
                                                           aCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAgg                                                                                                                       >  2:2161452/1‑118 (MQ=255)
                                                           aCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAgg                                                                                                                       <  1:2161452/118‑1 (MQ=255)
                                                            cAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTcaca                                                                                                                                             <  1:182100/95‑1 (MQ=38)
                                                            cAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTcaca                                                                                                                                             >  2:182100/1‑95 (MQ=38)
                                                             aGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                              >  1:1266721/1‑141 (MQ=255)
                                                             aGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                              <  2:1266721/141‑1 (MQ=255)
                                                               cTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCa                                                                                   <  1:122690/150‑1 (MQ=255)
                                                               cTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCa                                                                                   <  1:2362571/150‑1 (MQ=255)
                                                                tGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTc                                                                                                                             <  1:2058819/107‑1 (MQ=255)
                                                                tGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTc                                                                                                                             >  2:2058819/1‑107 (MQ=255)
                                                                  aTGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTa                                                                                <  2:139635/150‑1 (MQ=255)
                                                                      ccAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGtt                                                                                                <  2:1423204/130‑1 (MQ=255)
                                                                      ccAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGtt                                                                                                >  1:1423204/1‑130 (MQ=255)
                                                                      ccAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACg                                                                            >  2:1818799/1‑150 (MQ=255)
                                                                       cAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTa                                                                                                                                        <  1:666768/89‑1 (MQ=255)
                                                                       cAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTa                                                                                                                                        >  2:666768/1‑89 (MQ=255)
                                                                       cAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAg                                                                                                                                       >  2:689634/1‑90 (MQ=255)
                                                                       cAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAg                                                                                                                                       <  1:689634/90‑1 (MQ=255)
                                                                       cAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGc                                                                           <  2:1031136/150‑1 (MQ=255)
                                                                        aaGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCt                                                                          >  1:926014/1‑150 (MQ=255)
                                                                         aGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTc                                                                         <  1:1688493/150‑1 (MQ=255)
                                                                              cGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCct                                                                    <  1:145578/150‑1 (MQ=255)
                                                                               gCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCa                                                                                   >  1:172071/1‑134 (MQ=255)
                                                                               gCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCa                                                                                   <  2:172071/134‑1 (MQ=255)
                                                                               gCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCctc                                                                   >  1:147001/1‑150 (MQ=255)
                                                                                ccGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCctct                                                                  >  2:1865914/1‑150 (MQ=255)
                                                                                 cGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTg                                                                 <  2:1879898/150‑1 (MQ=255)
                                                                                   gATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCggg                                                                                     >  2:1049060/1‑128 (MQ=255)
                                                                                   gATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCggg                                                                                     <  1:1049060/128‑1 (MQ=255)
                                                                                   gATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGcc                                                               >  1:836368/1‑150 (MQ=255)
                                                                                     tCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCACCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGcccc                                                             >  1:251019/1‑150 (MQ=255)
                                                                                      cACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGATCGGGCATGCCTGTTCCCACAGTTTAGGTTTCGCCTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCa                                                            <  1:1367219/150‑1 (MQ=255)
                                                                                                ggCTGATCACTGCCATTCCTTGCCTCCGCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGTTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTc                                                  <  2:1159834/150‑1 (MQ=255)
                                                                                                      tCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCACCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCaa                                            <  2:251019/150‑1 (MQ=255)
                                                                                                        aTTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTCCCCCAGCTTCACTTCGCGCAACt                                          >  1:1456421/1‑150 (MQ=255)
                                                                                                         ttGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTc                                         >  1:1616212/1‑150 (MQ=255)
                                                                                                            cccTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGt                                      >  1:451909/1‑150 (MQ=255)
                                                                                                                  ccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTa                                                                                <  1:1106265/102‑1 (MQ=255)
                                                                                                                  ccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTa                                                                                >  2:1106265/1‑102 (MQ=255)
                                                                                                                   cTGCACGGTCCCTGTGGGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCAt                               >  1:628592/1‑150 (MQ=255)
                                                                                                                   cTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTc                                                                                                      >  1:1361432/1‑79 (MQ=255)
                                                                                                                   cTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTc                                                                                                      <  2:1361432/79‑1 (MQ=255)
                                                                                                                   cTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCAt                               >  1:2285146/1‑150 (MQ=255)
                                                                                                                        cGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCtt                          <  2:926014/150‑1 (MQ=255)
                                                                                                                          taaacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGc                        >  1:1393328/5‑150 (MQ=255)
                                                                                                                          gTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTg                                                                 >  1:1586132/1‑109 (MQ=255)
                                                                                                                          gTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTg                                                                 <  2:1586132/109‑1 (MQ=255)
                                                                                                                           tCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCa                       <  1:1865914/150‑1 (MQ=255)
                                                                                                                           tCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCa                       <  2:1848451/150‑1 (MQ=255)
                                                                                                                           aaacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCa                       <  2:1393328/147‑1 (MQ=255)
                                                                                                                            aacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGc                                                          <  1:401173/112‑1 (MQ=38)
                                                                                                                            aacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGc                                                          >  2:401173/3‑114 (MQ=38)
                                                                                                                              cTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgc                                                <  2:21146/122‑1 (MQ=255)
                                                                                                                              cTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgc                                                >  1:21146/1‑122 (MQ=255)
                                                                                                                              cTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACgg                    >  2:700838/1‑150 (MQ=255)
                                                                                                                               tgtgGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACt                                                    <  1:1672448/117‑1 (MQ=39)
                                                                                                                               tgtgGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACt                                                    >  2:1672448/1‑117 (MQ=39)
                                                                                                                               tgtgGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGc                   >  1:1259800/1‑150 (MQ=255)
                                                                                                                                   agagCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGcc               >  1:1569106/1‑150 (MQ=255)
                                                                                                                                     agCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCg             >  1:314250/1‑150 (MQ=255)
                                                                                                                                         ggCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCAc         >  1:1950587/1‑150 (MQ=255)
                                                                                                                                          gCATGCCTGCTCCCACAGTTTAAGGTTCGCTTCAACAGGTTCGCAGCGAACGTTTCACCGTTGGCAGCGCGGGCCGTACACGCTCCCCCTCTGCCCCAGCTTCACTTCGCGCAACCCACTGCCGCATACCTTGCACGGCAGCCCGCCCCg        >  2:799645/1‑150 (MQ=255)
                                                                                                                                            aTGCCTGCTCCCACAGTTTATGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGcc      >  2:404417/1‑150 (MQ=255)
                                                                                                                                               ccTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGt   >  1:745936/1‑150 (MQ=255)
                                                                                                                                                cTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTc                                                                         <  1:878355/79‑1 (MQ=255)
                                                                                                                                                cTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTc                                                                         >  2:878355/1‑79 (MQ=255)
                                                                                                                                                cTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTa  <  2:369159/150‑1 (MQ=255)
                                                                                                                                                cTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTa  >  1:1649006/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTA  >  NC_002947/5845883‑5846176

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