Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,846,246 G→A 10.2% G214G (GGC→GGT mutM ← formamidopyrimidine‑DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,846,2460GA10.2% 175.4 / 9.3 69G214G (GGC→GGTmutMformamidopyrimidine‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base A (4/3);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

TACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCCACCCCTACCAC  >  NC_002947/5846097‑5846395
                                                                                                                                                     |                                                                                                                                                     
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)
                        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                                                                         >  2:608693/1‑150 (MQ=255)
                                                                              tACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccgccc                                                                         >  1:4243/1‑150 (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                                                                          >  2:1168549/1‑135 (MQ=255)
                                                                                            ctgctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTAGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccgcc                                                                          <  1:1168549/135‑1 (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                                                                              >  2:628849/1‑118 (MQ=255)
                                                                                                         cccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcc                                                                              <  1:628849/118‑1 (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)
                                                                                                                                              cggtACCACCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCCAccc         <  1:549273/149‑1 (MQ=255)
                                                                                                                                                    cGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccgcccg                                                                        >  2:434274/1‑81 (MQ=255)
                                                                                                                                                    cGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccgcccg                                                                        <  1:434274/81‑1 (MQ=255)
                                                                                                                                                     gCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCCACCCCTaccac  <  1:608693/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCCACCCCTACCAC  >  NC_002947/5846097‑5846395

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: