Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,846,239 T→C 8.6% T217A (ACC→GCC)  mutM ← formamidopyrimidine‑DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,846,2390TC 8.6% 190.8 / 7.1 71T217A (ACC→GCC) mutMformamidopyrimidine‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base C (4/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

GGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCC  >  NC_002947/5846092‑5846384
                                                                                                                                                   |                                                                                                                                                 
ggCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGggtgg                                                                                                                                                 <  1:211413/150‑1 (MQ=255)
     tACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCg                                                                                                                                            >  1:999072/1‑150 (MQ=255)
        tgtactgCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGccc                                                                                                                                         <  2:6857/144‑1 (MQ=255)
          gCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGCGGTACCACCCTg                                                                                                                                       >  2:98553/1‑150 (MQ=255)
           cTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGCGGTACCACCCTGc                                                                                                                                      <  2:119976/150‑1 (MQ=255)
            tCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCt                                                                                                                                     >  2:1165023/1‑150 (MQ=255)
                ccTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGAt                                                                                                                                 <  2:999072/150‑1 (MQ=255)
                ccTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGAt                                                                                                                                 <  1:849815/150‑1 (MQ=255)
                    tGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGgcc                                                                                                                             <  2:941924/150‑1 (MQ=255)
                      ccccAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTa                                                                                                                                               <  1:420608/130‑1 (MQ=255)
                      ccccAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTa                                                                                                                                               >  2:420608/1‑130 (MQ=255)
                             ttCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGc                                                                                                                    >  2:441151/1‑150 (MQ=255)
                                       gcAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGa                                                                                                          <  1:1074425/150‑1 (MQ=255)
                                           cTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGAtctc                                                                                                      >  1:436027/1‑150 (MQ=255)
                                           cTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGAtctc                                                                                                      <  1:653808/150‑1 (MQ=255)
                                                   cGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCa                                                                                              >  1:104661/1‑150 (MQ=255)
                                                   cGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCa                                                                                              >  2:126583/1‑150 (MQ=255)
                                                    gCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAg                                                                                             >  1:663764/1‑150 (MQ=255)
                                                       tACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGCAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGGGGTACCGCCCTGCTAGATTGCCGCCCCCAGCAATCGTTTGATGTAGATCGCCCGttt                                                                                          >  2:539628/1‑150 (MQ=255)
                                                       tACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGttt                                                                                          >  1:799612/1‑150 (MQ=255)
                                                         ccTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCa                                                                                        <  2:109695/150‑1 (MQ=255)
                                                         ccTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCa                                                                                        <  1:441151/150‑1 (MQ=255)
                                                          tttacgcggcggCCCGCCACGCCCGTAGACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAg                                                                                       >  2:17717/7‑150 (MQ=255)
                                                          cTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAg                                                                                       >  1:631936/1‑150 (MQ=255)
                                                                    tGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGAACGCCAGTTACAGGTTGCTCGcc                                                                             >  2:274257/1‑150 (MQ=255)
                                                                        cGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCggg                                                                         >  1:472147/1‑150 (MQ=255)
                                                                            aCGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAAt                                                                     <  2:104661/150‑1 (MQ=255)
                                                                             cGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCgggcattg                                                                    <  2:1285714/150‑6 (MQ=255)
                                                                              gCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcc                                                                   <  2:663764/150‑1 (MQ=255)
                                                                                ccGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccg                                                                 >  2:202554/1‑150 (MQ=255)
                                                                                   tACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccgccc                                                              >  1:4243/1‑150 (MQ=255)
                                                                                   tACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccgccc                                                              >  2:608693/1‑150 (MQ=255)
                                                                                    acacAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTAc                                                                                                                                              >  2:1187621/1‑69 (MQ=255)
                                                                                    acacAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTAc                                                                                                                                              <  1:1187621/69‑1 (MQ=255)
                                                                                      acaAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCcgcccgcc                                                           >  2:245133/1‑150 (MQ=255)
                                                                                             gTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGc                                                    <  2:799612/150‑1 (MQ=255)
                                                                                              ttCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCg                                                   >  1:630409/1‑150 (MQ=255)
                                                                                               tCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGt                                                  <  1:202554/150‑1 (MQ=255)
                                                                                                ccTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTg                                                 <  1:1165023/150‑1 (MQ=255)
                                                                                                 ctgctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTgg                                                <  1:126583/150‑1 (MQ=255)
                                                                                                 ctgctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTAGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccgcc                                                               <  1:1168549/135‑1 (MQ=255)
                                                                                                 ctgctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTAGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccgcc                                                               >  2:1168549/1‑135 (MQ=255)
                                                                                                  tgctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCTCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCg                                                      <  2:872350/143‑1 (MQ=255)
                                                                                                  tgctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCg                                                      >  1:872350/1‑143 (MQ=255)
                                                                                                  tgctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTggg                                               >  2:92577/1‑150 (MQ=255)
                                                                                                       gAAGTACCCCGGCTGCGCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTTGag                                          <  2:78229/150‑2 (MQ=255)
                                                                                                       gAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGAt                                          <  1:274257/150‑1 (MQ=255)
                                                                                                         aGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATAc                                        >  1:206158/1‑150 (MQ=255)
                                                                                                            aCCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACcgg                                     >  2:404181/1‑150 (MQ=255)
                                                                                                              cccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcc                                                                   <  1:628849/118‑1 (MQ=255)
                                                                                                              cccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcc                                                                   >  2:628849/1‑118 (MQ=255)
                                                                                                              cccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACcggcg                                   <  1:245133/150‑1 (MQ=255)
                                                                                                               ccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACcggcgg                                  <  1:92577/150‑1 (MQ=255)
                                                                                                                   cTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCATGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCaaa                              <  1:307000/150‑1 (MQ=255)
                                                                                                                       ccGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGTTCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAaacaac                          <  1:404181/150‑1 (MQ=255)
                                                                                                                        cGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACg                         <  2:693853/150‑1 (MQ=255)
                                                                                                                             ccACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTTGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCg                    >  1:769425/1‑150 (MQ=255)
                                                                                                                             ccACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCg                    >  2:196354/1‑150 (MQ=255)
                                                                                                                             ccACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCCCg                    >  1:397731/1‑150 (MQ=255)
                                                                                                                              cACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCgg                   >  1:1232688/1‑150 (MQ=255)
                                                                                                                               aCCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGGGGTTCGGTACCGGCGGCAAACAACGCCTCGGt                  >  2:529944/1‑150 (MQ=255)
                                                                                                                                 cGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTgg                <  2:630409/150‑1 (MQ=255)
                                                                                                                                 cGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTgg                >  2:522077/1‑150 (MQ=255)
                                                                                                                                   acaaaGTCGCGCAGGGTCGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCNCCTCTCGGCNTGGGTCGANACCGGCGGCAAACAACGCCTCGGTGGcn              <  1:539628/148‑2 (MQ=255)
                                                                                                                                       catcgagcAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGa          >  1:854455/7‑150 (MQ=255)
                                                                                                                                       aGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTGGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGa          >  1:604328/1‑150 (MQ=255)
                                                                                                                                        atcgagcAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTggg                                               >  2:98678/6‑112 (MQ=38)
                                                                                                                                        atcgagcAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTggg                                               <  1:98678/107‑1 (MQ=38)
                                                                                                                                          cgagcAGGGCGGTACCACCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGt       >  2:549273/4‑150 (MQ=255)
                                                                                                                                           gcgcAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGtt      >  2:251209/1‑150 (MQ=255)
                                                                                                                                              cAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGcc   <  1:17717/150‑1 (MQ=255)
                                                                                                                                               aGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGccc  >  1:986648/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCC  >  NC_002947/5846092‑5846384

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: