Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,846,246 G→A 7.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,2460GA 7.2% 240.2 / 5.8 83G214G (GGC→GGTmutMformamidopyrimidine‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (37/40);  new base A (3/3);  total (40/43)
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.70e-01

TACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCCACCCCTACC  >  NC_002947/5846097‑5846393
                                                                                                                                                     |                                                                                                                                                    
tACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCg                                                                                                                                                      >  2:1317610/1‑150 (MQ=255)
   acGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGccc                                                                                                                                                   >  1:1808639/1‑150 (MQ=255)
     gCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTg                                                                                                                                                 <  1:1394660/150‑1 (MQ=255)
      cTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGCGGTACCACCCTGc                                                                                                                                                <  1:1324795/150‑1 (MQ=255)
       tCGCCGTCTGCCCCAGCTTCACTTCTCGCAACTCAGTGCCGCATACCTGGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCt                                                                                                                                               <  2:988077/150‑1 (MQ=255)
           ccTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGAt                                                                                                                                           <  2:222117/150‑1 (MQ=255)
           ccTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCCCGAt                                                                                                                                           >  2:82178/1‑150 (MQ=255)
            ctctGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATg                                                                                                                                          <  1:2209548/150‑1 (MQ=255)
              ctGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGgc                                                                                                                                        <  2:1950587/150‑1 (MQ=255)
                 ccccAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGgccgc                                                                                                                                     >  1:333017/1‑150 (MQ=255)
                   ccAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATg                                                                                                                                          <  2:2085218/143‑1 (MQ=255)
                   ccAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATg                                                                                                                                          >  1:2085218/1‑143 (MQ=255)
                     aGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATg                                                                                                                                          >  2:112150/1‑141 (MQ=255)
                     aGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATg                                                                                                                                          <  1:112150/141‑1 (MQ=255)
                             ttCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCg                                                                                                                         >  2:2048637/1‑150 (MQ=255)
                             ttCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCg                                                                                                                         >  1:2097323/1‑150 (MQ=255)
                               cgcgCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGtt                                                                                                                       <  2:2343668/150‑1 (MQ=255)
                                gcgcAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGttt                                                                                                                      <  2:2097323/150‑1 (MQ=255)
                                    aaCTCAGTTCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACGAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGGCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGAtc                                                                                                                  <  2:1456421/150‑1 (MQ=255)
                                      cTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGAtctc                                                                                                                >  2:1616393/1‑150 (MQ=255)
                                       tCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCg                                                                                                               >  1:1502207/1‑150 (MQ=255)
                                        cAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGa                                                                                                              <  1:2048637/150‑1 (MQ=255)
                                         aGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGAt                                                                                                             >  1:705431/1‑150 (MQ=255)
                                              ggtaTACCGTGCACGGCGGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCa                                                                                                        <  2:745936/147‑1 (MQ=255)
                                                   aCCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTc                                                                                                   <  1:1631634/150‑1 (MQ=255)
                                                     cTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGGTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGACGCACAGGGTGGTACCGCCCGACTCGAGAGCCCCCGCCAGAACCCGATTGATGTCGAGCGCCAGTTTCAg                                                                                                 >  2:719810/1‑150 (MQ=255)
                                                      ttGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGttt                                                                                                                      <  2:2114484/128‑1 (MQ=255)
                                                      ttGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGttt                                                                                                                      >  1:2114484/1‑128 (MQ=255)
                                                          aCGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGgcc                                                                                                                                       >  1:667176/1‑107 (MQ=255)
                                                          aCGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGgcc                                                                                                                                       <  2:667176/107‑1 (MQ=255)
                                                          aCGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAgc                                                                                            <  2:314250/150‑1 (MQ=255)
                                                            ggCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGcgc                                                                                                                                              >  1:2299205/1‑96 (MQ=255)
                                                            ggCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTc                                                                                                                                              <  2:2299205/98‑1 (MQ=255)
                                                                  ccGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGGCCGCCGCCACCGATACAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCTCGCCCgg                                                                                    <  1:265066/150‑1 (MQ=255)
                                                                  ccGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCg                                                                                                               <  1:2195376/123‑1 (MQ=255)
                                                                  ccGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCg                                                                                                               >  2:2195376/1‑123 (MQ=255)
                                                                  ccGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCgg                                                                                    <  2:1003940/150‑1 (MQ=255)
                                                                  ccGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGACCTCGATCGCCAGTTTCAGCTAGCGCGCCCgg                                                                                    >  2:2195945/1‑150 (MQ=255)
                                                                   cGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGAtct                                                                                                                 >  2:1710763/1‑120 (MQ=255)
                                                                   cGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGAtct                                                                                                                 <  1:1710763/120‑1 (MQ=255)
                                                                    gCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGa                                                                                  <  2:1808639/150‑1 (MQ=255)
                                                                        cGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATc                                                                              <  2:333017/150‑1 (MQ=255)
                                                                           ccGTACACAAACAGTTCCTGCTGGACGGANGCCGGCTGCCCGTCGNCNCCGATAAAGTCGCGCAGGGTGGTANCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccg                                                                           >  1:572926/1‑150 (MQ=255)
                                                                             gTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcc                                                                             <  1:1898005/146‑1 (MQ=255)
                                                                             gTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcc                                                                             >  2:1898005/1‑146 (MQ=255)
                                                                                   aaaCAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCctc                                                                   >  2:1488440/1‑150 (MQ=255)
                                                                                      cAGTTCCTGCTGGAACTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCg                                                                <  1:2044059/150‑1 (MQ=255)
                                                                                              gctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGt                                                        >  1:589857/1‑150 (MQ=255)
                                                                                               ctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGt                                                        >  1:2189955/1‑149 (MQ=255)
                                                                                               ctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGt                                                        <  2:2189955/149‑1 (MQ=255)
                                                                                               ctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCCCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTc                                                       <  1:82178/150‑1 (MQ=255)
                                                                                                  gAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGAt                                                    >  2:917979/1‑150 (MQ=255)
                                                                                                      tACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACcg                                                <  1:1092780/150‑1 (MQ=255)
                                                                                                         cccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACcggcg                                             <  1:1317610/150‑1 (MQ=255)
                                                                                                          ccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTa                                                                                              <  1:886586/100‑1 (MQ=255)
                                                                                                          ccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTa                                                                                              >  2:886586/1‑100 (MQ=255)
                                                                                                          ccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGc                                                              >  2:822922/1‑132 (MQ=255)
                                                                                                          ccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGc                                                              <  1:822922/132‑1 (MQ=255)
                                                                                                          ccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACcggcgg                                            >  2:512724/1‑150 (MQ=255)
                                                                                                            ggCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCa                                          >  2:1248774/1‑150 (MQ=255)
                                                                                                              cTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCaaa                                        <  1:1478664/150‑1 (MQ=255)
                                                                                                                 cccGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTgg                                                          >  1:1002924/1‑129 (MQ=255)
                                                                                                                 cccGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTgg                                                          <  2:1002924/129‑1 (MQ=255)
                                                                                                                 cccGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGTTCGATACCGGCGGCAaacag                                     <  2:705431/150‑2 (MQ=255)
                                                                                                                 cccGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAaacaa                                     <  1:917979/150‑1 (MQ=255)
                                                                                                                    gTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGc                                  <  1:1248774/150‑1 (MQ=255)
                                                                                                                        ccACCGTGAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCg                              <  2:1812470/150‑1 (MQ=255)
                                                                                                                        ccACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCg                              >  1:1735917/1‑150 (MQ=255)
                                                                                                                          aCCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGt                            >  2:655390/1‑150 (MQ=255)
                                                                                                                          aCCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCACCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGt                            >  1:154123/1‑150 (MQ=255)
                                                                                                                                  catCGAGCAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGa                    <  2:2407501/148‑1 (MQ=255)
                                                                                                                                  catCGAGCAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGa                    >  1:2407501/3‑150 (MQ=255)
                                                                                                                                   gTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGAt                   >  1:1317876/1‑150 (MQ=255)
                                                                                                                                    tcgagcAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTg                           <  2:1587015/137‑1 (MQ=255)
                                                                                                                                    tCGAGCAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTg                           >  1:1587015/1‑141 (MQ=255)
                                                                                                                                      gcgcAGGGGGGTACCGCCCGGCTCGATGGCCGCCCCCAGCGCCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAAGCCCGCCCGCCTCTCGGCGTGGGTCGAGACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGtt                <  1:799645/150‑1 (MQ=255)
                                                                                                                                        gcAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGc              >  2:888574/1‑150 (MQ=255)
                                                                                                                                         cAGGGTGATACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGGTTGGAGGCAGTGGTCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTGTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGcc             <  2:572926/150‑1 (MQ=255)
                                                                                                                                                gtACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAACCAACGCCTCGGTGGCGTAGATGTTGCCCACCCCt     >  2:901524/1‑150 (MQ=255)
                                                                                                                                                gtACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTACCCACCCCt     >  2:2334503/1‑150 (MQ=255)
                                                                                                                                                 tACCGCCCTGCTCGATGGCCGCCGACACTGCGCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACCACGCCTCGGTGGCGTAGATGTTGCCGCAcccta     >  1:63556/1‑148 (MQ=255)
                                                                                                                                                   ccGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTc                                                                                                   >  1:914842/1‑54 (MQ=255)
                                                                                                                                                   ccGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTc                                                                                                   <  2:914842/54‑1 (MQ=255)
                                                                                                                                                   ccGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGCTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCCACCCCTacc  >  2:1401305/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCCACCCCTACC  >  NC_002947/5846097‑5846393

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