Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,016 T→C 17.0% intergenic (+32/‑28) secA → / → mutT preprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,0160TC17.0% 257.3 / 49.0 124intergenic (+32/‑28)secA/mutTpreprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (57/46);  new base C (8/13);  total (65/59)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-01

Only 100 of 129 total aligned reads displayed.

GACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCA  >  NC_000913/110917‑111115
                                                                                                   |                                                                                                   
gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAtt                                                                                                     >  2:1166887/1‑100 (MQ=255)
gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAtt                                                                                                     >  1:1112568/1‑100 (MQ=255)
gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAtt                                                                                                     >  1:2636073/1‑100 (MQ=255)
 aCGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTc                                                                                                    >  2:2324958/1‑100 (MQ=255)
 aCGTAACGATCCTTGCACGTGCTGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGAAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTc                                                                                                    <  1:810355/100‑1 (MQ=255)
  cGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCt                                                                                                   >  1:2927570/1‑100 (MQ=255)
    tAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                                 >  2:936374/1‑100 (MQ=255)
    tAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                                 <  2:1791277/100‑1 (MQ=255)
    tAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCAGCCTGCAATCAAAGCGAACTGTTGGAGTAAAAGGCGCAGGATTCTg                                                                                                  >  2:2045674/1‑99 (MQ=38)
     aaCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCCATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgcg                                                                                                <  2:1652468/100‑1 (MQ=255)
     aaCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgcg                                                                                                >  1:1499797/1‑100 (MQ=255)
       cGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGcc                                                                                              <  2:3030478/100‑1 (MQ=255)
        gATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCt                                                                                             <  1:1177784/100‑1 (MQ=255)
         aTCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtn                                                                                            <  2:164773/100‑2 (MQ=255)
          ttcTTGTCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtn                                                                                            <  1:1328740/97‑2 (MQ=255)
          tCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCTGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttt                                                                                           <  1:975747/100‑1 (MQ=255)
          tCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttt                                                                                           <  2:2012764/100‑1 (MQ=255)
           ccTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtttt                                                                                          <  2:2147389/100‑1 (MQ=255)
            cTTGTCCGTGCGGTTTTGGTAAAAAATACAAGTAGTGTCATGGCCGCCTGCAATAAAAGCTAACTGTTGTAGTAAAAGGCGCAGAATCCTGCGCCtttt                                                                                          <  1:2761252/99‑1 (MQ=38)
            cTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAATCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttttt                                                                                         <  2:401731/100‑1 (MQ=255)
            cTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCACCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttttt                                                                                         <  1:866249/100‑1 (MQ=255)
             ttGCCCGGGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAATGCTAACTGTTGAAGTAAAAGGCGTAGAATCCTGCGCCttttac                                                                                        <  2:1337756/100‑3 (MQ=255)
              tGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTta                                                                                       >  1:849921/1‑100 (MQ=255)
              tGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTta                                                                                       <  1:2194623/100‑1 (MQ=255)
              tGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttttact                                                                                       <  2:3035227/100‑4 (MQ=39)
               gCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAACAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttttactt                                                                                      <  2:2993969/100‑5 (MQ=38)
                cccTTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtat                                                                                      <  1:2598515/99‑1 (MQ=255)
                cccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttttttn                                                                                       <  1:2396240/98‑2 (MQ=255)
                cccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtata                                                                                     <  1:2980362/100‑1 (MQ=255)
                cccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttttacttc                                                                                     <  1:1746066/100‑6 (MQ=38)
                 ccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtatan                                                                                    <  2:2095459/100‑2 (MQ=255)
                 ccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAg                                                                                    <  2:1956077/100‑1 (MQ=255)
                  cGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAgg                                                                                   <  1:2366896/100‑1 (MQ=255)
                   gTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGTTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtat                                                                                      >  2:1202716/1‑96 (MQ=255)
                      cGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTa                                                                               <  2:2636041/100‑1 (MQ=255)
                        gTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAg                                                                             >  1:2402536/1‑100 (MQ=255)
                        gTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGATTAAg                                                                             >  1:542182/1‑100 (MQ=255)
                          tCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                           <  2:556938/100‑1 (MQ=255)
                          tCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                           >  2:1641498/1‑100 (MQ=255)
                           cTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACa                                                                          <  2:544396/100‑1 (MQ=255)
                            tGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACaa                                                                         <  2:735970/100‑1 (MQ=255)
                               tAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGa                                                                      <  2:1785959/100‑1 (MQ=255)
                                aaaaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaa                                                                     >  1:2193985/1‑100 (MQ=255)
                                   aaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaaaaa                                                                  >  1:842721/1‑100 (MQ=255)
                                    aaTACAAGCAGTGCCATGGCCTCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAg                                                                 <  1:2385381/100‑1 (MQ=255)
                                      tACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCt                                                               <  1:2099415/100‑1 (MQ=255)
                                        cAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGc                                                             >  2:811761/1‑100 (MQ=255)
                                         aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAAGAGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCa                                                            >  1:787970/1‑100 (MQ=255)
                                          aGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAATTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCa                                                            <  2:712582/99‑1 (MQ=255)
                                          aGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGACTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCaa                                                           >  1:708708/1‑100 (MQ=255)
                                            cAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAAt                                                         >  2:514300/1‑100 (MQ=255)
                                                gCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGca                                                     >  1:317384/1‑99 (MQ=255)
                                                gCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                     >  1:337488/1‑100 (MQ=255)
                                                gCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                     <  2:2984543/100‑1 (MQ=255)
                                                gCCATGGCCGCCTGCAATAAAAGCCAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAAt                                                         >  2:1840048/1‑96 (MQ=255)
                                                 ccATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCgg                                                    >  1:1034414/1‑100 (MQ=255)
                                                  cATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAATACAATGAAAAAGCTGCAAATTGCggt                                                   <  2:1968680/100‑1 (MQ=255)
                                                  cATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCgg                                                    <  2:2688795/99‑1 (MQ=255)
                                                    tGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAATGCGCAGGATTCTGCGACTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggta                                                  <  1:2563116/99‑1 (MQ=255)
                                                    tGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtag                                                 <  2:2248391/100‑1 (MQ=255)
                                                     ggCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtagg                                                >  1:2012089/1‑100 (MQ=255)
                                                      gccgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTAGAGGTTTCAGACAATGAAAAAGCTGCAAATTGCggtaggt                                               >  2:2132925/1‑100 (MQ=255)
                                                       ccgccTGCAATAAAAGCTAACAGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGATTTAAGACAATGAAAAAGCTGAAAATTGCggtaggta                                              >  2:1010275/1‑100 (MQ=255)
                                                          ccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtatta                                           >  1:3022618/1‑100 (MQ=255)
                                                            tGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAAAATGAAAAAGCTGCAAATTGCGGTAGGTattatt                                         >  1:1583441/1‑100 (MQ=255)
                                                              cAATAAAATCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCg                                       <  2:600642/100‑1 (MQ=255)
                                                              cAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCg                                       <  2:945433/100‑1 (MQ=255)
                                                                aTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCa                                     >  1:2734535/1‑100 (MQ=255)
                                                                 tAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCaa                                    >  1:1855518/1‑100 (MQ=255)
                                                                  aaaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAAc                                   <  1:102522/100‑1 (MQ=255)
                                                                     aGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACgag                                >  2:1420474/1‑100 (MQ=255)
                                                                         caatGTTGAAGTACAAGGAGAAGGATTATGCGCATTTTTCGTTTGGTTAACACAACGTAAAAGCTGCAAATGGCGTTCCTTATTATTCGCAACGATAACa                            >  2:2846082/4‑100 (MQ=255)
                                                                            tGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATATGTTTAAGAAAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATg                         <  2:429896/100‑1 (MQ=255)
                                                                            tGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACaa                           <  2:2448203/98‑1 (MQ=255)
                                                                             acctataaaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATAGAGGTAGGTATTATTCACAACAAGAACAAtca                        >  1:2877977/9‑98 (MQ=21)
                                                                              cctataaaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAAAAATGaa                       >  1:650298/8‑100 (MQ=37)
                                                                                tataaaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAAt                     >  1:1301087/6‑100 (MQ=38)
                                                                                gAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAAt                     >  2:2553826/1‑100 (MQ=255)
                                                                                gAAGTAAAAGGCGCAGGATTCTGCCCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAAt                     <  2:2337896/100‑1 (MQ=255)
                                                                                 aaGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATAAAAAAGCTACAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATc                    >  1:765163/1‑100 (MQ=255)
                                                                                   gTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAAtat                   >  1:2201115/1‑97 (MQ=255)
                                                                                   gTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCtt                  >  2:720861/1‑100 (MQ=255)
                                                                                   gTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCtt                  >  1:1522971/1‑100 (MQ=255)
                                                                                    aaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                 >  1:1783654/2‑100 (MQ=255)
                                                                                     aaaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGGCAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                <  1:647463/100‑1 (MQ=255)
                                                                                     aaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                >  1:2513188/1‑100 (MQ=255)
                                                                                     aaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                <  2:1157598/100‑1 (MQ=255)
                                                                                      aaattcgcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATTAAAAAGCTTCAAATTGCGGTAGGTAATATTTTCAACGAGAACAATGAAATCtttct               >  2:413433/6‑98 (MQ=255)
                                                                                      aaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTACGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtat               >  1:2664977/1‑100 (MQ=255)
                                                                                       aaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtata              >  1:2505755/1‑100 (MQ=255)
                                                                                        aGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATaa             <  2:1205241/100‑1 (MQ=255)
                                                                                        aGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATaa             <  1:720865/100‑1 (MQ=255)
                                                                                         ggCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATAAAATCTTTATAAc            >  1:1610274/1‑100 (MQ=255)
                                                                                          gcgcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcg           >  1:2102326/1‑100 (MQ=255)
                                                                                          gcgcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcg           >  2:2697306/1‑100 (MQ=255)
                                                                                            gcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcgcg         >  1:1105019/1‑100 (MQ=255)
                                                                                             cAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGt        <  1:2768635/100‑1 (MQ=255)
                                                                                                gATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTcg      <  1:2108073/99‑1 (MQ=255)
                                                                                                  ttCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCgcgc   >  1:2460877/1‑100 (MQ=255)
                                                                                                   tCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgca  <  2:2668961/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCA  >  NC_000913/110917‑111115

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