Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,032,336 G→C 19.1% 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,3360GC19.1% 241.0 / 55.4 106intergenic (‑19/+48)dcm/yedJDNA cytosine methyltransferase/putative HD superfamily phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (38/47);  new base C (9/11);  total (47/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.09e-01

Only 100 of 113 total aligned reads displayed.

AATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTC  >  NC_000913/2032237‑2032435
                                                                                                   |                                                                                                   
aaTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTg                                                                                                     <  2:2266563/100‑1 (MQ=255)
 aTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTga                                                                                                    <  1:282961/100‑1 (MQ=255)
 aTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTga                                                                                                    <  2:2485916/100‑1 (MQ=255)
  tttGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTgag                                                                                                   <  2:2134662/100‑1 (MQ=255)
  gatctagcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTgag                                                                                                   <  2:387987/95‑1 (MQ=255)
    tGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGAt                                                                                                 <  2:2258408/100‑1 (MQ=255)
     gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATg                                                                                                <  1:2493185/100‑1 (MQ=255)
       agcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTgagagatc                                                                                              >  1:387985/1‑96 (MQ=255)
        gcagTTTCTCCAGCATTGCCTGTGCGGTATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGcc                                                                                             <  1:504079/100‑1 (MQ=255)
         cagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCg                                                                                            <  1:1825691/100‑1 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGTTGCCGTCTGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGATATGGcc                                                                                             >  1:354576/1‑98 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                           <  2:234598/100‑2 (MQ=255)
          agTTTCTCCAGCATGGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTGACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                           <  2:214112/100‑2 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                          <  1:1113640/100‑3 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                          <  1:443514/100‑3 (MQ=255)
            tttCTCCATCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtg                                                                                         <  1:1271395/100‑4 (MQ=255)
            tttCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtg                                                                                         <  1:2591836/100‑4 (MQ=255)
            tttCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtg                                                                                         <  2:2159619/100‑4 (MQ=255)
             ttCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTGTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtg                                                                                         <  1:1095060/99‑4 (MQ=255)
              tctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtgaa                                                                                       <  1:700243/100‑6 (MQ=255)
              tctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtgaa                                                                                       <  2:1778846/100‑6 (MQ=255)
               ctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                      <  1:2171144/100‑1 (MQ=255)
               ctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                      <  1:2169981/100‑1 (MQ=255)
                tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                     <  2:1818565/100‑1 (MQ=255)
                tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATGTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                     <  2:1964448/100‑1 (MQ=255)
                 ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTa                                                                                    >  2:1272054/1‑100 (MQ=255)
                 ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTa                                                                                    <  2:800355/100‑1 (MQ=255)
                 ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtgaaatc                                                                                    <  2:1824703/100‑9 (MQ=39)
                   aGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                  >  2:1780916/1‑100 (MQ=255)
                   aGCATTGCCTGTGCGGCAGTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                  <  2:48721/100‑1 (MQ=255)
                     cATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACAt                                                                                <  2:340421/100‑1 (MQ=255)
                      aTTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATa                                                                               >  1:1868498/1‑100 (MQ=255)
                        tGCCTGTGCGGCATTCCCGTTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAt                                                                             <  1:1155191/100‑1 (MQ=255)
                         gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAtt                                                                            >  1:2491919/1‑100 (MQ=255)
                         gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAtt                                                                            <  2:1315881/100‑1 (MQ=255)
                          ccTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATGTATTTTCCTGCATAGATCTACCCAGCCATCCGAGCAGGcc                                                                                             >  1:1413736/1‑82 (MQ=255)
                           cTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTAc                                                                          >  1:2073705/1‑100 (MQ=255)
                                cGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCa                                                                     >  2:134389/1‑100 (MQ=255)
                                       nccGGTGCTGTATGAATCGGATACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACCGGGTAACATAAATACAGCCATTTAg                                                                >  1:1074543/2‑98 (MQ=255)
                                       cccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAggg                                                              >  2:1129180/1‑100 (MQ=255)
                                       cccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAggg                                                              <  2:2212255/100‑1 (MQ=255)
                                              cTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcgcg                                                       >  2:69982/1‑100 (MQ=255)
                                               tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                      <  1:690867/100‑1 (MQ=255)
                                               tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                      >  1:1520562/1‑100 (MQ=255)
                                               tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                      >  1:2148187/1‑100 (MQ=255)
                                                gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCt                                                     >  2:2187753/1‑100 (MQ=255)
                                                gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCt                                                     >  2:249520/1‑100 (MQ=255)
                                                gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCt                                                     <  1:1050624/100‑1 (MQ=255)
                                                   tGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcag                                                  <  2:736849/100‑1 (MQ=255)
                                                    gAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagc                                                 <  2:442420/100‑1 (MQ=255)
                                                     aaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagcc                                                <  1:1864329/100‑1 (MQ=255)
                                                     aaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagcc                                                <  2:2407511/100‑1 (MQ=255)
                                                         ggTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                            <  2:427488/100‑1 (MQ=255)
                                                            tACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGAt                                         <  1:678541/100‑1 (MQ=255)
                                                            tACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATCGCCGACAGGTTAACATAATTACAGCCGTTTAGCGCCACGCGCTCAGCCAGCGGAt                                         >  1:2524212/1‑100 (MQ=255)
                                                             aCTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATg                                        <  1:25121/100‑1 (MQ=255)
                                                              cTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGa                                       <  2:2540523/100‑1 (MQ=255)
                                                               tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGCTGaa                                      >  1:1538245/1‑100 (MQ=255)
                                                                gATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaa                                     >  2:292865/1‑100 (MQ=255)
                                                                gATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaa                                     >  1:361315/1‑100 (MQ=255)
                                                                gATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaa                                     >  2:1640865/1‑100 (MQ=255)
                                                                gATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaa                                     >  1:406320/1‑100 (MQ=255)
                                                                  tatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAc                                   >  2:2519058/1‑100 (MQ=255)
                                                                  tatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCCCTCAGCCAGCGGATGAAAAc                                   <  1:1582098/100‑1 (MQ=255)
                                                                   atatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACg                                  <  1:1441949/100‑1 (MQ=255)
                                                                    tatTTTCCTTCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTCGGGCCACGCGCTCAGCCAGAGGATGAAAACGc                                 >  1:1350048/1‑100 (MQ=255)
                                                                         tCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTa                            >  2:700076/1‑100 (MQ=255)
                                                                          ccTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTAt                           >  2:2112999/1‑100 (MQ=255)
                                                                          ccTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGGTAACATAATTACAGCCATTTAGGGCCACGCGCCCAGCCAGCGGATGAAAGCGCATCTAt                           >  1:428913/1‑100 (MQ=255)
                                                                           cTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATc                          <  1:2397263/100‑1 (MQ=255)
                                                                           cTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATc                          <  1:1684898/100‑1 (MQ=255)
                                                                             gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                        >  1:1759449/1‑100 (MQ=255)
                                                                             gCATAGATTTCACCGGCCATCTGAGATGGCAGACAGGATAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                        >  1:800744/1‑100 (MQ=255)
                                                                              cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                       >  2:509498/1‑100 (MQ=255)
                                                                                 aGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtg                    >  1:1462365/1‑100 (MQ=255)
                                                                                  gATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgt                   >  2:1442439/1‑100 (MQ=255)
                                                                                   aTTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcaaggtg                  <  2:64194/100‑7 (MQ=255)
                                                                                    tttCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATCTAGGGCCACGCGCTCAGCCAGCGGATGCAAACGCATCTATCACCTTGTGa                 >  1:498558/1‑100 (MQ=255)
                                                                                    tttCACCGGCCATCTAAGAAGGACGACAAGTTAACATAATTACATTCATTTAGTGCCACGCGCTCAACCAGCAAATGAAAACGCATCTATCACCTTGTGa                 <  1:1617984/100‑1 (MQ=255)
                                                                                     ttCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt                >  1:1085519/1‑100 (MQ=255)
                                                                                     cctgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt                >  2:1374064/6‑100 (MQ=255)
                                                                                      tCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc               <  1:14334/100‑1 (MQ=255)
                                                                                      tCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc               <  2:2165017/100‑1 (MQ=255)
                                                                                       tgtcGGCCATCTCAGATGGCCGGCAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg              >  1:906017/4‑100 (MQ=255)
                                                                                       tgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg              >  1:1632682/4‑100 (MQ=255)
                                                                                       cACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAACCAGCGGATGAAAACGCATCTATCACCTTGTGATCg              <  1:1129176/100‑1 (MQ=255)
                                                                                        gtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa             >  1:1722051/3‑100 (MQ=255)
                                                                                         cgggCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac            <  1:1569546/98‑1 (MQ=255)
                                                                                          cGGCCATCTGGGATGGCCGACAGGTTAACATATTTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg              <  1:776348/97‑1 (MQ=255)
                                                                                          cGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGaca           <  2:2359624/100‑1 (MQ=255)
                                                                                          cGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGaca           >  2:1193909/1‑100 (MQ=255)
                                                                                           ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGtcac          >  1:985029/1‑100 (MQ=255)
                                                                                            gCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACAcc         <  1:2548897/100‑1 (MQ=255)
                                                                                            gCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACAcc         >  2:2270078/1‑100 (MQ=255)
                                                                                              cATCTGAGATGGCCGACAGGTTAACATAATTAGAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACAcc         <  1:589764/98‑1 (MQ=255)
                                                                                              cATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCtt       <  1:1143911/100‑1 (MQ=255)
                                                                                                tCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCa     <  2:2293048/100‑1 (MQ=255)
                                                                                                  tGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAtt   >  2:1903583/1‑100 (MQ=255)
                                                                                                  tGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAtt   >  1:2014888/1‑100 (MQ=255)
                                                                                                   gagaTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATtc  <  2:162823/100‑1 (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: