Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,496,430 T→C I322I (ATT→ATC cysG → uroporphyrinogen III methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,496,4300TC99.7% 963.3 / ‑5.6 305I322I (ATT→ATCcysGuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (131/173);  total (131/174)
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.97e-01

Only 100 of 323 total aligned reads displayed.

AGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAG  >  NC_000913/3496340‑3496519
                                                                                          |                                                                                         
 ggCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCg                                                                                        >  2:6238104/1‑93 (MQ=255)
    aaaCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTtc                                                                                     <  2:5853318/93‑1 (MQ=255)
     aaCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTtct                                                                                    <  1:1790713/93‑1 (MQ=255)
     aaCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTtct                                                                                    >  1:4441734/1‑93 (MQ=255)
      aCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTtctc                                                                                   >  2:4071495/1‑93 (MQ=255)
      aCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTtctc                                                                                   >  1:6045444/1‑93 (MQ=255)
       cgcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCg                                                                                  <  2:4270785/93‑1 (MQ=255)
        gcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCgg                                                                                 >  1:1440300/1‑93 (MQ=255)
        gcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCgg                                                                                 <  2:1210385/93‑1 (MQ=255)
             gtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGtt                                                                            <  2:7112138/93‑1 (MQ=255)
             gtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGtt                                                                            >  2:1889143/1‑93 (MQ=255)
             gtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGtt                                                                            >  1:5457568/1‑93 (MQ=255)
              tgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTc                                                                           <  2:6738341/93‑1 (MQ=255)
                 ggCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCgg                                                                        <  2:609549/93‑1 (MQ=255)
                  gCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCggg                                                                       <  1:5701826/93‑1 (MQ=255)
                    tGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCgggca                                                                     <  2:3240204/93‑3 (MQ=255)
                    tGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCgggca                                                                     >  2:2772956/1‑91 (MQ=255)
                     gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAt                                                                    >  2:1864919/1‑93 (MQ=255)
                     gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAt                                                                    >  1:3995696/1‑93 (MQ=255)
                     gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAt                                                                    >  1:2210032/1‑93 (MQ=255)
                     gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAt                                                                    >  1:4965171/1‑93 (MQ=255)
                      aaaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAtt                                                                   <  2:1770183/93‑1 (MQ=255)
                      aaaGGTGGCGATCAGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAtt                                                                   >  1:836549/1‑93 (MQ=255)
                       aaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTa                                                                  >  1:7765278/1‑93 (MQ=255)
                       aaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTa                                                                  >  2:6656031/1‑93 (MQ=255)
                        aGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAc                                                                 >  2:1099000/1‑93 (MQ=255)
                        aGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAc                                                                 >  2:4813761/1‑93 (MQ=255)
                         ggTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAcc                                                                <  2:135444/93‑1 (MQ=255)
                          gTGGAGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCg                                                               >  1:7836703/1‑93 (MQ=255)
                            ggCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGAATCCCGTTCTCGGTGGTTCCGGGCATTACCGCa                                                             >  1:3491079/1‑93 (MQ=255)
                              cGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGc                                                           <  2:4182713/93‑1 (MQ=255)
                               gATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCt                                                          >  2:5514193/1‑93 (MQ=255)
                                 tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                        >  1:2802279/1‑93 (MQ=255)
                                 tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                        <  2:4887677/93‑1 (MQ=255)
                                 tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                        >  2:2261087/1‑93 (MQ=255)
                                   cGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTg                                                      <  2:3485817/93‑1 (MQ=255)
                                    gTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTgg                                                     <  2:7070675/93‑1 (MQ=255)
                                     tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGt                                                    <  1:6734428/93‑1 (MQ=255)
                                     tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGt                                                    <  1:7253200/93‑1 (MQ=255)
                                     tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGt                                                    <  1:6196934/93‑1 (MQ=255)
                                     tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGt                                                    >  2:5114797/1‑93 (MQ=255)
                                     tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGt                                                    >  1:4396434/1‑93 (MQ=255)
                                      ttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGtt                                                   >  2:5695852/1‑93 (MQ=255)
                                      ttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGtt                                                   <  1:3622699/93‑1 (MQ=255)
                                      ttATTTTTGGCCGAGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCGTCGGGtt                                                   >  2:1486234/1‑93 (MQ=255)
                                        aTTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGc                                                 <  2:4379087/93‑1 (MQ=255)
                                        aTTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGc                                                 <  2:6978913/93‑1 (MQ=255)
                                        aTTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGc                                                 <  1:1322915/93‑1 (MQ=255)
                                          ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctc                                               >  1:7549166/1‑93 (MQ=255)
                                          ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctc                                               <  1:1599765/93‑1 (MQ=255)
                                              ggCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGcc                                           <  1:1129274/93‑1 (MQ=255)
                                                 cGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAt                                        <  2:4484362/93‑1 (MQ=255)
                                                  gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAtt                                       >  2:5339349/1‑93 (MQ=255)
                                                  gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAtt                                       >  2:6324226/1‑93 (MQ=255)
                                                  gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAtt                                       <  1:3548998/93‑1 (MQ=255)
                                                   tggtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTc                                      >  1:1006254/1‑93 (MQ=255)
                                                   tggtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTc                                      <  1:6688576/93‑1 (MQ=255)
                                                    ggtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCg                                     <  2:4660868/93‑1 (MQ=255)
                                                     gtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCgg                                    <  2:1985701/93‑1 (MQ=255)
                                                      tggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCggg                                   >  2:2585850/1‑93 (MQ=255)
                                                      tggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCggg                                   >  1:5863248/1‑93 (MQ=255)
                                                         cGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAt                                <  1:5305275/93‑1 (MQ=255)
                                                         cGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAt                                <  2:2154897/93‑1 (MQ=255)
                                                         cGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAt                                <  1:6942000/93‑1 (MQ=255)
                                                          gAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAtt                               <  1:6346377/93‑1 (MQ=255)
                                                           aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                              >  1:3256571/1‑93 (MQ=255)
                                                            agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                             >  2:4719484/1‑93 (MQ=255)
                                                            agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                             >  2:860127/1‑93 (MQ=255)
                                                            agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                             <  2:3625018/93‑1 (MQ=255)
                                                            agagCTGGAAACACTGTGCAACGCCGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                             <  1:6278588/93‑1 (MQ=255)
                                                              agCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCAc                           <  1:6826072/93‑1 (MQ=255)
                                                                  ggAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCAc                       <  2:7905279/93‑1 (MQ=255)
                                                                  ggAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCAc                       >  1:4768415/1‑93 (MQ=255)
                                                                   gAAACACTGTGCAACGCGGGTATCTCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACg                      >  2:3348686/1‑93 (MQ=255)
                                                                   gAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACg                      <  2:5018670/93‑1 (MQ=255)
                                                                    aaaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGc                     <  2:1449414/93‑1 (MQ=255)
                                                                     aaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCa                    >  1:2444449/1‑93 (MQ=255)
                                                                     aaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCa                    >  1:6439797/1‑93 (MQ=255)
                                                                      acacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCAt                   <  1:695550/93‑1 (MQ=255)
                                                                       cacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATc                  >  2:2774301/1‑93 (MQ=255)
                                                                       cacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATc                  <  2:5764207/93‑1 (MQ=255)
                                                                       cacTGTGCAACGCAGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATc                  <  1:6993060/93‑1 (MQ=255)
                                                                        acTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcg                 <  2:7651071/93‑1 (MQ=255)
                                                                              cAACGCTGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa           <  1:197015/93‑1 (MQ=255)
                                                                              cAACGCGGGTATCCGGTTCTCGGTGTGTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa           <  2:4216896/93‑1 (MQ=255)
                                                                              cAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa           <  2:3163916/93‑1 (MQ=255)
                                                                              cAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa           <  2:3839649/93‑1 (MQ=255)
                                                                               aaCGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTAt          <  2:1670303/93‑1 (MQ=255)
                                                                                aCGCCGGTATCCCGGTCTCGGTGGGTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGGATTCCACTCACGCAGCGCGATTATg         <  2:4516952/93‑1 (MQ=255)
                                                                                 cgcgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGc        <  2:4627722/93‑1 (MQ=255)
                                                                                   cgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGccc      <  2:2409757/93‑1 (MQ=255)
                                                                                    gggTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCa     <  1:427924/93‑1 (MQ=255)
                                                                                    gggTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCa     <  2:6350319/93‑1 (MQ=255)
                                                                                      gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCaga   >  1:4551304/1‑93 (MQ=255)
                                                                                      gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCaga   >  1:5469130/1‑93 (MQ=255)
                                                                                      gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTACTCCACTCACGCATCGCGATTATGCCCaga   <  1:6508614/93‑1 (MQ=255)
                                                                                      gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGATCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCaga   >  2:7575759/1‑93 (MQ=255)
                                                                                       tATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCagag  <  1:4700695/93‑1 (MQ=255)
                                                                                       tATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCagag  <  2:2607144/93‑1 (MQ=255)
                                                                                       tATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCagag  <  2:7336034/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAG  >  NC_000913/3496340‑3496519

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