Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,223,855 G→A 18.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,8550GA18.5% 252.5 / 44.2 113A10T (GCG→ACG) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (44/48);  new base A (10/11);  total (54/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 = 5.33e-01

Only 100 of 119 total aligned reads displayed.

CACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTT  >  NC_000913/4223756‑4223952
                                                                                                   |                                                                                                 
    gATGCGTTTTTGCCAGTATGGTTTGTTGAATTTTTATGAAATCTGAGGTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGGGGAACAACTGCGTGCGCa                                                                                               <  1:2051119/100‑1 (MQ=255)
    gATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                                                                                               >  2:1767080/1‑100 (MQ=255)
    gATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                                                                                               >  2:166295/1‑100 (MQ=255)
       gCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGtt                                                                                            <  1:2840288/100‑1 (MQ=255)
        cGTTTATGTTAGTATGGTGTGTTGAATTTTTCTTAAATCTGGGTTCTGCGGGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGcgcacgcaatag                                                                                           <  1:1887196/100‑5 (MQ=255)
        cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                                                                                           <  2:2880496/100‑1 (MQ=255)
        cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                                                                                           >  1:2014428/1‑100 (MQ=255)
          tttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaaa                                                                                         <  2:1399570/100‑1 (MQ=255)
          tttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgtt                                                                                         <  2:1315625/100‑4 (MQ=255)
           ttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt                                                                                        <  2:2246854/100‑1 (MQ=255)
           ttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTCGGTTGAGCGTGTCGGTAGCAACTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt                                                                                        <  2:44970/100‑1 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  2:2008125/100‑6 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  1:2127199/100‑6 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  2:1190999/100‑6 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  2:1274672/100‑6 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  1:1686323/100‑6 (MQ=255)
             aTGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcca                                                                                      <  1:2559605/100‑7 (MQ=255)
              tGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa                                                                                     <  2:2367375/100‑1 (MQ=255)
               gCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc                                                                                    >  2:2836668/1‑100 (MQ=255)
               gCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc                                                                                    <  1:2073924/100‑1 (MQ=255)
               gCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttccact                                                                                    <  1:1638008/100‑9 (MQ=255)
                ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                   <  2:39381/100‑1 (MQ=255)
                       ggTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg                                                                            >  2:69642/1‑100 (MQ=255)
                        gTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg                                                                           <  1:2440795/100‑1 (MQ=255)
                          ttgttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTg                                                                         >  1:1944382/1‑100 (MQ=255)
                            gttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCt                                                                       >  1:305129/1‑100 (MQ=255)
                             ttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTg                                                                      >  1:2533061/1‑100 (MQ=255)
                                aaTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACATCTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTggcc                                                                   >  1:2290617/1‑98 (MQ=255)
                                aaTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAc                                                                   >  1:185360/1‑100 (MQ=255)
                                 aTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcg                                                                  >  2:451076/1‑100 (MQ=255)
                                   ttttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCTCATTTAAATGAACGTATTCTGGTGCTGGAcggc                                                                <  1:1921100/100‑1 (MQ=255)
                                   ttttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcggc                                                                >  1:1170674/1‑100 (MQ=255)
                                      tattaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGt                                                             >  2:903276/1‑100 (MQ=255)
                                       attaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTa                                                            >  2:2753883/1‑100 (MQ=255)
                                       attaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTa                                                            >  1:371525/1‑100 (MQ=255)
                                          aaaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATgg                                                         <  2:399981/100‑1 (MQ=255)
                                          aaaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATgg                                                         <  2:1111076/100‑1 (MQ=255)
                                          aaaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAGTCTGGTGCTGGACGGCGGTATgg                                                         <  1:1699535/100‑1 (MQ=255)
                                           aaTCTGGGTTGAGCGTGTCGGGAGCGAGTGTGAGCAGCACAGTGGAACCACTGCGAGCGCAGTTAAATGAACGTATTATGGTGCTGGAc                                                                   >  1:2459000/1‑89 (MQ=255)
                                           aaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATggg                                                        <  2:641578/100‑1 (MQ=255)
                                            aTCTGGGTTGAGCGTGTCGGAAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGc                                                       <  2:2344229/100‑1 (MQ=255)
                                             tCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGCACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCa                                                      >  1:773171/1‑100 (MQ=255)
                                             tCTGGGTCGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCa                                                      <  2:1932422/100‑1 (MQ=255)
                                              cTGGGTTGGGCGTGTCGGGAGCAAGTGTGAGCGGCAAGGTGGAACAACTGCGTGGGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGc                                                       <  1:1681888/98‑1 (MQ=255)
                                               tGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                    <  2:904115/100‑1 (MQ=255)
                                                 ggTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCAt                                                  <  1:2204324/100‑1 (MQ=255)
                                                   ttGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGa                                                <  2:1579121/100‑1 (MQ=255)
                                                    tGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGAt                                               <  2:1500574/100‑1 (MQ=255)
                                                        cGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCag                                           <  1:2227035/100‑1 (MQ=255)
                                                        cGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGGATGGGCACCATGATCCag                                           >  1:2852813/1‑100 (MQ=255)
                                                         gtgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCaga                                          <  1:451074/100‑1 (MQ=255)
                                                          tgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                         >  1:705070/1‑100 (MQ=255)
                                                          tgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                         <  1:2754617/100‑1 (MQ=255)
                                                          tgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                         <  1:1683426/100‑1 (MQ=255)
                                                          tgtCGGGAGCAAGTGTGAGCAGCAAAGTGCAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                         >  2:125647/1‑100 (MQ=255)
                                                               ggAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATc                                    <  1:2033737/100‑1 (MQ=255)
                                                                gAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCg                                   <  2:1490351/100‑1 (MQ=255)
                                                                    aaGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTg                               <  1:873684/100‑1 (MQ=255)
                                                                     aGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGa                              >  2:2680624/1‑100 (MQ=255)
                                                                       tgtgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAAc                            >  1:1650117/1‑100 (MQ=255)
                                                                       tgtgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAAc                            <  1:2551283/100‑1 (MQ=255)
                                                                       tgtgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAAc                            <  2:2619847/100‑1 (MQ=255)
                                                                        gtgAGCAGCAAAGTGGAACAACTGGGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGa                              >  1:770613/1‑97 (MQ=255)
                                                                          gAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGaa                         >  1:12628/1‑100 (MQ=255)
                                                                          gAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGaa                         >  1:20338/1‑100 (MQ=255)
                                                                          gAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGaa                         >  2:1745800/1‑100 (MQ=255)
                                                                           agcagcAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAg                        >  1:484482/1‑100 (MQ=255)
                                                                           agagggaaaGTGGAACAACTGGGTGGGCAGTTAGATGAACGTAGTCTGGTGCTGGGCGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAg                        <  1:2821886/94‑1 (MQ=255)
                                                                            gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGc                       >  1:1000731/1‑100 (MQ=255)
                                                                            gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGc                       >  2:1674720/1‑100 (MQ=255)
                                                                              agcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCg                     <  2:2533076/100‑1 (MQ=255)
                                                                               gcaAAGTGGAACAACTGCGTGCGCAGTTAAATGCACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGa                    >  1:1916150/1‑100 (MQ=255)
                                                                               gcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCg                     >  2:2513501/1‑99 (MQ=255)
                                                                               gcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGa                    >  1:2179871/1‑100 (MQ=255)
                                                                                 aaaGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAtt                  >  1:1036666/1‑100 (MQ=255)
                                                                                   gttcatttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACTAAGCCGAtttt                >  1:73103/9‑100 (MQ=255)
                                                                                     tcatttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCg              >  2:299342/7‑100 (MQ=255)
                                                                                      ggAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgt             <  2:2014439/100‑1 (MQ=255)
                                                                                      ggAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgt             <  2:2660763/100‑1 (MQ=255)
                                                                                      catttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGAACCAGAGTTATCGACTGAACGAAGCCGATTTTCgt             >  1:820086/6‑100 (MQ=255)
                                                                                       gAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtn            <  1:621958/100‑2 (MQ=255)
                                                                                       gAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg            <  2:2080956/100‑1 (MQ=255)
                                                                                       gAACAACCGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGTTATGGGCACCATGATCCAGAGTTATTGACTGAACGAAGCCGATTTTCg              <  1:296538/98‑1 (MQ=255)
                                                                                       atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg            >  2:234418/5‑100 (MQ=255)
                                                                                       atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg            >  2:2912886/5‑100 (MQ=255)
                                                                                       atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg            >  2:1625997/5‑100 (MQ=255)
                                                                                       atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCCACTGAACGAAGCCGATTTTCgtg            >  1:2762712/5‑100 (MQ=255)
                                                                                        tttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg           >  1:400451/4‑100 (MQ=255)
                                                                                        tttaacTGCGCACGCAGTTAAATGAACGTATTCTAGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAATCCGATTTTCgtgg           >  2:1369875/4‑100 (MQ=255)
                                                                                        aacaacTGCGTGCGCAGTTAAATGAACGTTTTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCTGAGTTATCGACTGAACGAAGCCGATTTTCgtgg           <  1:1745792/100‑1 (MQ=255)
                                                                                        aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg           <  1:551139/100‑1 (MQ=255)
                                                                                        aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg           >  2:1394183/1‑100 (MQ=255)
                                                                                         acaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtggt          <  2:2858710/100‑1 (MQ=255)
                                                                                          caacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtggtg         <  2:22970/100‑1 (MQ=255)
                                                                                           aacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGa        <  2:2302018/100‑1 (MQ=255)
                                                                                            acTGCGTGCGCAGTTAAATGAACGCATTCTGATGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGACGATTTTCGTGGTGa        >  1:580543/1‑99 (MQ=255)
                                                                                             cTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAAc      <  2:672975/100‑1 (MQ=255)
                                                                                              tgcgtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACg     <  2:388283/100‑1 (MQ=255)
                                                                                              tgcgtgcgCAGGTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGCTTTTCGTGGTGAACg     >  1:2235742/1‑100 (MQ=255)
                                                                                                 gtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCtt  >  1:57379/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTT  >  NC_000913/4223756‑4223952

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: