Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,223,855 G→A 17.5% A10T (GCG→ACG)  metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,8550GA17.5% 369.7 / 64.5 161A10T (GCG→ACG) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (60/72);  new base A (15/13);  total (76/85)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

Only 100 of 177 total aligned reads displayed.

CACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTG  >  NC_000913/4223756‑4223954
                                                                                                   |                                                                                                   
cacaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                                                                                     <  2:3956291/100‑1 (MQ=255)
cacaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                                                                                     <  2:15972/100‑1 (MQ=255)
 acaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGAGTGAGCAGCAAAGTGGAACAACtgcgtgc                                                                                                    >  1:56589/1‑100 (MQ=255)
  caGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtgcg                                                                                                   >  2:1147273/1‑100 (MQ=255)
    gATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                                                                                                 >  2:291227/1‑100 (MQ=255)
     aTGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                                                                                                 <  1:2437367/99‑1 (MQ=255)
     aTGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAg                                                                                                <  2:1732763/100‑1 (MQ=255)
     aTGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAg                                                                                                <  2:340853/100‑1 (MQ=255)
      tGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGt                                                                                               >  2:202176/1‑100 (MQ=255)
       gCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGtt                                                                                              <  2:4140222/100‑1 (MQ=255)
       gCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGtt                                                                                              <  2:3833284/100‑1 (MQ=255)
        cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                                                                                             >  2:48385/1‑100 (MQ=255)
        cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                                                                                             <  1:2869225/100‑1 (MQ=255)
        cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                                                                                             >  2:345780/1‑100 (MQ=255)
         gTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgt                                                                                            <  1:926955/100‑3 (MQ=255)
          tttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaaa                                                                                           >  2:4225361/1‑100 (MQ=255)
           ttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt                                                                                          <  2:3914889/100‑1 (MQ=255)
           ttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttc                                                                                          <  2:1304190/100‑5 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGTGTTGAGCGTATCTGGAGCAAGTATGAGAAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt                                                                                          >  1:261267/1‑99 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGCGCAACTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAgctgttcc                                                                                         <  1:2411462/100‑8 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                         <  1:1687580/100‑6 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                         <  2:3625257/100‑6 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                         <  2:1391851/100‑6 (MQ=255)
               gCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc                                                                                      >  2:3667102/1‑100 (MQ=255)
                ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                     <  2:3386382/100‑1 (MQ=255)
                     aTGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTc                                                                                <  2:3096511/100‑1 (MQ=255)
                      tGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCt                                                                               <  1:2144249/100‑1 (MQ=255)
                       ggTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg                                                                              >  2:1160448/1‑100 (MQ=255)
                        gTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg                                                                             <  2:905806/100‑1 (MQ=255)
                           tgttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                          <  1:3319052/100‑1 (MQ=255)
                           tgttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                          <  1:822737/100‑1 (MQ=255)
                            gttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCt                                                                         >  2:1798187/1‑100 (MQ=255)
                              tgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTgg                                                                       >  1:1149103/1‑100 (MQ=255)
                              tgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTg                                                                        >  1:768147/1‑99 (MQ=255)
                                 aTTTTTATTAAATCTGGGTTGAGCGTGTCGGTAGCAAGTGTGAGCAGAAAAGTGGAACGACTGCGTTCGCAGTTAAGTGAACGTATTCTGGTGCTgg                                                                       >  1:1196267/1‑97 (MQ=255)
                                 aTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGGTAAATGAACGTATTCTGGTGCTGGAc                                                                     >  1:3953605/1‑99 (MQ=255)
                                   ttttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATGCTGGTGCTGGAcggc                                                                  >  1:2707068/1‑100 (MQ=255)
                                      tattaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGt                                                               <  1:1505987/100‑1 (MQ=255)
                                         taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATg                                                            >  1:3637204/1‑100 (MQ=255)
                                         taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATg                                                            <  2:4131662/100‑1 (MQ=255)
                                          aaaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATgg                                                           <  2:2735078/100‑1 (MQ=255)
                                          aaaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATgg                                                           >  1:618864/1‑100 (MQ=255)
                                          aaaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATgg                                                           <  2:2732571/100‑1 (MQ=255)
                                          aaaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATg                                                            >  1:2858586/1‑99 (MQ=255)
                                           aaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATggg                                                          <  2:2262811/100‑1 (MQ=255)
                                            aTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGc                                                         <  2:2845527/100‑1 (MQ=255)
                                             tCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCa                                                        <  2:110700/100‑1 (MQ=255)
                                             tCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCa                                                        >  1:3807101/1‑100 (MQ=255)
                                              cTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAc                                                       >  2:277124/1‑100 (MQ=255)
                                               tGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                      <  1:1014397/100‑1 (MQ=255)
                                               tGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                      <  2:144925/100‑1 (MQ=255)
                                                 ggTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                      <  1:2681618/98‑1 (MQ=255)
                                                   ttGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTAGTGCTGGACGGCGGTATGGGCACCATGa                                                  >  2:1476051/1‑100 (MQ=255)
                                                    tGAGCGTGTCGGGAGCAACTGTGAGCAGAAAAGCGGAACAACTGCGTGCGCAGCTACATGAACGTATTCTGGTGCTGGACGGCGGTATGGTCACCATg                                                   >  1:1329205/1‑98 (MQ=255)
                                                     gAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATc                                                <  2:2132105/100‑1 (MQ=255)
                                                      aGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGCGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGAGCTGGACGGCGGTATGGGCACCATGATcc                                               >  1:99451/1‑100 (MQ=255)
                                                       gCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCa                                              <  1:2926234/100‑1 (MQ=255)
                                                         gtgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCaga                                            <  2:2957267/100‑1 (MQ=255)
                                                          tgtCGGGAGCACGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGTCGTCGGTATGGGCACCATGAt                                                 <  1:1523096/94‑1 (MQ=255)
                                                          tgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                           <  2:1666702/100‑1 (MQ=255)
                                                          tgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                           <  1:1807048/100‑1 (MQ=255)
                                                          tgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                           <  2:2591951/100‑1 (MQ=255)
                                                           gtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGt                                          <  2:2643243/100‑1 (MQ=255)
                                                               ggAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTAt                                       <  1:2384106/99‑1 (MQ=255)
                                                               ggAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATc                                      >  2:1875466/1‑100 (MQ=255)
                                                               ggAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATc                                      <  1:1630563/100‑1 (MQ=255)
                                                                gAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCg                                     <  1:2500299/100‑1 (MQ=255)
                                                                gAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCg                                     <  2:3777266/100‑1 (MQ=255)
                                                                gAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCg                                     <  2:3297774/100‑1 (MQ=255)
                                                                  gCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGAc                                   <  2:2444208/100‑1 (MQ=255)
                                                                    aaGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTg                                 <  2:2526901/100‑1 (MQ=255)
                                                                      gtgtGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGaa                               >  1:4183662/1‑100 (MQ=255)
                                                                        gtgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATAGACTGAACg                             >  1:3452547/1‑100 (MQ=255)
                                                                         tgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGa                            >  2:297258/1‑100 (MQ=255)
                                                                         tgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGa                            <  1:2373630/100‑1 (MQ=255)
                                                                         tgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGa                            >  2:4047029/1‑100 (MQ=255)
                                                                           agcagcAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAg                          <  2:65068/100‑1 (MQ=255)
                                                                            gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGc                         <  2:3267953/100‑1 (MQ=255)
                                                                               gcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGa                      >  2:670229/1‑100 (MQ=255)
                                                                               gcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAAGGCGGTATGAGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGa                      >  1:2973971/1‑100 (MQ=255)
                                                                                 aaaGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAtt                    <  2:4333045/100‑1 (MQ=255)
                                                                                  cgttcatttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAttt                   >  2:4201746/10‑100 (MQ=255)
                                                                                  aaGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAttt                   <  2:3482895/100‑1 (MQ=255)
                                                                                     tcatttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCg                >  2:1568770/7‑100 (MQ=255)
                                                                                      ggAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgt               >  1:703727/1‑100 (MQ=255)
                                                                                      catttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgt               >  2:854786/6‑100 (MQ=255)
                                                                                      catttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATAGACTGAACGAAGCCGATTTTCgt               >  1:3630603/6‑100 (MQ=255)
                                                                                       atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg              >  1:3468984/5‑100 (MQ=255)
                                                                                       atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg              >  2:706284/5‑100 (MQ=255)
                                                                                       atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg              >  1:3893587/5‑100 (MQ=255)
                                                                                        tttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg             >  2:1620026/4‑100 (MQ=255)
                                                                                        aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg             >  1:534741/1‑100 (MQ=255)
                                                                                        aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg             >  1:3672314/1‑100 (MQ=255)
                                                                                          caacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtggtg           >  1:321827/1‑100 (MQ=255)
                                                                                            acTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGaa         >  1:3100585/1‑100 (MQ=255)
                                                                                              tgcgtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtggtagatc       >  1:1338863/1‑95 (MQ=255)
                                                                                               gcgtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGc      <  1:2733211/100‑1 (MQ=255)
                                                                                                cgtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCt     >  1:454701/1‑100 (MQ=255)
                                                                                                  tgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCttt   >  2:3643580/1‑100 (MQ=255)
                                                                                                   gcgcAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTg  <  1:1088298/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTG  >  NC_000913/4223756‑4223954

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: