Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,032,336 G→C 15.4% 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.4% 279.0 / 40.0 105intergenic (‑19/+48)dcm/yedJDNA cytosine methyltransferase/putative HD superfamily phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (43/45);  new base C (7/9);  total (51/54)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

Only 100 of 111 total aligned reads displayed.

AATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATT  >  NC_000913/2032237‑2032434
                                                                                                   |                                                                                                  
aaTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTg                                                                                                    <  1:498955/100‑1 (MQ=255)
 aTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTga                                                                                                   <  1:929035/100‑1 (MQ=255)
 aTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTga                                                                                                   <  2:1279519/100‑1 (MQ=255)
   ttGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTgaga                                                                                                 <  1:102931/100‑1 (MQ=255)
   ttGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTgaga                                                                                                 <  1:2316568/100‑1 (MQ=255)
     gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATg                                                                                               <  1:1166243/100‑1 (MQ=255)
       agcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGc                                                                                             <  2:77989/100‑1 (MQ=255)
        gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCCGTATGAATCGGTTACTGATATATTTTCCTGCATTGATTCCACCGGCCATCTgaga                                                                                                 >  2:1827346/1‑95 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                          <  1:2458683/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                          <  1:1341738/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                          <  1:2045183/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                          <  2:196474/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                          <  2:1880831/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                          <  2:1421294/100‑2 (MQ=255)
          agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTACATGCATAGATTTCACCGGCCATCTGAGATGGc                                                                                             >  2:273824/1‑97 (MQ=255)
           gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                         <  1:1328782/100‑3 (MQ=255)
                ttccgatctGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                    <  2:1460735/92‑1 (MQ=255)
                tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtgaaat                                                                                    <  1:1367418/100‑8 (MQ=39)
                 ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGGAGGATTCGTTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGGGATGGCCAACAGGt                                                                                     >  1:377647/1‑98 (MQ=255)
                   aGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                 <  1:460403/100‑1 (MQ=255)
                    gCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACa                                                                                <  1:1025164/100‑1 (MQ=255)
                        tGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAt                                                                            <  2:573931/100‑1 (MQ=255)
                         gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTagatcggaa                                                                           >  1:1460746/1‑92 (MQ=255)
                         gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAtt                                                                           >  2:2259497/1‑100 (MQ=255)
                         gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAtt                                                                           <  1:2485019/100‑1 (MQ=255)
                          ccTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTa                                                                          <  2:1118929/100‑1 (MQ=255)
                           cTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTAc                                                                         <  1:2372380/100‑1 (MQ=255)
                           cTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTAc                                                                         <  2:422693/100‑1 (MQ=255)
                               gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGcc                                                                     >  1:2361871/1‑100 (MQ=255)
                               gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGcc                                                                     >  2:2132079/1‑100 (MQ=255)
                               gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGcc                                                                     <  1:1890941/100‑1 (MQ=255)
                               gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGcc                                                                     >  1:2230839/1‑100 (MQ=255)
                                   cATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCAttt                                                                 >  2:568469/1‑100 (MQ=255)
                                    aTTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTa                                                                <  1:2163053/100‑1 (MQ=255)
                                      tCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAgg                                                              <  1:619236/100‑1 (MQ=255)
                                       cccGTTGCTGTATTATTCGGTTACTGATATATTTTACTACATAAATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAggg                                                             <  2:1593885/100‑1 (MQ=255)
                                       cccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAggg                                                             <  2:698088/100‑1 (MQ=255)
                                         cGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAggg                                                             <  2:698110/98‑1 (MQ=255)
                                               tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                     >  2:147501/1‑100 (MQ=255)
                                                gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                     >  2:1141621/1‑99 (MQ=255)
                                                  tctaaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTca                                                  <  1:667452/97‑1 (MQ=255)
                                                    gAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagc                                                <  2:1168800/100‑1 (MQ=255)
                                                     aaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagcc                                               <  2:525589/100‑1 (MQ=255)
                                                     aaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcaaga                                               >  2:667445/1‑97 (MQ=255)
                                                      aTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagcc                                               <  2:186639/99‑1 (MQ=255)
                                                        cGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagccagc                                            >  1:1747053/1‑100 (MQ=255)
                                                         ggTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGCTGACCGACAGATTAACATAAGAACAACCATAAACGCCCACGCGCTCAGCCAACg                                           >  2:850860/1‑100 (MQ=255)
                                                         ggTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagccagc                                            >  1:2440253/1‑99 (MQ=255)
                                                         ggTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                           >  2:451384/1‑100 (MQ=255)
                                                         ggTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                           >  1:2447065/1‑100 (MQ=255)
                                                             aCTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATg                                       >  2:1945295/1‑100 (MQ=255)
                                                             aCTGATATATTTTCCTGCATAGATTTCACCGGCCACCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATg                                       <  1:2538767/100‑1 (MQ=255)
                                                              cTGATATCTTTTCCTGCATAGATTTAACCGGCCATCTAATATGGCCGACAGGTTAACATAATTAAAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATg                                       >  2:1428471/1‑99 (MQ=255)
                                                              cTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGa                                      >  1:930124/1‑100 (MQ=255)
                                                              cTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGa                                      >  2:1407822/1‑100 (MQ=255)
                                                               tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaa                                     <  2:2220852/100‑1 (MQ=255)
                                                                gATATATTTTCCTGCATAGATTTCACCGGCCCTGTGAGATGTCCGACAGGTTAACATAATTACAGCCATTTAGGGACACGCGCTCAGCCAGCGGATGaaa                                    >  2:671272/1‑100 (MQ=255)
                                                                gATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaa                                    >  1:1883354/1‑100 (MQ=255)
                                                                 atgtatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaaa                                   >  1:1157460/4‑100 (MQ=255)
                                                                  tatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAc                                  <  1:891757/100‑1 (MQ=255)
                                                                    tatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACg                                 >  2:376347/1‑99 (MQ=255)
                                                                    tatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGc                                <  1:1644712/100‑1 (MQ=255)
                                                                     atTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCa                               <  1:1999267/100‑1 (MQ=255)
                                                                     atTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCa                               <  2:181321/100‑1 (MQ=255)
                                                                       tttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATc                             <  2:2545015/100‑1 (MQ=255)
                                                                       tttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATc                             <  1:47112/100‑1 (MQ=255)
                                                                        ttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCt                            >  1:2299707/1‑100 (MQ=255)
                                                                         tCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTa                           >  1:1217703/1‑100 (MQ=255)
                                                                         tCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTa                           <  2:1093046/100‑1 (MQ=255)
                                                                           cTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCCTTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTAt                          >  2:49636/1‑99 (MQ=255)
                                                                           cTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATc                         >  1:43546/1‑100 (MQ=255)
                                                                             gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                       <  2:214602/100‑1 (MQ=255)
                                                                             gCATAGAGTTCTCCGGCCATCTGAGATGGCCGACAAGCTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCa                        <  2:342099/99‑1 (MQ=255)
                                                                              cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                      <  2:1982014/100‑1 (MQ=255)
                                                                              cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                      >  2:1414410/1‑100 (MQ=255)
                                                                               aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                      >  2:456708/1‑99 (MQ=255)
                                                                                  gATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgt                  <  2:928832/100‑1 (MQ=255)
                                                                                  gATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgt                  >  2:1905923/1‑100 (MQ=255)
                                                                                  gATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgt                  <  2:1285294/100‑1 (MQ=255)
                                                                                  gATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgt                  <  2:2447051/100‑1 (MQ=255)
                                                                                  gATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgt                  >  1:2268934/1‑100 (MQ=255)
                                                                                     cctgtcGGCCATCTCAGATGGCCGGCAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt               >  1:2264313/6‑100 (MQ=255)
                                                                                     cctgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGCAAACGCATCTATCACCTTGTGAt               >  2:1748098/6‑100 (MQ=255)
                                                                                      ctgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc              >  1:1504899/5‑100 (MQ=255)
                                                                                       tgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg             >  1:1826473/4‑100 (MQ=255)
                                                                                       tgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg             >  1:226207/4‑100 (MQ=255)
                                                                                       cACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg             <  1:911335/100‑1 (MQ=255)
                                                                                       cACCGGCCATCTGACATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg             <  2:1747041/100‑1 (MQ=255)
                                                                                        gtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa            >  1:1055404/3‑100 (MQ=255)
                                                                                        gtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa            >  2:2499251/3‑100 (MQ=255)
                                                                                         ccGGCCTTCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCg             <  2:1403934/98‑1 (MQ=255)
                                                                                         ccGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac           >  1:1250032/1‑100 (MQ=255)
                                                                                         ccGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac           >  1:2535872/1‑100 (MQ=255)
                                                                                          cGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGaca          <  1:258392/100‑1 (MQ=255)
                                                                                           ggCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGacac         >  1:491512/1‑100 (MQ=255)
                                                                                                tCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCa    <  2:1007197/100‑1 (MQ=255)
                                                                                                tCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCCATCACCTTGTGATCGACACCTTCa    <  2:1960547/100‑1 (MQ=255)
                                                                                                 cTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAt   <  1:1959154/100‑1 (MQ=255)
                                                                                                  tGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAt   >  2:1366565/1‑99 (MQ=255)
                                                                                                  tGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGATCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAtt  >  2:1739851/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: