Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,013 G→A 17.2% intergenic (+29/‑31) 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,0130GA17.2% 299.5 / 45.4 128intergenic (+29/‑31)secA/mutTpreprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (61/45);  new base A (8/14);  total (69/59)
Fisher's exact test for biased strand distribution p-value = 9.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 130 total aligned reads displayed.

TAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCG  >  NC_000913/110914‑111111
                                                                                                   |                                                                                                  
  ggACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGTTAACTGTTGAAGTAAAAGGCGCAGAAt                                                                                                  <  2:2146465/100‑1 (MQ=255)
  ggACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAt                                                                                                  <  1:2228479/100‑1 (MQ=255)
   gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAtt                                                                                                 >  1:1112568/1‑100 (MQ=255)
   gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAtt                                                                                                 >  2:1166887/1‑100 (MQ=255)
   gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAtt                                                                                                 >  1:2636073/1‑100 (MQ=255)
    aCGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTc                                                                                                <  1:9582/100‑1 (MQ=255)
    aCGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTc                                                                                                >  2:2324958/1‑100 (MQ=255)
    aCGTAACGATCCTTGCACGTGCTGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGAAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTc                                                                                                <  1:810355/100‑1 (MQ=255)
     cGTAACGATCCTTGCTCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTc                                                                                                <  1:2154323/99‑1 (MQ=255)
       tAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                             >  2:936374/1‑100 (MQ=255)
       tAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                             <  2:1791277/100‑1 (MQ=255)
        aaCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCCATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgcg                                                                                            <  2:1652468/100‑1 (MQ=255)
        aaCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgcg                                                                                            >  1:1499797/1‑100 (MQ=255)
        aaCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgcg                                                                                            >  1:1499651/1‑100 (MQ=255)
          cGATCCTTGCCCTTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgcgca                                                                                          <  2:1351740/100‑2 (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)
             tCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGTAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttt                                                                                       <  1:240640/100‑1 (MQ=255)
             tCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCTGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttt                                                                                       <  1:975747/100‑1 (MQ=255)
             tCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttt                                                                                       <  2:2012764/100‑1 (MQ=255)
              ccTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtttt                                                                                      <  2:2147389/100‑1 (MQ=255)
               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)
                 tGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttttact                                                                                   <  2:3035227/100‑4 (MQ=39)
                  gCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAACAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttttactt                                                                                  <  2:2993969/100‑5 (MQ=38)
                   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)
                     cGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAgg                                                                               <  1:2366896/100‑1 (MQ=255)
                         cGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTa                                                                           >  1:1686720/1‑100 (MQ=255)
                         cGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTa                                                                           <  2:2636041/100‑1 (MQ=255)
                          ggTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTaa                                                                          >  1:2361247/1‑100 (MQ=255)
                           gTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAg                                                                         >  1:2402536/1‑100 (MQ=255)
                           gTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGATTAAg                                                                         >  1:542182/1‑100 (MQ=255)
                             tCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                       >  2:1641498/1‑100 (MQ=255)
                             tCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                       <  2:556938/100‑1 (MQ=255)
                              cTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACa                                                                      <  2:544396/100‑1 (MQ=255)
                               tGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACaa                                                                     <  2:735970/100‑1 (MQ=255)
                                   aaaaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaa                                                                 >  1:2193985/1‑100 (MQ=255)
                                         tACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCt                                                           <  1:2099415/100‑1 (MQ=255)
                                           cAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGc                                                         >  2:811761/1‑100 (MQ=255)
                                            aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCa                                                        <  2:1664015/100‑1 (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)
                                                   gCCATGGCCGCCTGCAATAATAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                 >  2:1844506/1‑100 (MQ=255)
                                                   gCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                 >  1:2531150/1‑100 (MQ=255)
                                                   gCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                 >  1:337488/1‑100 (MQ=255)
                                                   gCCATGGCCGCCTGCAATAAAAGCCAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAAt                                                     >  2:1840048/1‑96 (MQ=255)
                                                    ccATGGCCGCCTGCAATAAAATCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCgg                                                <  1:631713/100‑1 (MQ=255)
                                                    ccATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCgg                                                >  2:2021173/1‑100 (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)
                                                             ccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtatta                                       >  1:3022618/1‑100 (MQ=255)
                                                               tGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTattatt                                     >  2:1190173/1‑100 (MQ=255)
                                                                gCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTc                                    >  1:2997048/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)
                                                                     aaaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAAc                               <  1:102522/100‑1 (MQ=255)
                                                                            caatGTTGAAGTACAAGGAGAAGGATTATGCGCATTTTTCGTTTGGTTAACACAACGTAAAAGCTGCAAATGGCGTTCCTTATTATTCGCAACGATAACa                        >  2:2846082/4‑100 (MQ=255)
                                                                               tGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATATGTTTAAGAAAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATg                     <  2:429896/100‑1 (MQ=255)
                                                                               tGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACaa                       <  2:2448203/98‑1 (MQ=255)
                                                                                gTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGa                    >  1:2702608/1‑100 (MQ=255)
                                                                                gTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGa                    >  1:1976487/1‑100 (MQ=255)
                                                                                acctataaaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATAGAGGTAGGTATTATTCACAACAAGAACAAtca                    >  1:2877977/9‑98 (MQ=21)
                                                                                   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)
                                                                                     aGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCt               >  2:2953751/1‑100 (MQ=255)
                                                                                     aGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCt               >  1:2834571/1‑100 (MQ=255)
                                                                                     aGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCt               >  1:714158/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)
                                                                                           aGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATaa         <  1:720865/100‑1 (MQ=255)
                                                                                           aGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATaa         <  2:1205241/100‑1 (MQ=255)
                                                                                            ggCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATAAAATCTTTATAAc        >  1:1610274/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)
                                                                                                   |                                                                                                  
TAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCG  >  NC_000913/110914‑111111

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: