Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,223,854 T→C 18.9% R9R (CGT→CGC 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,8540TC18.9% 225.7 / 50.0 113R9R (CGT→CGCmetHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base T (43/48);  new base C (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 = 7.16e-01

Only 100 of 120 total aligned reads displayed.

CCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTT  >  NC_000913/4223755‑4223952
                                                                                                   |                                                                                                  
ccACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgt                                                                                                    >  1:344732/1‑100 (MQ=255)
 cacaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                                                                                   <  2:1789162/100‑1 (MQ=255)
     gATGCGTTTTTGCCAGTATGGTTTGTTGAATTTTTATGAAATCTGAGGTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGGGGAACAACTGCGTGCGCa                                                                                               <  1:2051119/100‑1 (MQ=255)
     gATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                                                                                               >  2:1767080/1‑100 (MQ=255)
     gATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                                                                                               >  2:166295/1‑100 (MQ=255)
         cGTTTATGTTAGTATGGTGTGTTGAATTTTTCTTAAATCTGGGTTCTGCGGGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGcgcacgcaatag                                                                                           <  1:1887196/100‑5 (MQ=255)
         cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAATGGGGAACGACTGCCTGCGCAGtt                                                                                            >  1:491274/1‑99 (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:15888/1‑100 (MQ=255)
            ttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt                                                                                        <  2:2246854/100‑1 (MQ=255)
            ttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttc                                                                                        <  1:1681236/100‑5 (MQ=255)
            ttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTCGGTTGAGCGTGTCGGTAGCAACTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt                                                                                        <  2:44970/100‑1 (MQ=255)
             tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  1:2127199/100‑6 (MQ=255)
             tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  2:238747/100‑6 (MQ=255)
             tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  2:1274672/100‑6 (MQ=255)
             tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  2:2008125/100‑6 (MQ=255)
             tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                       <  1:1686323/100‑6 (MQ=255)
              aTGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcca                                                                                      <  2:2279453/100‑7 (MQ=255)
              aTGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcca                                                                                      <  1:2559605/100‑7 (MQ=255)
               tGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa                                                                                     <  2:2367375/100‑1 (MQ=255)
                gCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc                                                                                    >  2:2836668/1‑100 (MQ=255)
                gCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttccact                                                                                    <  1:1638008/100‑9 (MQ=255)
                 ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                   <  2:39381/100‑1 (MQ=255)
                     tATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAt                                                                                >  1:722166/1‑99 (MQ=255)
                        ggTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg                                                                            >  2:69642/1‑100 (MQ=255)
                         gTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg                                                                           <  1:2440795/100‑1 (MQ=255)
                            tgttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                        >  1:2889279/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)
                                      ttattaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcggcgg                                                              >  2:1367924/1‑100 (MQ=255)
                                       tattaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGt                                                             >  2:903276/1‑100 (MQ=255)
                                        attaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTa                                                            >  2:2753883/1‑100 (MQ=255)
                                           aaaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATgg                                                         <  2:1111076/100‑1 (MQ=255)
                                           aaaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATgg                                                         <  2:399981/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)
                                               cTGGGTTGGGCGTGTCGGGAGCAAGTGTGAGCGGCAAGGTGGAACAACTGCGTGGGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGc                                                       <  1:1681888/98‑1 (MQ=255)
                                                tGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                    <  2:904115/100‑1 (MQ=255)
                                                  ggTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCAt                                                  <  1:2204324/100‑1 (MQ=255)
                                                   gTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATg                                                 >  1:2637880/1‑100 (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:1683426/100‑1 (MQ=255)
                                                           tgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                         <  1:2754617/100‑1 (MQ=255)
                                                           tgtCGGGAGCAAGTGTGAGCAGCAAAGTGCAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                         >  2:125647/1‑100 (MQ=255)
                                                            gtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGt                                        <  1:1810827/100‑1 (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                            <  2:2619847/100‑1 (MQ=255)
                                                                           gAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACg                           >  1:2405867/1‑98 (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)
                                                                                gcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGa                    >  1:2179871/1‑100 (MQ=255)
                                                                                  aaaGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAtt                  >  1:2006388/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)
                                                                                        atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCCACTGAACGAAGCCGATTTTCgtg            >  1:2762712/5‑100 (MQ=255)
                                                                                         tttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg           >  1:400451/4‑100 (MQ=255)
                                                                                         tttaacTGCGCACGCAGTTAAATGAACGTATTCTAGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAATCCGATTTTCgtgg           >  2:1369875/4‑100 (MQ=255)
                                                                                         aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg           <  1:551139/100‑1 (MQ=255)
                                                                                          acaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtggt          <  2:2858710/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)
                                                                                              cTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAAc      >  1:876247/1‑100 (MQ=255)
                                                                                               tgcgtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACg     <  2:388283/100‑1 (MQ=255)
                                                                                                gcgtgagCAGTTAAATGAACGGCTTCTGATGCTGGACGGAGGTATGGGCACCTTGATCCAGAGTTATCGACTGAACGAAGCCGATTTTAGTGGTGATCGc    <  1:213209/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTT  >  NC_000913/4223755‑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: