Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,223,854 T→C 15.4% 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,8540TC15.4% 352.4 / 51.9 153R9R (CGT→CGCmetHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base T (56/72);  new base C (11/12);  total (67/86)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 168 total aligned reads displayed.

CCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTT  >  NC_000913/4223755‑4223953
                                                                                                   |                                                                                                   
ccACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgt                                                                                                     >  1:2507555/1‑100 (MQ=255)
ccACAGATGCGTTTATGCCAGTATGGATTGTTGAATTTTTATTACATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgt                                                                                                     >  2:202716/1‑100 (MQ=255)
 cacaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                                                                                    <  1:2613194/100‑1 (MQ=255)
 cacaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                                                                                    <  2:3595033/100‑1 (MQ=255)
 cacaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                                                                                    >  2:3389161/1‑100 (MQ=255)
 cacaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                                                                                    >  1:2123736/1‑100 (MQ=255)
   caGATGCGATTATGCCAGTATGGTCTGTAGAATTTTTATTAAATCTGGGTTAAGCGTTTTGGGAGCAAGTGTAAGCAGCAAAGTGGAACCACttcgtg                                                                                                    >  2:684018/1‑98 (MQ=255)
     gATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                                                                                                >  1:475082/1‑100 (MQ=255)
      aTGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAg                                                                                               <  1:3265214/100‑1 (MQ=255)
        gCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAg                                                                                               >  2:250024/1‑98 (MQ=255)
         cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                                                                                            <  1:2820085/100‑1 (MQ=255)
         cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                                                                                            <  1:2153010/100‑1 (MQ=255)
         cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                                                                                            <  1:3280281/100‑1 (MQ=255)
         cGTTTATGCCAGTATGGTCTGTTGAATTTTTATTCAATCTGGGTTGAGCGTGACGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                                                                                            >  2:3506658/1‑100 (MQ=255)
           tttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGTGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGAAAAGTGGAACAACTGCGCACGCAGttgtt                                                                                          <  2:182634/100‑4 (MQ=255)
           tttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgtt                                                                                          <  1:1456131/100‑4 (MQ=255)
             tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                        <  1:3557273/100‑6 (MQ=255)
             tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                        <  2:1123662/100‑6 (MQ=255)
             tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                        <  1:3048729/100‑6 (MQ=255)
               tGCCAGTCTGGTTTGTTAAATTTTAATTAAATCTGGGTTCAGCGTGTCCGGAGAAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttccac                                                                                      <  1:165768/100‑8 (MQ=255)
               tGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGTGTTTAGCGTGTCGGGAGTAAGTGTGAGCCGCAAAGTGGAACCACTGCGCACGCAGttgttccac                                                                                      <  2:3364636/100‑8 (MQ=255)
                gCCAGTCTGGCTTGTAGAATTTTTATTAAATCTGGGTTCAGCGAGTCCGGAGCAAGTGTGAGCACCAAAGTGGAACAACTGCGCACGCAGttgttccac                                                                                      <  1:165783/99‑8 (MQ=255)
                gCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa                                                                                      >  2:1486428/1‑99 (MQ=255)
                 ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                    >  1:2560587/1‑100 (MQ=255)
                 ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                    >  1:718206/1‑100 (MQ=255)
                 ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                    <  1:3520191/100‑1 (MQ=255)
                 ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                    <  1:1597774/100‑1 (MQ=255)
                 ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                    <  1:1522409/100‑1 (MQ=255)
                  cAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGt                                                                                   >  2:2292629/1‑100 (MQ=255)
                   aGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTa                                                                                  <  2:196817/100‑1 (MQ=255)
                      aTGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAttn                                                                               <  2:1826023/100‑2 (MQ=255)
                        ggTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg                                                                             >  1:2183622/1‑100 (MQ=255)
                        ggTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg                                                                             <  2:891930/100‑1 (MQ=255)
                         gTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg                                                                            <  1:957764/100‑1 (MQ=255)
                         gTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg                                                                            <  1:1521272/100‑1 (MQ=255)
                           ttgttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACGGCGTGCGCAGTTAAATGAACGTATTCTGGTg                                                                          <  2:2933303/100‑1 (MQ=255)
                            tgttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                         <  1:1748278/100‑1 (MQ=255)
                            atttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                         <  1:2141320/98‑1 (MQ=255)
                             gttgAATTTTTATTAAATCTGGGTTGCGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                         >  2:2646468/1‑99 (MQ=255)
                               tgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAACAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTgg                                                                      >  2:2644523/1‑100 (MQ=255)
                                 aaTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAc                                                                    <  1:1928771/100‑1 (MQ=255)
                                 aaTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTACGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAc                                                                    >  2:991663/1‑100 (MQ=255)
                                    ttttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcggc                                                                 >  1:2342474/1‑100 (MQ=255)
                                    ttttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcggc                                                                 <  1:3501600/100‑1 (MQ=255)
                                     tttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAGTGAACGTATTCTGGTGCTGGAcggcg                                                                >  1:2734022/1‑100 (MQ=255)
                                     tttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcggcg                                                                <  2:1128024/100‑1 (MQ=255)
                                      ttattaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcggcgg                                                               <  1:1384045/100‑1 (MQ=255)
                                        attaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTa                                                             <  2:1217316/100‑1 (MQ=255)
                                        attaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACCGCGGTa                                                             >  2:1647690/1‑100 (MQ=255)
                                         ttaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTAt                                                            <  1:1033387/100‑1 (MQ=255)
                                         ttaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTAt                                                            >  1:2248436/1‑100 (MQ=255)
                                          taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATg                                                           >  2:662056/1‑100 (MQ=255)
                                          taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTCCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTAt                                                            >  2:2479875/1‑99 (MQ=255)
                                            aaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATggg                                                         >  1:2534884/1‑100 (MQ=255)
                                             aTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATgg                                                          >  2:1255256/1‑98 (MQ=255)
                                             aTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGc                                                        <  1:61240/100‑1 (MQ=255)
                                              tCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCa                                                       <  1:1300618/100‑1 (MQ=255)
                                                tGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                     <  1:1807894/100‑1 (MQ=255)
                                                tGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                     <  1:2644541/100‑1 (MQ=255)
                                                 gggTTGAGCGGGTCGGGAGCAAGTGTATGAAGCAAAGTGGAACAACTGCGAGCGCAGTTAAATGAACGAATTCTGGTGCTGGACGGCGGTATGGGCACCa                                                    <  2:526309/100‑1 (MQ=255)
                                                  ggTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCAt                                                   <  1:1020829/100‑1 (MQ=255)
                                                   gTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATg                                                  >  1:2033066/1‑100 (MQ=255)
                                                    ttGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGa                                                 <  2:2560571/100‑1 (MQ=255)
                                                     tGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGAt                                                >  1:818955/1‑100 (MQ=255)
                                                     tGAGCGTGTCGAGAGCAAGTGTGAGCAGCAAAGTGGAAAAACTGCGTGAGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGAt                                                <  1:1511293/100‑1 (MQ=255)
                                                       aGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATcc                                              <  2:1758997/100‑1 (MQ=255)
                                                       aGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATcc                                              >  2:3060307/1‑100 (MQ=255)
                                                       aGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATcc                                              <  1:140231/100‑1 (MQ=255)
                                                        gCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCa                                             >  2:2781224/1‑100 (MQ=255)
                                                         cGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCag                                            <  2:425466/100‑1 (MQ=255)
                                                           tgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagcg                                          <  1:2451230/100‑3 (MQ=255)
                                                           tgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                          <  2:2776686/100‑1 (MQ=255)
                                                             tCGGGAGCCAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGCACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGtt                                        >  2:725199/1‑100 (MQ=255)
                                                               gggAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACAATGATCCAGGGtt                                        >  2:2948514/1‑98 (MQ=255)
                                                                 gAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCg                                    <  1:2427281/100‑1 (MQ=255)
                                                                  aGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGa                                   >  1:2217323/1‑100 (MQ=255)
                                                                   gCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGAc                                  <  1:2781242/100‑1 (MQ=255)
                                                                     aaGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGAcag                                >  2:2204225/1‑98 (MQ=255)
                                                                      aGTGTGCGCAGAAACGTGGAACCACTGCGAGTGCACTTAAATGAACGTATTCTGTTGCTGGACGGCGGTTTGGGCACCATTATCTAGAGTTATCGACTg                                <  2:141396/99‑1 (MQ=255)
                                                                             gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCATAGTTATCGACTGAACGAAGc                        <  2:2223514/100‑1 (MQ=255)
                                                                             gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGCGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGc                        >  2:1820354/1‑100 (MQ=255)
                                                                                gcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGa                     <  2:2740351/100‑1 (MQ=255)
                                                                                 caAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAt                    <  1:992544/100‑1 (MQ=255)
                                                                                     ttcatttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTc                >  1:3145289/8‑100 (MQ=255)
                                                                                      tcatttaacTGCGCACGCAGTTAAATGAACGTCTTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCg               >  1:3456026/7‑100 (MQ=255)
                                                                                      tGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCg               <  1:2580773/100‑1 (MQ=255)
                                                                                      tGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCg               <  1:1183284/100‑1 (MQ=255)
                                                                                       ggAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgt              >  2:826665/1‑100 (MQ=255)
                                                                                       ggAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgt              <  2:207546/100‑1 (MQ=255)
                                                                                       cctttcacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGCGCTATCGCCTGAACGAAGCCGATTTTCgt              >  2:1806825/7‑100 (MQ=255)
                                                                                       catttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGTGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgt              >  2:1219439/6‑100 (MQ=255)
                                                                                        atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg             >  2:831364/5‑100 (MQ=255)
                                                                                        atttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg             >  1:726643/5‑100 (MQ=255)
                                                                                         aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg            <  1:3131647/100‑1 (MQ=255)
                                                                                          acaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtggt           <  1:813493/100‑1 (MQ=255)
                                                                                              cTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAAc       >  1:1100985/1‑100 (MQ=255)
                                                                                               tgcgtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTAGGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACg      <  2:377699/100‑1 (MQ=255)
                                                                                               tgcgtgcgCAGTTAAAGGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTGATCGACTGAACGAAGCCGATTTTCgtggtg          <  2:2420694/96‑1 (MQ=255)
                                                                                                gcgtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGc     >  1:3205305/1‑100 (MQ=255)
                                                                                                  gtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCtt   >  2:2098441/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTT  >  NC_000913/4223755‑4223953

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: