Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,845,972 T→C 8.1% 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,845,9720TC 8.1% 211.4 / 6.8 74intergenic (+44/+47)yfhL/mutM4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base C (3/3);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGT  >  NC_002947/5845824‑5846119
                                                                                                                                                    |                                                                                                                                                   
tGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCt                                                                                                                                                     <  3:54444/149‑1 (MQ=255)
    agcagGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTc                                                                                                                                                 >  2:316430/1‑149 (MQ=255)
     gcagGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTcc                                                                                                                                                <  2:82598/149‑1 (MQ=255)
       agGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCca                                                                                                                                              <  6:73461/149‑1 (MQ=255)
        ggTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTcac                                                                                                                                             <  8:349063/149‑1 (MQ=255)
          tCTGCCCGGTCTACTGGATCCCGCTGGATTAAGCGCTTCCGTAACCCGAAGAACAGCTGATGGCCAAGTATCGCCGGATGACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTcaca                                                                                                                                            <  3:291192/148‑1 (MQ=255)
          tCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTCACAg                                                                                                                                           <  8:129288/149‑1 (MQ=255)
              cccGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGtt                                                                                                                                         >  7:2041/1‑147 (MQ=255)
              cccGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTa                                                                                                                                       <  7:236525/149‑1 (MQ=255)
                cGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCCCAGTTTAgg                                                                                                                                     >  8:371033/1‑149 (MQ=255)
                cGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAgg                                                                                                                                     <  6:13207/149‑1 (MQ=255)
                cGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAgg                                                                                                                                     <  6:45319/149‑1 (MQ=255)
                     gACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGtttag                                                                                                                                >  2:127777/1‑147 (MQ=255)
                      aCTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCCGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGc                                                                                                                               <  2:324918/149‑1 (MQ=255)
                        tGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCtt                                                                                                                             >  3:341557/1‑149 (MQ=255)
                          cATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTc                                                                                                                            >  5:277642/1‑148 (MQ=255)
                           aTCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCaa                                                                                                                          >  5:147016/1‑149 (MQ=255)
                              ccGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGtt                                                                                                                                         >  2:275490/1‑131 (MQ=255)
                              ccGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTACCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGtt                                                                                                                                         <  1:275490/131‑1 (MQ=255)
                               cGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAgg                                                                                                                      >  7:43034/1‑149 (MQ=255)
                                    gATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTc                                                                                                                   <  6:277642/147‑1 (MQ=255)
                                           cgcATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAAc                                                                                                          >  7:212019/1‑149 (MQ=255)
                                            gcATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACg                                                                                                         >  5:67939/1‑149 (MQ=255)
                                             cATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGt                                                                                                        <  6:264436/149‑1 (MQ=255)
                                             cATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGt                                                                                                        >  4:297924/1‑149 (MQ=255)
                                             cATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGt                                                                                                        <  1:316430/149‑1 (MQ=255)
                                                ccGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGttt                                                                                                      >  3:64826/1‑148 (MQ=255)
                                                  ggAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCa                                                                                                    >  4:114985/1‑148 (MQ=255)
                                                   gAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATc                                                                                                  <  5:81588/149‑1 (MQ=255)
                                                      aCCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGtt                                                                                               <  8:170569/149‑1 (MQ=255)
                                                       ccGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTg                                                                                              <  6:73217/149‑1 (MQ=255)
                                                        cGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTTAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                             <  4:45952/149‑1 (MQ=255)
                                                        cGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTa                                                                                                                                       >  4:217929/1‑107 (MQ=255)
                                                        cGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                             >  1:202178/1‑149 (MQ=255)
                                                         gaagaaCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGc                                                                                            <  8:2041/149‑1 (MQ=255)
                                                         gaagaaCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTGCAACAGGGTCGCAGCGACGGGTTCAGCGTTgg                                                                                             >  2:372246/1‑148 (MQ=255)
                                                              aCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAgcgc                                                                                       <  3:181502/149‑1 (MQ=255)
                                                                  cTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTTCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                             <  5:395200/139‑1 (MQ=255)
                                                                  cTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                             >  6:395200/1‑139 (MQ=255)
                                                                        gCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTac                                                                              <  5:342178/148‑1 (MQ=255)
                                                                         ccAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTACTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTg                                                                                              >  3:140224/1‑131 (MQ=255)
                                                                         ccAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTACTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTg                                                                                              <  4:140224/131‑1 (MQ=255)
                                                                             gTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTc                                                                        >  1:98930/1‑149 (MQ=255)
                                                                              tATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCg                                                                       <  4:385274/149‑1 (MQ=255)
                                                                                tCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCTGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGc                                                                      <  2:98930/148‑1 (MQ=255)
                                                                                   ccGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCctc                                                                  >  6:315258/1‑149 (MQ=255)
                                                                                          aCTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCctct                                                                 >  8:234022/1‑143 (MQ=255)
                                                                                          aCTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCctct                                                                 <  7:234022/143‑1 (MQ=255)
                                                                                           cTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAg                                                          >  2:348057/1‑149 (MQ=255)
                                                                                              gCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCTTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCt                                                        >  5:120327/1‑148 (MQ=255)
                                                                                              gCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCtt                                                       <  8:43034/149‑1 (MQ=255)
                                                                                                aaGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCa                                                     <  6:147016/149‑1 (MQ=255)
                                                                                                        aTCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGTCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgc                                               <  6:67939/147‑1 (MQ=255)
                                                                                                        aTCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgc                                               <  1:348057/147‑1 (MQ=255)
                                                                                                          cATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgcg                                              >  3:191538/1‑146 (MQ=255)
                                                                                                          cATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgcg                                              <  4:191538/146‑1 (MQ=255)
                                                                                                           aTTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCaa                                           >  3:228096/1‑148 (MQ=255)
                                                                                                              gCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTc                                        <  6:173010/148‑1 (MQ=255)
                                                                                                                ccTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGt                                     >  1:221678/1‑149 (MQ=255)
                                                                                                                    cccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGcc                                  <  7:371033/148‑1 (MQ=255)
                                                                                                                      cTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTc                                                                                                     <  4:324338/79‑1 (MQ=255)
                                                                                                                      cTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTc                                                                                                     >  3:324338/1‑79 (MQ=255)
                                                                                                                      cTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCa                               >  5:399560/1‑149 (MQ=255)
                                                                                                                       tGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCa                               >  8:243644/1‑148 (MQ=255)
                                                                                                                           cctaaacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTTGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATGCCt                          >  4:8243/7‑149 (MQ=38)
                                                                                                                           cGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCt                          <  2:202178/149‑1 (MQ=255)
                                                                                                                              aaacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTg                        >  2:29065/4‑148 (MQ=255)
                                                                                                                                ccTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCAc                     >  1:242296/1‑149 (MQ=255)
                                                                                                                                  tgtgGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACgg                   >  7:296043/1‑149 (MQ=255)
                                                                                                                                    tgAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACgg                   <  5:315258/147‑1 (MQ=255)
                                                                                                                                           gggCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCTCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGcc          <  4:228096/149‑1 (MQ=255)
                                                                                                                                                 gCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCAt                              >  5:190772/1‑123 (MQ=255)
                                                                                                                                                 gCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCAt                              <  6:190772/123‑1 (MQ=255)
                                                                                                                                                   cTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTGCACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCg   >  3:57933/1‑148 (MQ=255)
                                                                                                                                                   cTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGt  <  8:212019/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGT  >  NC_002947/5845824‑5846119

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