Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,016 T→C 8.3% intergenic (+32/‑28) secA → / → mutT protein translocation ATPase/8‑oxo‑dGTP diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,0160TC 8.3% 155.5 / 3.6 65intergenic (+32/‑28)secA/mutTprotein translocation ATPase/8‑oxo‑dGTP diphosphatase
Reads supporting (aligned to +/- strand):  ref base T (32/26);  new base C (4/2);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

CTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCG  >  NC_000913/110881‑111155
                                                                                                                                       |                                                                                                                                           
cTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGATAACTGTTGAAGTAAAAGGGGCAGGATTCTGCGcc                                                                                                                                      >  2:238050/1‑143 (MQ=255)
       ggcggcGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTta                                                                                                                               <  1:11881/143‑1 (MQ=255)
       ggcggcGCAAACCGGAGAGCGCAAAGTAGGACGTAACGAGCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTta                                                                                                                               <  2:578196/143‑1 (MQ=255)
            cgcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCGTTGCCCGTGCGGATCTGGTAAGAAATACAAGCAGTGCCGTGGCCGCATGCAATAAAGGCTAACTGTTGAAGTAAAAGGCGCAGGATGCGGCGGCtttattgtcggg                                                                                                                          >  2:455717/1‑137 (MQ=255)
             gcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGtt                                                                                                                         >  1:684471/1‑143 (MQ=255)
             gcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGtt                                                                                                                         >  1:806633/1‑143 (MQ=255)
             gcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGtt                                                                                                                         <  1:884929/143‑1 (MQ=255)
              cAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGACCCTGCGCCCTTTTTTTACGttt                                                                                                                        >  1:567329/1‑143 (MQ=255)
                aaCCGGCGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTaa                                                                                                                      <  1:113077/143‑1 (MQ=255)
                aaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCt                                                                                                                                     >  1:830008/1‑128 (MQ=255)
                aaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCt                                                                                                                                     <  2:830008/128‑1 (MQ=255)
                  ccGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGa                                                                                                                    >  1:394555/1‑143 (MQ=255)
                           gcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATCCAAGCAGTGCCATGGCAGCCTGCAATAAAAGCTAACTGCTGGAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTCTAGGTTTAAGACAATGaaat                                                                                                           >  1:233766/1‑142 (MQ=255)
                            cAAAGTAGGCCGTCGCGTTGCTGGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaaaaa                                                                                                          <  2:789899/143‑1 (MQ=255)
                             aaaGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAg                                                                                                         <  2:522009/143‑1 (MQ=255)
                                 tAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGTCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGTAGGAGTCTGCGCCTTTTTTGTAGGTTCAAGACTATGAAAAAGCTGc                                                                                                     >  1:600798/1‑143 (MQ=255)
                                    gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCCGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTTGAGGTTTAACACCATTAAAAAGCTGCgaa                                                                                                  >  1:337072/1‑143 (MQ=255)
                                    gACGTAACGATCCGTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCGGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCaaa                                                                                                  <  1:373507/143‑1 (MQ=255)
                                       gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTGTGTCAATGAAACACCTCCAAATTg                                                                                               >  1:761791/1‑143 (MQ=255)
                                         aaCGATCCCGGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGCGCCATGTCCGCCTGCAATAAAAGCTTCCTGTTGAAGTAAAAGGCGCAGGACTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                                                             <  1:715301/143‑1 (MQ=255)
                                             aTCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGACGTGCCATGGCCGCCTGCAATAAAAGCTAACTGCGGAAGTAAAGGGCGCAGGAATCTGCACCTTTTTTATAGGATGAAGACAAAGAAACAGCTGAGAAATGCgggag                                                                                         >  2:623531/1‑143 (MQ=255)
                                               ccTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtaggt                                                                                       <  2:828553/143‑1 (MQ=255)
                                                 ttGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtcggtac                                                                                     >  1:466041/1‑142 (MQ=255)
                                                   tcccGTGCGGGTCTTATAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAATTAAAAGACTCAAGATTTTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtatta                                                                                   <  2:366746/142‑1 (MQ=255)
                                                          cGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCaa                                                                            >  1:700293/1‑143 (MQ=255)
                                                            gTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACAGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTAGAGGTTGAAGACAATGAAAAAATTGCAAACTGCGGTAGTGATAATTCGCAACg                                                                          >  2:399488/1‑143 (MQ=255)
                                                             ttCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACgc                                                                         >  1:27407/1‑142 (MQ=255)
                                                                tGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGaa                                                                      >  2:709990/1‑143 (MQ=255)
                                                                 ggTAAAAAATACAAGCACTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAAc                                                                     <  1:149802/143‑1 (MQ=255)
                                                                     aaaaaTACAAGCAGCGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATg                                                                 <  2:250022/143‑1 (MQ=255)
                                                                       aaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAGGCTGCAAATTGCGGTAGGGATTATTCGCAACGAGCACaacgca                                                               >  2:613949/1‑139 (MQ=255)
                                                                       aaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaa                                                               <  2:690577/143‑1 (MQ=255)
                                                                       aaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaa                                                               >  1:283133/1‑143 (MQ=255)
                                                                        aaTACAAGCAGTGCCATGGCCGCCTGCAATCAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATCGGTTTAAGACAAGGAAAACGCTGCACAGCGCGGTAGGTATGATTCGCAACGAGACGAATGaca                                                              >  2:267294/1‑141 (MQ=255)
                                                                        aaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaaa                                                              >  1:55419/1‑143 (MQ=255)
                                                                          tACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATc                                                            <  1:432982/143‑1 (MQ=255)
                                                                             aaGCACTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                         <  2:198002/143‑1 (MQ=255)
                                                                                cAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtata                                                      >  2:159383/1‑143 (MQ=255)
                                                                                   tGCCATTGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcg                                                   <  1:429215/143‑1 (MQ=255)
                                                                                     ccATGGCCGCCTGCAATAAAAGCTAACTGTGGGAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTATAAGACAATGAAAAAGCTGCAAATTGCGGTAGGGATTCTTCGCAACGCGAACGAGGAAATCTGTATAAcgcg                                                 >  2:234921/1‑143 (MQ=255)
                                                                                       aTGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTc                                               <  2:806633/143‑1 (MQ=255)
                                                                                          gccgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGTAGGATCCTGCGCCTTTGTGATAGGTTTAAGACAATGAATAGGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTcgcg                                            <  2:435846/143‑1 (MQ=255)
                                                                                          gccgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTcgcg                                            <  1:374554/143‑1 (MQ=255)
                                                                                            cgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgca                                          <  2:330823/143‑1 (MQ=255)
                                                                                             gccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgcag                                         <  1:439459/143‑1 (MQ=255)
                                                                                              ccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTACGCCTTTTTTATAGGATTAAGACAAGCAAAAAGCTGCCGCTTGCGGTAGGTATTATTAGCAACGAGAAGAATGAAATATTAATAGCACGTCGcacaaa                                        >  2:602897/1‑141 (MQ=255)
                                                                                                   aaTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATg                                   >  1:726069/1‑143 (MQ=255)
                                                                                                    aTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGGGTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATGATTCGCAAGGAGAACAATGAAGTCTTTATAACGCGTCGCGcagcaggggc                                  >  1:31916/1‑139 (MQ=255)
                                                                                                     tAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCAGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGCTgcg                                 >  2:411887/1‑143 (MQ=255)
                                                                                                       acaGCTGACGGTTGAAGGGAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGca                               <  2:344998/141‑1 (MQ=255)
                                                                                                          gCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAAAACAATGAAATCTTTATAACGCGTCGCGAAGCAGATGCGCACAt                            >  2:784523/1‑143 (MQ=255)
                                                                                                                tggtGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGaa                      <  1:898404/140‑1 (MQ=255)
                                                                                                                 gTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAAt                     >  2:911572/1‑143 (MQ=255)
                                                                                                                  ttGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGCATa                    >  2:323302/1‑143 (MQ=255)
                                                                                                                       gTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTCTTATTCTCAACGAGCCGCGTCATTTCGTTAAAATGCTTCGCGCAGCATATACGCACATGGCGAAGAAACta               >  1:385448/1‑142 (MQ=255)
                                                                                                                        tAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAGTCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTgg              >  1:752058/1‑143 (MQ=255)
                                                                                                                        aaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTgg              >  1:248069/2‑143 (MQ=255)
                                                                                                                        aaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGACAAATAAACTgc              >  2:714291/2‑142 (MQ=255)
                                                                                                                         aaaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa             <  2:73830/143‑1 (MQ=255)
                                                                                                                          aaaGGCGCAGGATTTTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAg            <  2:722460/143‑1 (MQ=255)
                                                                                                                               cgcAGGAGTCTGCGCCTTTTTTATGGGTGTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATGCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTcc       <  2:427552/143‑1 (MQ=255)
                                                                                                                                 cAGGATTCTGCGCCTTTTTTATTGGTTTAAGACAATGAAAAAGCTGCAATTTTCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATACGCACATGTCGAATAAACTACAGGTTCCcg     >  2:619884/1‑143 (MQ=255)
                                                                                                                                    gATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGACATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcggcg  <  1:469662/143‑1 (MQ=255)
                                                                                                                                    gATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcggcg  >  2:778279/1‑143 (MQ=255)
                                                                                                                                    gATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGAGCACATGGCGAATAAACTGGAGATTCCcggcg  >  1:243510/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCG  >  NC_000913/110881‑111155

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