Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,032,336 G→C 20.0% 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,3360GC20.0% 389.5 / 93.0 165intergenic (‑19/+48)dcm/yedJDNA cytosine methyltransferase/putative HD superfamily phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (68/64);  new base C (18/15);  total (86/79)
Fisher's exact test for biased strand distribution p-value = 8.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 174 total aligned reads displayed.

AATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTC  >  NC_000913/2032237‑2032435
                                                                                                   |                                                                                                   
   ttGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTgaga                                                                                                  <  2:3911996/100‑1 (MQ=255)
     gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATg                                                                                                >  1:4014984/1‑100 (MQ=255)
     gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATg                                                                                                >  2:2274865/1‑100 (MQ=255)
     gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCGGCATAGATTTCACCGGCCATCTGAGATg                                                                                                <  2:1815127/100‑1 (MQ=255)
      cagcagTTTTTCCGGTATTGCCTGTTCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACAAGCCATCTGAGATgg                                                                                               <  1:3542124/100‑1 (MQ=255)
      cagcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATgg                                                                                               >  1:2142139/1‑100 (MQ=255)
        gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGcc                                                                                             >  1:3342767/1‑100 (MQ=255)
        gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGcc                                                                                             >  2:2066784/1‑100 (MQ=255)
        gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGcc                                                                                             <  2:550823/100‑1 (MQ=255)
        gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGcc                                                                                             <  1:3553300/100‑1 (MQ=255)
        gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGcc                                                                                             <  2:4130565/100‑1 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                           <  2:2798345/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                           <  2:196349/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                           <  2:1769390/100‑2 (MQ=255)
           gtgtctcCAGCATGGCGTGTGCGGCATTCCCGGGGATGTATGAATCGGTTACTGCTATACTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                          <  1:2566486/97‑3 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGAc                                                                                          <  2:1146902/100‑1 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGAc                                                                                          >  1:776795/1‑100 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                          <  2:3905540/100‑3 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                          <  2:2652472/100‑3 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                          <  2:2658417/100‑3 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                          <  2:2828203/100‑3 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                          <  1:4263046/100‑3 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTAATGATATATTTTCCTGCATAGATTTCACCGACCATCTGAAATGGCCGAc                                                                                          >  1:3925013/1‑100 (MQ=255)
             ttCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAg                                                                                        <  1:2320044/100‑1 (MQ=255)
              tctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAgg                                                                                       >  1:1172705/1‑100 (MQ=255)
              tctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAgg                                                                                       <  2:2042183/100‑1 (MQ=255)
              tctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtgaa                                                                                       <  2:1270491/100‑6 (MQ=255)
               ntcCAGCATTTCCTGTGCGGCATTCCCGGTGCCGTCTGAATCGGTTACTGATATATGTCCCTGCATAGATTTTACCGGCCATCTgcgctg                                                                                                >  1:3215142/2‑85 (MQ=255)
               ctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAgg                                                                                       >  1:980920/1‑99 (MQ=255)
               ctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                      >  2:2605186/1‑100 (MQ=255)
               ctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                      >  1:1469/1‑100 (MQ=255)
               ctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                      >  2:3850387/1‑100 (MQ=255)
                 ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtgaaatc                                                                                    <  2:3684068/100‑9 (MQ=39)
                   aGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTCATATCTTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTa                                                                                    <  1:2991375/98‑1 (MQ=255)
                    gCATTGCCTGTGCGGCATTCCCGTTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACa                                                                                 <  1:3491190/100‑1 (MQ=255)
                    gCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACa                                                                                 >  1:114423/1‑100 (MQ=255)
                     cATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACAt                                                                                <  1:125962/100‑1 (MQ=255)
                     atctGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACAt                                                                                <  1:472318/97‑1 (MQ=255)
                         gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATagat                                                                            >  2:472320/1‑97 (MQ=255)
                          ccTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTa                                                                           >  2:2492953/1‑100 (MQ=255)
                            tgtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                         >  1:330537/1‑100 (MQ=255)
                             gtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAg                                                                        >  2:3478584/1‑100 (MQ=255)
                                cGGCATTCCCGGTGCTGTATTAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGc                                                                       <  1:4213402/98‑1 (MQ=255)
                                  gCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCAtt                                                                   >  1:2855948/1‑100 (MQ=255)
                                    aTTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTa                                                                 <  2:2079818/100‑1 (MQ=255)
                                       cccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAggg                                                              <  1:1005089/100‑1 (MQ=255)
                                         cGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGcc                                                            <  2:4139646/100‑1 (MQ=255)
                                            tGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcg                                                         <  2:2908132/100‑1 (MQ=255)
                                              cTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcgcg                                                       <  2:4362039/100‑1 (MQ=255)
                                              cTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcgcg                                                       >  1:1561480/1‑100 (MQ=255)
                                               tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                      <  2:1305561/100‑1 (MQ=255)
                                                gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCt                                                     >  2:3423682/1‑100 (MQ=255)
                                                 tATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTc                                                    <  2:658749/100‑1 (MQ=255)
                                                  aTGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTca                                                   >  1:1113384/1‑100 (MQ=255)
                                                   tGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcag                                                  <  2:902228/100‑1 (MQ=255)
                                                    gAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagc                                                 >  2:3555111/1‑100 (MQ=255)
                                                     aaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagcc                                                >  1:1664959/1‑100 (MQ=255)
                                                         ggTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                            <  2:1357546/100‑1 (MQ=255)
                                                         ggTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                            <  1:3782850/100‑1 (MQ=255)
                                                          gTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCgg                                           <  1:2211016/100‑1 (MQ=255)
                                                          gTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCgg                                           >  1:3858304/1‑100 (MQ=255)
                                                            tACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGAt                                         <  2:230753/100‑1 (MQ=255)
                                                              cTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGTGATGGCCGACTTTTTAACATAATTACTGCCATTTAGGGCCACGCGCTCAGCCAGCGGATg                                        <  1:2618700/99‑1 (MQ=255)
                                                                   atatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACg                                  <  1:4152525/100‑1 (MQ=255)
                                                                    tatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGc                                 >  2:4059179/1‑100 (MQ=255)
                                                                      ttttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCAt                               >  1:3825329/1‑100 (MQ=255)
                                                                       tttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATc                              >  1:1859372/1‑100 (MQ=255)
                                                                       tttCCTGCATACATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCa                                <  1:2051864/98‑1 (MQ=255)
                                                                            tGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCa                         <  2:590127/100‑1 (MQ=255)
                                                                            tGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCa                         >  2:2670695/1‑100 (MQ=255)
                                                                              cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                       <  2:215568/100‑1 (MQ=255)
                                                                              cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                       <  1:3346079/100‑1 (MQ=255)
                                                                              cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                       >  1:2316727/1‑100 (MQ=255)
                                                                                 aGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtg                    <  2:2533282/100‑1 (MQ=255)
                                                                                     cctgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt                >  2:1571090/6‑100 (MQ=255)
                                                                                     cctgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt                >  2:339930/6‑100 (MQ=255)
                                                                                     cctgacGGCCGTCGCAGATGGCCGACAGGTTAACATGCGGAGAGGCAAGTAGGGCCACGAGCTCAGCAGGCGGATGACAACGCATCTATCAGCAtgtgc                 >  1:289065/6‑98 (MQ=255)
                                                                                     catgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCGACGCGCTCAGACAACGGAAGAAAACGCATCTATCACGTTGTGAt                >  1:2584851/6‑100 (MQ=255)
                                                                                      tCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc               <  2:3967639/100‑1 (MQ=255)
                                                                                      tCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc               <  1:1349605/100‑1 (MQ=255)
                                                                                      ctgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc               >  2:1163111/5‑100 (MQ=255)
                                                                                       cACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTGAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg              <  1:2713318/100‑1 (MQ=255)
                                                                                       cACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg              >  2:375215/1‑100 (MQ=255)
                                                                                       cACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg              >  1:3424951/1‑100 (MQ=255)
                                                                                        gtcGGCCATCTCAGATGGCCGACAGGTTCACATAATTACAGCCATTTATCGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCGa             >  1:2355420/3‑100 (MQ=255)
                                                                                        gtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa             >  1:1055585/3‑100 (MQ=255)
                                                                                         ccGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc               <  1:2441795/97‑1 (MQ=255)
                                                                                          cGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGaca           >  2:1670807/1‑100 (MQ=255)
                                                                                           ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGaca           >  2:3430024/1‑99 (MQ=255)
                                                                                           ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCACGCAGCCAGCGGATGAAAACACATCTATCCCCTTATGATCGacac          >  1:3739965/1‑100 (MQ=255)
                                                                                            gCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACAcc         <  2:1208287/100‑1 (MQ=255)
                                                                                              tctCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCtt       <  2:2600135/98‑1 (MQ=255)
                                                                                              tctCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCtt       <  2:2454389/98‑1 (MQ=255)
                                                                                              cATCTGAGATGGCTGACAGGTTAACATACTTACAGGCATTGAGGGTCACGCGCCCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc               <  1:2766853/92‑1 (MQ=255)
                                                                                                tCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCa     <  2:3643049/100‑1 (MQ=255)
                                                                                                 cTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCttaga    >  1:2454377/1‑97 (MQ=255)
                                                                                                 cTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAt    >  2:1035678/1‑100 (MQ=255)
                                                                                                 cTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAt    <  2:3042843/100‑1 (MQ=255)
                                                                                                   gagaTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGAATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATtc  >  1:2038685/1‑100 (MQ=255)
                                                                                                   gagaTGGCCGACAGGGTCACATAAATACAGCCATTTAGGGCTACGCGCTCAGCCAGCGTAAGAAACCGCATCTATCACCTCGCGATCGACACCTTCATtc  >  1:901350/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTC  >  NC_000913/2032237‑2032435

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: