Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,223,855 G→A 15.1% 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,8550GA15.1% 385.5 / 47.8 153A10T (GCG→ACG) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (55/74);  new base A (11/12);  total (67/86)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 168 total aligned reads displayed.

CACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTG  >  NC_000913/4223756‑4223954
                                                                                                   |                                                                                                   
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)
 acaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtgc                                                                                                    <  2:2701586/100‑1 (MQ=255)
  caGATGCGATTATGCCAGTATGGTCTGTAGAATTTTTATTAAATCTGGGTTAAGCGTTTTGGGAGCAAGTGTAAGCAGCAAAGTGGAACCACttcgtg                                                                                                     >  2:684018/1‑98 (MQ=255)
    gATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                                                                                                 >  1:475082/1‑100 (MQ=255)
     aTGCTTTTATGCCATTATGGTTTGTTGAGGTTTTATTAAGCCTGAGTTCAGCGTGTCGGGAGCAAGTGTGCGCAGCAAAGTTGATCCACTGCGTGCGCAg                                                                                                <  2:2479019/100‑1 (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)
          tttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgtt                                                                                           <  1:2397851/100‑4 (MQ=255)
           ttATGCCAGTATGGTTTGTTGAATTTTTAGTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttc                                                                                          <  2:2572319/100‑5 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                         <  1:3557273/100‑6 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                         <  1:3048729/100‑6 (MQ=255)
            tATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgttcc                                                                                         <  2:1123662/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)
                ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGTAACAACTGCGTGCGCAGTTAAATGAAc                                                                                      >  2:2090645/1‑99 (MQ=255)
                ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                     >  1:718206/1‑100 (MQ=255)
                ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                     >  1:2560587/1‑100 (MQ=255)
                ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                     <  1:1522409/100‑1 (MQ=255)
                ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                     <  1:1597774/100‑1 (MQ=255)
                ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                     <  1:3520191/100‑1 (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)
                           tgttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                          <  1:1748278/100‑1 (MQ=255)
                           atttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                          <  1:2141320/98‑1 (MQ=255)
                            gttgAATTTTTATTAAATCTGGGTTGCGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                          >  2:2646468/1‑99 (MQ=255)
                            gttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCt                                                                         <  1:2342169/100‑1 (MQ=255)
                              tgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTgg                                                                       >  2:2298688/1‑100 (MQ=255)
                                aaTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAc                                                                     <  1:1928771/100‑1 (MQ=255)
                                aaTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTACGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAc                                                                     >  2:991663/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)
                                           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)
                                               tGGGTTGAGCGTGTCGGGAGCAAGTGTGCGCAGCAAAGTGCAACAACTGCGTGCGCAGTTAAATGAACGGATTCTGGTGCTGGACGGCGGTATGGCCAcc                                                      <  2:2511018/100‑1 (MQ=255)
                                               tGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGTACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                      <  2:2947323/100‑1 (MQ=255)
                                               tGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                      <  1:2644541/100‑1 (MQ=255)
                                               tGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCAcc                                                      <  1:1807894/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)
                                                    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)
                                                           gtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGCGt                                          >  2:2784351/1‑100 (MQ=255)
                                                            tCGGGAGCCAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGCACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGtt                                         >  2:725199/1‑100 (MQ=255)
                                                               ggAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGCACGTATTCTGGTGCTGGACGGCGGTATAGGCACCATGATCCAGAGTTATc                                      >  2:2245941/1‑100 (MQ=255)
                                                                gAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCg                                     <  1:2427281/100‑1 (MQ=255)
                                                                 aGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGa                                    >  1:2217323/1‑100 (MQ=255)
                                                                  gCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGAc                                   <  1:2781242/100‑1 (MQ=255)
                                                                     aGTGTGCGCAGAAACGTGGAACCACTGCGAGTGCACTTAAATGAACGTATTCTGTTGCTGGACGGCGGTTTGGGCACCATTATCTAGAGTTATCGACTg                                 <  2:141396/99‑1 (MQ=255)
                                                                            gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCATAGTTATCGACTGAACGAAGc                         <  2:2223514/100‑1 (MQ=255)
                                                                            gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGCGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGc                         >  2:1820354/1‑100 (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:207546/100‑1 (MQ=255)
                                                                                      ggAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgt               >  2:826665/1‑100 (MQ=255)
                                                                                      ggAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCAAGAGTTATCGACTGAACGAAGCCGATTTTc                 >  2:2753885/1‑98 (MQ=255)
                                                                                      cctttcacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGCGCTATCGCCTGAACGAAGCCGATTTTCgt               >  2:1806825/7‑100 (MQ=255)
                                                                                      catttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGTGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgt               >  2:1219439/6‑100 (MQ=255)
                                                                                       gAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg              <  1:2493597/100‑1 (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)
                                                                                               gcgtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAAc        >  2:2118589/1‑98 (MQ=255)
                                                                                               gcgtgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGc      >  1:3205305/1‑100 (MQ=255)
                                                                                                   gcgcAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTg  >  2:3068587/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: