Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,032,336 G→C 15.3% intergenic (‑19/+48) dcm ← / ← yedJ DNA cytosine methyltransferase/putative HD superfamily phosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,032,3360GC15.3% 391.4 / 58.9 144intergenic (‑19/+48)dcm/yedJDNA cytosine methyltransferase/putative HD superfamily phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (61/60);  new base C (11/11);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 150 total aligned reads displayed.

AATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATT  >  NC_000913/2032237‑2032434
                                                                                                   |                                                                                                  
aaTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTg                                                                                                    <  1:3334161/100‑1 (MQ=255)
aaTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTg                                                                                                    <  2:1974060/100‑1 (MQ=255)
aaTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTg                                                                                                    <  1:1599810/100‑1 (MQ=255)
     gcaggaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATg                                                                                               <  1:3621084/100‑1 (MQ=255)
      cagcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATgg                                                                                              >  2:1194893/1‑100 (MQ=255)
      cagcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATCTCACCGGCCCTCTgagct                                                                                                >  1:3558655/1‑96 (MQ=255)
        gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGcc                                                                                            >  2:2094409/1‑100 (MQ=255)
         cagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCg                                                                                           <  2:3708670/100‑1 (MQ=255)
         cagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCg                                                                                           >  2:82408/1‑100 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGa                                                                                          >  2:1916268/1‑100 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGa                                                                                          <  2:284849/100‑1 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                          <  1:882193/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                          <  1:163614/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                          <  2:1171039/100‑2 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGAc                                                                                         <  2:303796/100‑1 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                         <  1:3513943/100‑3 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                         <  1:2827427/100‑3 (MQ=255)
            tttCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtg                                                                                        <  1:3246113/100‑4 (MQ=255)
             ttCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAg                                                                                       >  1:1653887/1‑100 (MQ=255)
             ttCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtga                                                                                       <  1:727505/100‑5 (MQ=255)
             ttCTCCAGCATTGCCTGTACGACATTACCTGTGCTTTATGAATCGGTTACTTATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACa                                                                                        <  2:2689651/99‑1 (MQ=255)
              tctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtgaa                                                                                      <  1:2754271/100‑6 (MQ=255)
                tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                    >  2:2989761/1‑100 (MQ=255)
                 ncAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTa                                                                                   >  1:2400650/2‑100 (MQ=255)
                 ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTa                                                                                   <  1:3014375/100‑1 (MQ=255)
                 ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtgaaatc                                                                                   <  1:2122491/100‑9 (MQ=39)
                  cAGCATTGCCTGTGCGGCAATCCCAGTGCTGTATGAAGCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                    <  1:3798763/98‑1 (MQ=255)
                   aGCATTGCCTGTGCGGCATTCCGGGTGCGGTATGAATCGGTTACTGATACATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                 <  2:3382993/100‑1 (MQ=255)
                   aGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                 >  2:3400800/1‑100 (MQ=255)
                   aGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                 <  1:722988/100‑1 (MQ=255)
                    gCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                 >  1:1887689/1‑99 (MQ=255)
                    gCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACa                                                                                <  1:1624495/100‑1 (MQ=255)
                    gCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACa                                                                                >  2:3253214/1‑100 (MQ=255)
                    gCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACa                                                                                <  1:1973136/100‑1 (MQ=255)
                        tGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAt                                                                            >  1:2191706/1‑100 (MQ=255)
                        tGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTGACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAt                                                                            <  2:3007004/100‑1 (MQ=255)
                         gCCTGTGCGGCATTCCCGTTGCTGTATGAAGCGGTTACTGATCTATTTTCCTGCATTGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACTTAAtt                                                                           >  1:10322/1‑100 (MQ=255)
                         gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAtt                                                                           >  1:2741560/1‑100 (MQ=255)
                         gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAtt                                                                           >  2:2199018/1‑100 (MQ=255)
                             gtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATGTCACCGGCCATGTGAGATGGCCGACAGGTTAACATAATTACa                                                                        >  2:755377/1‑99 (MQ=255)
                              tgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGc                                                                      <  1:2548512/100‑1 (MQ=255)
                               gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGcc                                                                     >  1:1513797/1‑100 (MQ=255)
                                cGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCa                                                                    <  1:268925/100‑1 (MQ=255)
                                cGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCa                                                                    <  1:1357585/100‑1 (MQ=255)
                                   cATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGGTGGCCGACAGGTTAACATAATTACAGCCa                                                                    >  2:2334311/1‑97 (MQ=255)
                                   cATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCAttt                                                                 <  1:2794803/100‑1 (MQ=255)
                                     ttCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                               <  1:1545622/100‑1 (MQ=255)
                                     ttCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                               <  1:3006757/100‑1 (MQ=255)
                                            tGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcg                                                        >  2:1996341/1‑100 (MQ=255)
                                             gCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcgc                                                       <  2:2582581/100‑1 (MQ=255)
                                              cTGTATGAATCGGTTACTGATATATTTTCTTGCATCGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcg                                                        <  2:2572757/98‑1 (MQ=255)
                                               tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                     >  1:3637498/1‑100 (MQ=255)
                                               tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                     >  2:2842936/1‑100 (MQ=255)
                                               tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                     >  2:1826441/1‑100 (MQ=255)
                                                gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCt                                                    >  1:881355/1‑100 (MQ=255)
                                                gTAAGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCt                                                    >  1:2542924/1‑100 (MQ=255)
                                                 tATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTc                                                   >  2:53964/1‑100 (MQ=255)
                                                    gAATCAGTTACTGATATATTTTCCTGGATAGATTTAACCGGCCATCTTAGAGTGCCGACCGGTTTACATAATTACAACCATTTAGGGCCAAGCGCTca                                                  >  2:908102/1‑98 (MQ=255)
                                                    gAACCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagc                                                >  1:1119687/1‑100 (MQ=255)
                                                     aaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagcc                                               <  2:1570945/100‑1 (MQ=255)
                                                          gTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCgg                                          <  1:3389101/100‑1 (MQ=255)
                                                           ttACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGa                                         >  1:505568/1‑100 (MQ=255)
                                                              cTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGa                                      >  1:362927/1‑100 (MQ=255)
                                                                  tatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAc                                  >  1:572480/1‑100 (MQ=255)
                                                                    tatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGc                                <  2:330853/100‑1 (MQ=255)
                                                                     atTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCa                               >  2:2771032/1‑100 (MQ=255)
                                                                     atTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCa                               >  1:3259706/1‑100 (MQ=255)
                                                                       tttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATc                             <  2:3760670/100‑1 (MQ=255)
                                                                       tttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATc                             <  2:1863676/100‑1 (MQ=255)
                                                                       tttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATc                             <  1:3203296/100‑1 (MQ=255)
                                                                         tCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTa                           <  1:1056630/100‑1 (MQ=255)
                                                                           cTGCATGGATTTCACCGGCCATCTGAGATGGCCGACACGTTAACACAATTACAGCCACTTAGGGCCACCCGCTCAGCCAGCTGATGACAACCCAACCATc                         >  2:3271787/1‑100 (MQ=255)
                                                                             gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                       <  1:1847105/100‑1 (MQ=255)
                                                                             gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                       <  1:2703323/100‑1 (MQ=255)
                                                                             gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                       >  2:679077/1‑100 (MQ=255)
                                                                             gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                       <  1:364920/100‑1 (MQ=255)
                                                                             gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                       >  2:1932372/1‑100 (MQ=255)
                                                                              cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                      >  1:218101/1‑100 (MQ=255)
                                                                                tAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCtt                    <  2:234675/100‑1 (MQ=255)
                                                                                 aGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtg                   <  2:1383385/100‑1 (MQ=255)
                                                                                 aGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtg                   <  2:1614759/100‑1 (MQ=255)
                                                                                  gATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgt                  <  1:1175259/100‑1 (MQ=255)
                                                                                     ttCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt               <  2:1086200/100‑1 (MQ=255)
                                                                                     cctgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt               >  1:2310660/6‑100 (MQ=255)
                                                                                       tgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg             >  1:1897128/4‑100 (MQ=255)
                                                                                       tgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAACTTGTGATCg             >  1:3178983/4‑100 (MQ=255)
                                                                                       cACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg             <  1:977853/100‑1 (MQ=255)
                                                                                        gtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg             >  2:518629/3‑99 (MQ=255)
                                                                                        gtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa            >  2:1457570/3‑100 (MQ=255)
                                                                                         tcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac           >  2:3361254/2‑100 (MQ=255)
                                                                                         ccGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac           >  1:3770003/1‑100 (MQ=255)
                                                                                         ccGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac           >  1:1945561/1‑100 (MQ=255)
                                                                                          cGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGaca          <  1:1704022/100‑1 (MQ=255)
                                                                                          cGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac           <  2:3302334/99‑1 (MQ=255)
                                                                                           ggCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGacac         >  2:2175533/1‑100 (MQ=255)
                                                                                           ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGaca          >  2:3440945/1‑99 (MQ=255)
                                                                                            gCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACAcc        <  2:3155610/100‑1 (MQ=255)
                                                                                              cATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCtt      >  1:1821586/1‑100 (MQ=255)
                                                                                                tCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCGCGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCa    >  1:2202536/1‑100 (MQ=255)
                                                                                                  tGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAtt  >  2:2095384/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
AATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATT  >  NC_000913/2032237‑2032434

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: