Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,574,062 A→T 49.0% C273S (TGC→AGC)  asd ← aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,574,0620AT49.0% 1.5 / 297.0 198C273S (TGC→AGC) asdaspartate‑semialdehyde dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (42/59);  new base T (50/47);  total (92/106)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 201 total aligned reads displayed.

GAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCTTCGCG  >  NC_000913/3573963‑3574161
                                                                                                   |                                                                                                   
gAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGCATGCCTGGCTGTGGCa                                                                                                     >  1:2495754/1‑100 (MQ=255)
 aaCGACTTTCGTCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAg                                                                                                    >  1:2176141/1‑100 (MQ=255)
 aaCGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAg                                                                                                    >  1:3760481/1‑100 (MQ=255)
    gACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTgcgc                                                                                                 <  2:349000/100‑1 (MQ=255)
    gACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAgcgc                                                                                                 <  1:2379803/100‑1 (MQ=255)
      cTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCaa                                                                                               >  2:3293838/1‑100 (MQ=255)
       tttCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAAt                                                                                              >  1:513340/1‑100 (MQ=255)
        ttCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATg                                                                                             <  2:4103072/100‑1 (MQ=255)
        ttCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATg                                                                                             <  1:607087/100‑1 (MQ=255)
          cGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGcc                                                                                           >  2:3681270/1‑100 (MQ=255)
           gCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGccc                                                                                          >  2:2462829/1‑100 (MQ=255)
           gCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGccc                                                                                          <  2:406264/100‑1 (MQ=255)
            cccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATCGACACATCTTTTTTCAATTTAATAGGGAATGCCTGGCTGTGGCAGCGCAATGcccc                                                                                         >  1:1287291/1‑100 (MQ=255)
             ccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCg                                                                                        <  1:2676976/100‑1 (MQ=255)
               aCGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGac                                                                                      <  1:954271/100‑1 (MQ=255)
                  gATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACg                                                                                   <  2:3209712/100‑1 (MQ=255)
                    ttGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGca                                                                                 <  1:3293827/100‑1 (MQ=255)
                     tgtgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGcac                                                                                <  1:3222840/100‑1 (MQ=255)
                       tgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcacac                                                                              >  2:3467243/1‑100 (MQ=255)
                           cagccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATaa                                                                          >  2:295899/1‑100 (MQ=255)
                            agccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATaaa                                                                         <  2:1373282/100‑1 (MQ=255)
                             gccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAAc                                                                        <  1:3873888/100‑1 (MQ=255)
                              ccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAAcc                                                                       <  2:2795948/100‑1 (MQ=255)
                               cagcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCa                                                                      >  2:1951809/1‑100 (MQ=255)
                               cagcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCa                                                                      >  1:3417991/1‑100 (MQ=255)
                                agcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTTTTTTTTTTTGTGAATGCCTGGCTGTGGCTGCGATCTGCTTCTACATGCGCACATAAACCAt                                                                     >  1:3519892/1‑100 (MQ=255)
                                agcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCAt                                                                     >  2:2165046/1‑100 (MQ=255)
                                  cagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCt                                                                   <  1:241406/100‑1 (MQ=255)
                                   agTTCTTCCACGGTCGGAATAGACACATCTTTTTTTTTTTTTTTTTTGAATGCCTGGCTGTGGCTGCGCAATGCCCTGACACGCACACATAAACCATCTa                                                                  >  1:1731844/1‑100 (MQ=255)
                                     ttcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTAcc                                                                <  2:650928/100‑1 (MQ=255)
                                     ttcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCACAATGCCCCGACACGCACACATAAACCATCTAcc                                                                <  2:3106468/100‑1 (MQ=255)
                                        ttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGa                                                             <  1:3706659/100‑1 (MQ=255)
                                         tcCACGGTCGGAATAGACACATCTTTTTTTTTTTTTTTTTTAAATGCCTAGCTGTTGCTGCTCAATTCTCTGACACGCACACATAAACCATCTAcc                                                                >  1:1261029/1‑96 (MQ=255)
                                          ccACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAAt                                                           <  1:1017773/100‑1 (MQ=255)
                                           cACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATACCCCGACACGCACACATAAACCATCTACCGGAAtt                                                          <  2:3012504/100‑1 (MQ=255)
                                             cGGGCGGAAGAGACACATCTTTTTTCAATTTAGTAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGaa                                                            <  1:1212150/96‑1 (MQ=255)
                                               gTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACgg                                                      >  2:2554107/1‑100 (MQ=255)
                                               gTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACgg                                                      >  1:3078373/1‑100 (MQ=255)
                                                tCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGa                                                     <  2:1712705/100‑1 (MQ=255)
                                                tCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCATGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACTGAAATTACGGa                                                     >  1:2208833/1‑100 (MQ=255)
                                                   gAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCACAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGa                                                  >  1:3146830/1‑100 (MQ=255)
                                                   gAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGa                                                  <  1:2462823/100‑1 (MQ=255)
                                                     aTAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtg                                                <  2:4234869/100‑1 (MQ=255)
                                                     aTAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtg                                                <  1:1819446/100‑1 (MQ=255)
                                                       aGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtggt                                              <  1:137492/100‑3 (MQ=38)
                                                       aGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtgtg                                              <  2:2437817/100‑1 (MQ=255)
                                                        gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATgtgt                                             <  2:3988442/100‑1 (MQ=255)
                                                        gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCCATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtgtg                                              >  2:2475273/1‑99 (MQ=255)
                                                        gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACAAAAAAACCATCTACCGGAATGACGGAAGATgtgt                                             >  2:2073432/1‑100 (MQ=255)
                                                          cacaTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTg                                           >  2:3608032/1‑100 (MQ=255)
                                                          cacaTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTg                                           <  2:273782/100‑1 (MQ=255)
                                                          cacaTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTg                                           <  1:879538/100‑1 (MQ=255)
                                                            caTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtggtttatg                                         <  2:3060569/100‑8 (MQ=25)
                                                            caTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtggtttat                                          <  1:1989466/99‑7 (MQ=37)
                                                            caTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAg                                         >  1:2303547/1‑100 (MQ=255)
                                                            caTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAg                                         >  2:1787103/1‑100 (MQ=255)
                                                              tCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtggtttatgtg                                       <  2:1079921/100‑10 (MQ=25)
                                                               cTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAgg                                        >  2:1279106/1‑98 (MQ=255)
                                                                 tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                                    >  2:3601299/1‑100 (MQ=255)
                                                                 tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                                    <  2:3875839/100‑1 (MQ=255)
                                                                 tttttCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                                    <  1:1315436/100‑1 (MQ=255)
                                                                   tttCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttg                                  >  2:3721246/1‑100 (MQ=255)
                                                                   tttCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttg                                  >  1:1301317/1‑100 (MQ=255)
                                                                     tCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgtt                                <  2:1508748/100‑1 (MQ=255)
                                                                     tCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgtt                                >  2:394135/1‑100 (MQ=255)
                                                                      cAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgttg                               >  1:2031786/1‑100 (MQ=255)
                                                                      cAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgttg                               >  2:1536570/1‑100 (MQ=255)
                                                                       aaTTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAAGCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTgg                              <  1:4057589/100‑1 (MQ=255)
                                                                       aaTTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTgg                              >  1:2950370/1‑100 (MQ=255)
                                                                        aTTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGt                             >  1:1390532/1‑100 (MQ=255)
                                                                         tttAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGtt                            <  2:576482/100‑1 (MQ=255)
                                                                          ttAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGTTCTTGTTGGttt                           <  2:2815766/100‑1 (MQ=255)
                                                                          ttAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGttt                           <  1:2100232/100‑1 (MQ=255)
                                                                           tgaTAGTGATTGCCTGGCTGTGGCTGCGCAATGCCCCGAAACGCACACATAAACCTTCTACTGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTc                          <  1:658283/98‑1 (MQ=255)
                                                                             aTAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGCGGATCTCGTTGGTTTCCg                        >  2:1248366/1‑100 (MQ=255)
                                                                              tAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGc                       >  1:2116594/1‑100 (MQ=255)
                                                                              tAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGc                       <  2:3603455/100‑1 (MQ=255)
                                                                              tAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACAGAAACCATCTACCGGAATTACGGTAGATGTGTTGAGGATCTTGTTGGTTTCCGc                       <  1:2806824/100‑1 (MQ=255)
                                                                               actGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGcc                      <  2:2202612/98‑1 (MQ=255)
                                                                               aGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGcc                      >  1:4001505/1‑100 (MQ=255)
                                                                               aGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGcc                      >  1:1102872/1‑100 (MQ=255)
                                                                               aGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGcc                      >  2:275980/1‑100 (MQ=255)
                                                                                gTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCt                     >  1:726286/1‑100 (MQ=255)
                                                                                 tGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTg                    >  2:1696819/1‑100 (MQ=255)
                                                                                 gatctGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTg                    <  1:1754368/96‑1 (MQ=255)
                                                                                  gAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGc                   >  2:843282/1‑100 (MQ=255)
                                                                                  gAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGc                   >  1:1152780/1‑100 (MQ=255)
                                                                                     tGCCTGGCTGTGAAAGCGCAATGCCCCGACACGCACACATAGACCCGCTACCGGAATTACGGAAGATGTGTTGAGGATCTTATTGTTTTCCGCCTGccc                 <  1:3817380/99‑1 (MQ=255)
                                                                                      gCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTgagatc               >  2:1754373/1‑95 (MQ=255)
                                                                                       ccTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCttt              <  1:1795509/100‑1 (MQ=255)
                                                                                         tGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTcc            <  2:3760498/100‑1 (MQ=255)
                                                                                          ggCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCa           >  1:2298913/1‑100 (MQ=255)
                                                                                           gCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCAc          <  2:1331333/100‑1 (MQ=255)
                                                                                           gCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGCTCCTGTTGGTTTCCGCCTGCCCTTTCCAc          >  2:1571234/1‑100 (MQ=255)
                                                                                             tgtgGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCActc        <  1:1905403/100‑1 (MQ=255)
                                                                                             tgtgGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCActc        <  1:1528684/100‑1 (MQ=255)
                                                                                                ggCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCTTc     >  2:162770/1‑100 (MQ=255)
                                                                                                ggCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCTTc     >  1:3797804/1‑100 (MQ=255)
                                                                                                ggCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCTTc     >  1:3676655/1‑100 (MQ=255)
                                                                                                  cAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCTTcgc   <  1:4253658/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCTTCGCG  >  NC_000913/3573963‑3574161

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: