Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,223,595 C→T A730V (GCC→GTC)  metH → B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,5950CT99.4% 601.9 / NA 178A730V (GCC→GTC) metHB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (102/75);  minor base G (1/0);  total (103/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 191 total aligned reads displayed.

GGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGCCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCAC  >  NC_000913/4223504‑4223683
                                                                                           |                                                                                        
 gCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCt                                                                                        >  1:1510462/1‑93 (MQ=255)
 gCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCt                                                                                        >  1:4263477/1‑93 (MQ=255)
  cATGAATGTGGTGGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTa                                                                                       <  2:1481950/93‑1 (MQ=255)
   aTGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTAc                                                                                      >  2:1645923/1‑93 (MQ=255)
   aTGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTAc                                                                                      >  1:1701752/1‑93 (MQ=255)
   aTGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTAc                                                                                      >  1:445068/1‑93 (MQ=255)
    tGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTAcc                                                                                     >  1:4292376/1‑93 (MQ=255)
    tGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTAcc                                                                                     >  2:768049/1‑93 (MQ=255)
     gAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCt                                                                                    >  2:1665529/1‑93 (MQ=255)
      aaTGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTc                                                                                   >  1:3735732/1‑93 (MQ=255)
      aaTGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTc                                                                                   >  1:1380160/1‑93 (MQ=255)
       aTGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCg                                                                                  >  1:3357410/1‑93 (MQ=255)
        tgtgGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGa                                                                                 <  2:3095352/93‑1 (MQ=255)
           ggTCGGTGACCTGCTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGTGCGCGTCATGAAACAGGCGGGGGTCTACCTCGAAcc                                                                              >  2:10370/1‑93 (MQ=255)
            gTCGGCGACCTGTTTGGCGAAGGGAAAATGTTTCTGCCACAGGTGGTCAAATAGGCGCGCGTCATGAAACAGGCGGGGGTCTACCTCGAACCg                                                                             >  2:1301322/1‑93 (MQ=255)
            gTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCg                                                                             >  2:3578665/1‑93 (MQ=255)
              cGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGtt                                                                           <  2:3261651/93‑1 (MQ=255)
              cGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGtt                                                                           >  1:715243/1‑93 (MQ=255)
              cGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCAACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGtt                                                                           >  2:3594156/1‑93 (MQ=255)
                gtgACGTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTa                                                                         <  1:3419474/91‑1 (MQ=255)
                gCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTa                                                                         <  2:3837001/93‑1 (MQ=255)
                gCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTa                                                                         <  1:1459545/93‑1 (MQ=255)
                  gACCGGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTAtt                                                                       <  1:3919116/93‑1 (MQ=255)
                   aCCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTg                                                                      <  2:3140369/93‑1 (MQ=255)
                     cTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGaa                                                                    >  2:1082416/1‑93 (MQ=255)
                     cTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGaa                                                                    <  2:2466206/93‑1 (MQ=255)
                      tGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAg                                                                   >  1:4394077/1‑93 (MQ=255)
                      tGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAg                                                                   >  2:2302554/1‑93 (MQ=255)
                      tGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAg                                                                   >  2:3341998/1‑93 (MQ=255)
                       gTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGc                                                                  >  1:998995/1‑93 (MQ=255)
                       gTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGCCGGTGGTCTACCTCGAACCGTTTATTGAAGc                                                                  >  1:1647626/1‑93 (MQ=255)
                        tttGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGcc                                                                 >  1:2791627/1‑93 (MQ=255)
                        ttgggCGAAGGGAAAATGTCCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGcc                                                                 <  2:722244/90‑1 (MQ=255)
                         ttGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCa                                                                <  1:835115/93‑1 (MQ=255)
                            gcaaaGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCa                                                             <  2:2283059/90‑1 (MQ=255)
                              gaggggAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGTGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCaaa                                                           <  1:4633867/90‑1 (MQ=255)
                              gaggggAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCaaa                                                           <  2:3373131/90‑1 (MQ=255)
                              gAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCaaa                                                           <  2:2436362/93‑1 (MQ=255)
                                aGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAaga                                                         <  1:2616338/93‑1 (MQ=255)
                                aGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAaga                                                         <  2:1621997/93‑1 (MQ=255)
                                aGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAaga                                                         >  2:4345576/1‑93 (MQ=255)
                                aGGGAAAAGGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAaga                                                         <  2:1879385/93‑1 (MQ=255)
                                 gggAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAagag                                                        >  1:291092/1‑93 (MQ=255)
                                 gggAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAagag                                                        >  1:1773826/1‑93 (MQ=255)
                                  ggAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGc                                                       >  1:1873395/1‑93 (MQ=255)
                                  ggAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGc                                                       >  1:1581294/1‑93 (MQ=255)
                                   gAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCa                                                      >  2:1711865/1‑93 (MQ=255)
                                    aaaaTTTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAg                                                     <  1:4552749/93‑1 (MQ=255)
                                    aaaaTGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAg                                                     <  1:2470808/93‑1 (MQ=255)
                                    aaaaTGTCCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAg                                                     <  1:4059074/93‑1 (MQ=255)
                                      aaTGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAggg                                                   >  2:678972/1‑93 (MQ=255)
                                      aaTGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAggg                                                   <  1:878445/93‑1 (MQ=255)
                                       aTGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGc                                                  >  2:1127109/1‑93 (MQ=255)
                                          ttCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCaaa                                               >  2:2555835/1‑93 (MQ=255)
                                             cTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAAcc                                            >  2:2831737/1‑93 (MQ=255)
                                               gCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCaa                                          <  1:1038643/93‑1 (MQ=255)
                                                ccACAGGGGGTCAAATCGGCGCGCGTCATGAAACAGGCGGGGGTCTACCTCGAACCGGTTATTGAAGCCAGCAAAGCGCAGCGCAAAACCAAc                                         >  2:2591209/1‑93 (MQ=255)
                                                 cacaGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACg                                        >  1:921988/1‑93 (MQ=255)
                                                 cacaGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACg                                        >  1:2205813/1‑93 (MQ=255)
                                                  acaGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACgg                                       >  1:1235041/1‑93 (MQ=255)
                                                   caGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGc                                      >  1:3835710/1‑93 (MQ=255)
                                                   caGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGc                                      <  1:3224535/93‑1 (MQ=255)
                                                    aGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCa                                     >  1:2092558/1‑93 (MQ=255)
                                                    aGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCa                                     >  1:2189386/1‑93 (MQ=255)
                                                    aGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCa                                     >  1:1712950/1‑93 (MQ=255)
                                                     ggtggtCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCaa                                    <  1:3296333/93‑1 (MQ=255)
                                                     ggtggtCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCaa                                    <  2:3433035/93‑1 (MQ=255)
                                                      gtggtCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAg                                   >  1:242869/1‑93 (MQ=255)
                                                        ggtCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGAt                                 <  2:2595895/93‑1 (MQ=255)
                                                         gtCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATg                                <  1:293117/93‑1 (MQ=255)
                                                         gtCAAATCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATg                                >  2:1679306/1‑93 (MQ=255)
                                                         gtCAAATCGGCGCGCGTCATGAAACAGGCGGTGGGCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATg                                >  2:408483/1‑93 (MQ=255)
                                                            aaaTCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTg                             <  1:1140712/93‑1 (MQ=255)
                                                             aaTCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGa                            <  2:1795540/93‑1 (MQ=255)
                                                             aaTCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGa                            <  2:1271129/93‑1 (MQ=255)
                                                             aaTCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGa                            >  1:303156/1‑93 (MQ=255)
                                                              aTCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGAt                           <  1:2359351/93‑1 (MQ=255)
                                                              aTCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGAt                           <  1:3828706/93‑1 (MQ=255)
                                                               tCGGCGCGCGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATc                          <  1:3113491/93‑1 (MQ=255)
                                                                  gcgcgcGTCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGcc                       <  2:468096/93‑1 (MQ=255)
                                                                         tCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGa                <  1:3080813/93‑1 (MQ=255)
                                                                         tCATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGa                <  1:950156/93‑1 (MQ=255)
                                                                          cATGAAGCAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGCGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGaa               <  1:634703/93‑1 (MQ=255)
                                                                          cATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGaa               <  2:1671094/93‑1 (MQ=255)
                                                                          cATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGaa               <  2:4323877/93‑1 (MQ=255)
                                                                          cATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGaa               <  1:790750/93‑1 (MQ=255)
                                                                          cATGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGaa               <  1:4481058/93‑1 (MQ=255)
                                                                           aTGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAg              >  1:2020833/1‑93 (MQ=255)
                                                                            tGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAgg             <  1:4996/93‑1 (MQ=255)
                                                                            tGAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAgg             <  1:1827678/93‑1 (MQ=255)
                                                                             gAAACAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGCAggg            >  2:2677622/1‑93 (MQ=255)
                                                                               aaCAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCg          >  2:4218780/1‑93 (MQ=255)
                                                                                aCAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGa         >  1:3555134/1‑93 (MQ=255)
                                                                                 cAGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGag        >  1:1316242/1‑92 (MQ=255)
                                                                                  aGGCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACg       >  1:1090582/1‑93 (MQ=255)
                                                                                   ggCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGt      >  1:3768074/1‑93 (MQ=255)
                                                                                    gCGGTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTc     >  1:1911518/1‑93 (MQ=255)
                                                                                      ggtggtCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCa   <  1:165933/93‑1 (MQ=255)
                                                                                      ggtggtCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCa   <  2:896913/93‑1 (MQ=255)
                                                                                       gTGGTCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCACCGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCAc  <  2:3539731/93‑1 (MQ=255)
                                                                                           |                                                                                        
GGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCACAGGTGGTCAAATCGGCGCGCGTCATGAAACAGGCGGTGGCCTACCTCGAACCGTTTATTGAAGCCAGCAAAGAGCAGGGCAAAACCAACGGCAAGATGGTGATCGCCACCGTGAAGGGCGACGTCCAC  >  NC_000913/4223504‑4223683

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