Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,496,295 A→G V277V (GTA→GTG cysG → uroporphyrinogen III methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,496,2950AG93.7% 969.0 / 36.8 351V277V (GTA→GTGcysGuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base A (14/8);  new base G (161/167);  total (175/176)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 368 total aligned reads displayed.

CGACCGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTT  >  NC_000913/3496205‑3496385
                                                                                          |                                                                                          
cGACCGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTAcc                                                                                          >  1:4546028/1‑93 (MQ=255)
 gACCGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGccc                                                                                         >  1:4848172/1‑93 (MQ=255)
 gACCGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGccc                                                                                         >  1:3043281/1‑93 (MQ=255)
 gACCGTCTGGTTTCTCACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGccc                                                                                         >  2:1024547/1‑93 (MQ=255)
  aCCGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGcccc                                                                                        >  2:7728550/1‑93 (MQ=255)
  aCCGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGcccc                                                                                        <  2:565633/93‑1 (MQ=255)
  aCCGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTAcccc                                                                                        <  2:2032299/93‑1 (MQ=255)
   ccGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCa                                                                                       >  2:3664792/1‑93 (MQ=255)
    cGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAg                                                                                      >  2:5433824/1‑93 (MQ=255)
    cGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAg                                                                                      >  1:2267954/1‑93 (MQ=255)
     gTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAgg                                                                                     >  2:6132176/1‑93 (MQ=255)
     gTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAgg                                                                                     <  2:1919058/93‑1 (MQ=255)
      tCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGa                                                                                    >  1:7852104/1‑93 (MQ=255)
       cTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGaa                                                                                   <  2:4636860/93‑1 (MQ=255)
       cTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGaa                                                                                   <  2:3019875/93‑1 (MQ=255)
       cTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGaa                                                                                   >  1:8722686/1‑93 (MQ=255)
        tGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAag                                                                                  <  2:9547538/93‑1 (MQ=255)
         gttttCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGGGTGCCCCAGGAaga                                                                                 <  2:7249580/91‑1 (MQ=255)
          gTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAagag                                                                                <  2:8547800/93‑1 (MQ=255)
           tttCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAgaga                                                                               <  2:9081866/93‑1 (MQ=255)
             tCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGAtt                                                                             >  1:8251063/1‑93 (MQ=255)
              cTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTa                                                                            <  1:4774380/93‑1 (MQ=255)
               tGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTaa                                                                           <  2:2687841/93‑1 (MQ=255)
               tGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTaa                                                                           >  2:8162513/1‑93 (MQ=255)
               tGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTaa                                                                           <  2:6810429/93‑1 (MQ=255)
                taCGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAAc                                                                          <  1:2532137/92‑1 (MQ=255)
                gACGATATTATGAATCTGGTACGCCGCGATGCGGCCCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAAc                                                                          >  2:8643409/1‑93 (MQ=255)
                 aCGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAAcc                                                                         <  2:4976529/93‑1 (MQ=255)
                 aCGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAAcc                                                                         >  1:7578277/1‑93 (MQ=255)
                 aCGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAAcc                                                                         <  2:3853958/93‑1 (MQ=255)
                  cGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCa                                                                        <  1:7815018/93‑1 (MQ=255)
                   gATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAg                                                                       <  1:5939807/93‑1 (MQ=255)
                   gATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAg                                                                       <  1:648085/93‑1 (MQ=255)
                   catatTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAg                                                                       <  2:8396199/92‑1 (MQ=255)
                    atatTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGa                                                                      <  1:1454506/93‑1 (MQ=255)
                    atatTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGa                                                                      >  2:4027244/1‑93 (MQ=255)
                    atatTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGa                                                                      >  1:2652881/1‑93 (MQ=255)
                    atatTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGa                                                                      >  1:4417224/1‑93 (MQ=255)
                    atatTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGa                                                                      <  2:2491068/93‑1 (MQ=255)
                       ttatGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATcc                                                                   <  2:7047884/93‑1 (MQ=255)
                        tatGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCct                                                                  <  1:2554473/93‑1 (MQ=255)
                         atGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCctg                                                                 <  1:1897106/93‑1 (MQ=255)
                           gAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCctgct                                                               >  2:6223900/1‑93 (MQ=255)
                             aTCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCtgctgc                                                             <  1:7315473/93‑1 (MQ=255)
                             aTCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCtgctgc                                                             <  1:4261235/93‑1 (MQ=255)
                              tCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCg                                                            >  2:6748108/1‑93 (MQ=255)
                              tCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCg                                                            <  2:5379663/93‑1 (MQ=255)
                               cTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCgg                                                           <  2:4298944/93‑1 (MQ=255)
                                tGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCggg                                                          >  1:353371/1‑93 (MQ=255)
                                tGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCggg                                                          >  1:1333687/1‑93 (MQ=255)
                                 ggTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGa                                                         >  1:6056773/1‑93 (MQ=255)
                                 ggTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGa                                                         <  2:2279815/93‑1 (MQ=255)
                                  gTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGaa                                                        <  1:7000484/93‑1 (MQ=255)
                                    aCGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAgc                                                      >  1:2970588/1‑93 (MQ=255)
                                       ccgcGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCa                                                   >  1:5522299/1‑93 (MQ=255)
                                         gcgATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCaaa                                                 >  2:3571617/1‑93 (MQ=255)
                                            aTGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAg                                              >  2:3434768/1‑93 (MQ=255)
                                            aTGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAg                                              >  1:3312688/1‑93 (MQ=255)
                                             tGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAgg                                             >  1:9870896/1‑93 (MQ=255)
                                              gCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGc                                            <  1:4153686/93‑1 (MQ=255)
                                                ggACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCaa                                          <  2:936401/93‑1 (MQ=255)
                                                 gACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCaaa                                         >  2:644043/1‑93 (MQ=255)
                                                 gACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCaaa                                         <  1:6327500/93‑1 (MQ=255)
                                                  aCCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAAc                                        >  1:102992/1‑93 (MQ=255)
                                                    cGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAATAGGTAAAcgc                                      >  2:8941224/1‑93 (MQ=255)
                                                    cGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAAcgc                                      >  2:7783836/1‑93 (MQ=255)
                                                    cGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAAcgc                                      >  1:1562360/1‑93 (MQ=255)
                                                      tgtTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCgt                                    >  1:3795288/1‑93 (MQ=255)
                                                      tgtTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCgt                                    >  2:4695589/1‑93 (MQ=255)
                                                      tgtTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCgt                                    <  2:8708560/93‑1 (MQ=255)
                                                      tgtTTTCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCgt                                    <  1:5358239/93‑1 (MQ=255)
                                                        ttttCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCgtgg                                  >  2:1307720/1‑93 (MQ=255)
                                                        ttttCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCgtgg                                  >  1:4645191/1‑93 (MQ=255)
                                                          ttCGTCGGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCgtggtg                                <  2:3186009/93‑1 (MQ=255)
                                                           tcgtcgGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGc                               <  1:9656300/93‑1 (MQ=255)
                                                              tcgGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGc                            <  1:2049406/93‑1 (MQ=255)
                                                              tcgGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGc                            <  1:8432111/93‑1 (MQ=255)
                                                              gcgGCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGc                            <  1:6800302/92‑1 (MQ=255)
                                                                 gCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGa                         >  1:615652/1‑93 (MQ=255)
                                                                 gCAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGa                         >  2:180592/1‑93 (MQ=255)
                                                                  cAAACGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGaa                        >  2:261593/1‑93 (MQ=255)
                                                                   aaaCGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGaaa                       <  2:4387953/93‑1 (MQ=255)
                                                                   aaaCGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGaaa                       >  2:7365230/1‑93 (MQ=255)
                                                                    aaCGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAg                      <  2:3827013/93‑1 (MQ=255)
                                                                    aaCGCGCGGGATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAg                      >  1:2763354/1‑93 (MQ=255)
                                                                           gggATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGGGGCGa               >  2:7266056/1‑93 (MQ=255)
                                                                           gggATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGCGCGGCTGAAAGGTGGCGa               <  1:5671638/93‑1 (MQ=255)
                                                                             gATACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGTTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATc             <  2:9888727/93‑1 (MQ=255)
                                                                              aTACCACTGCGTGCCCCAGGAAGAGATTACCCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATcc            <  1:9161207/93‑1 (MQ=255)
                                                                               tACCACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCg           >  2:9321211/1‑93 (MQ=255)
                                                                                 ccACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGtt         >  1:8090247/1‑93 (MQ=255)
                                                                                 ccACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGtt         >  1:10006545/1‑93 (MQ=255)
                                                                                 ccACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGGGCGGCTGAAAGGTGGCGATCCGtt         >  1:8971238/1‑93 (MQ=255)
                                                                                  cattGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGttt        <  2:1477155/90‑1 (MQ=255)
                                                                                  cACTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGttt        >  2:2860720/1‑93 (MQ=255)
                                                                                   aCTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTa       >  2:4091431/1‑93 (MQ=255)
                                                                                    cTGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTAt      <  2:662566/93‑1 (MQ=255)
                                                                                     tGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTagt     >  1:9436881/1‑91 (MQ=255)
                                                                                     tGCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTAtt     >  2:5884881/1‑93 (MQ=255)
                                                                                      gCGTGCCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTAttt    <  1:5495763/93‑1 (MQ=255)
                                                                                          |                                                                                          
CGACCGTCTGGTTTCTGACGATATTATGAATCTGGTACGCCGCGATGCGGACCGTGTTTTCGTCGGCAAACGCGCGGGATACCACTGCGTACCCCAGGAAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTT  >  NC_000913/3496205‑3496385

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