Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,958 G→T 15.3% intergenic (‑45/+28) ridA ← / ← pyrI enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9580GT15.3% 430.3 / 61.9 182intergenic (‑45/+28)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (76/75);  new base T (12/16);  total (90/92)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-02

Only 100 of 191 total aligned reads displayed.

GTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGG  >  NC_000913/4470859‑4471056
                                                                                                   |                                                                                                  
tccgatctggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAgg                                                                                                    <  2:454925/92‑1 (MQ=255)
 tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAggg                                                                                                   <  2:1655531/100‑1 (MQ=255)
  aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGa                                                                                                  >  1:619205/1‑100 (MQ=255)
   cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                                 <  1:996599/100‑1 (MQ=255)
     tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGcc                                                                                               <  1:1166810/100‑1 (MQ=255)
     tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGcc                                                                                               >  2:2980877/1‑100 (MQ=255)
     tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGc                                                                                                >  1:2599535/1‑99 (MQ=255)
       aGGACCGATAGCTGCCGGTGTATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGTCAg                                                                                             <  1:1743673/100‑1 (MQ=255)
       aGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAg                                                                                             <  1:83749/100‑1 (MQ=255)
        ggACCGATAGCTGCCGGTGCATTTTCCGTCGTGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                            <  1:3109891/100‑1 (MQ=38)
        ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAggagatcgga                                                                                            >  1:454931/1‑92 (MQ=255)
        ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                            >  2:1040512/1‑100 (MQ=255)
         gattgATAGCTGCCGGTGGATTTTTCGTCGCGATAGTTTTGCTCATTGTTTCTCCTTTATTACAGCGGTATAAAGTCTGGCTCATTATAGGGAGCCAGag                                                                                           <  2:446294/96‑2 (MQ=255)
         gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCGAGAc                                                                                           <  2:3031454/100‑1 (MQ=38)
         gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                           >  1:587242/1‑100 (MQ=255)
         gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAACGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                           >  1:3089442/1‑100 (MQ=255)
           ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACTGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                         <  2:1137692/100‑1 (MQ=37)
            tgATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACttt                                                                                        <  1:1415015/99‑1 (MQ=255)
            gtatgccTGCCAGTGTATTTTCCGTCGCGATAGTTGTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                         <  2:3853497/93‑1 (MQ=25)
            cGATAGCTGCCGGTTCATTTTCCGTCGTGCTAGTTTTTCGCATGATTTCTCCTTTATTACAGCGGTTAAAAGTCTGGCTTCCTATAATGAGCCAGACtt                                                                                         <  2:2545616/99‑1 (MQ=21)
               tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                     <  2:828955/100‑1 (MQ=255)
                 gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtttttacc                                                                                   <  2:273618/100‑5 (MQ=25)
                 gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                   <  1:3704659/100‑1 (MQ=255)
                 gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                   >  2:2104134/1‑100 (MQ=255)
                 gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTAGCTCATTATAGGGAGCCAGACTTttatta                                                                                   >  1:782792/1‑100 (MQ=255)
                  cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAc                                                                                  <  1:2786658/100‑1 (MQ=255)
                  cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAc                                                                                  >  1:319392/1‑100 (MQ=255)
                  cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAc                                                                                  <  2:2830451/100‑1 (MQ=255)
                   tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTCAAAAGTCTGGCTAATTATAGGGAGCCAGACTTttatt                                                                                    >  2:380140/1‑97 (MQ=255)
                   tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                                 >  2:3887350/1‑100 (MQ=255)
                    gCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCa                                                                                >  2:2997972/1‑100 (MQ=255)
                     ccGGTGCATTTTCCGTCGCGATAGTTTTACTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCaa                                                                               >  2:1484243/1‑100 (MQ=255)
                      cGGTTCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                              <  1:537995/100‑1 (MQ=255)
                      cGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                              <  2:619198/100‑1 (MQ=255)
                       ggTGCATTTTCCTTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                             <  1:1207011/100‑1 (MQ=255)
                         tGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                           >  2:222321/1‑100 (MQ=255)
                          gCATTTTAAGTCGCGATAGTATGGCTAATGATTGCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGGGCGAGACTTGTATTAACAACCgaa                                                                          >  2:312377/1‑98 (MQ=255)
                           cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaa                                                                         >  2:2620423/1‑100 (MQ=255)
                            aTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                        >  2:1930163/1‑100 (MQ=255)
                            aTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                        <  1:1990270/100‑1 (MQ=255)
                             ttttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatt                                                                       >  1:2492968/1‑100 (MQ=255)
                                tCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                    <  2:2967402/100‑1 (MQ=255)
                                  cGTCTCGATAGTTTTGCTCATGTTTTCTCCTTTATTACATCGGTAAAATGTCTGTCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                   <  1:3274355/99‑1 (MQ=255)
                                   gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCACTTAATTgg                                                                 >  1:220539/1‑100 (MQ=255)
                                     cgcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGc                                                                >  2:1797715/1‑99 (MQ=255)
                                      gcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCa                                                              <  1:609283/100‑1 (MQ=255)
                                       cgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                             <  1:767690/100‑1 (MQ=255)
                                        gATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGc                                                            >  2:2957920/1‑100 (MQ=255)
                                         aTAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGAGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGca                                                           >  1:1017569/1‑100 (MQ=255)
                                          tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGTAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                          >  2:897120/1‑100 (MQ=255)
                                              tttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACcaca                                                      >  1:2431503/1‑100 (MQ=255)
                                                 gCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTa                                                   <  1:1105134/100‑1 (MQ=255)
                                                    cATGATTTCTTCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                <  1:2181461/100‑1 (MQ=255)
                                                    cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                >  1:2296490/1‑100 (MQ=255)
                                                    cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                >  2:2101173/1‑100 (MQ=255)
                                                     aTGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATggg                                               <  2:561870/100‑1 (MQ=255)
                                                      tGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTTATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGa                                              <  1:1530125/100‑1 (MQ=255)
                                                       gATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                             >  2:1701048/1‑100 (MQ=255)
                                                         tttCTCCTTTATTACAGCGGTAAAAAGTCTTGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaaa                                           <  1:26898/100‑1 (MQ=255)
                                                           tctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAct                                         >  1:210224/1‑100 (MQ=255)
                                                            ctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                        >  2:3188356/1‑100 (MQ=255)
                                                             tcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActct                                       >  1:3844003/1‑100 (MQ=255)
                                                               cTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                     >  2:2903685/1‑100 (MQ=255)
                                                                tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                    >  2:3532668/1‑100 (MQ=255)
                                                                 ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                                   <  1:3353690/100‑1 (MQ=255)
                                                                 ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                                   <  1:2980904/100‑1 (MQ=255)
                                                                  tattaCAGCGGTAAAAAGTCTGGCTAATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                  <  2:1541742/100‑1 (MQ=255)
                                                                   attaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                 <  1:1598422/100‑1 (MQ=255)
                                                                    ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                <  2:808463/100‑1 (MQ=255)
                                                                    ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                <  2:2451231/100‑1 (MQ=255)
                                                                     taCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                               <  1:1040523/100‑1 (MQ=255)
                                                                     taCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                               <  1:1825963/100‑1 (MQ=255)
                                                                      aCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                              <  2:2473589/100‑1 (MQ=255)
                                                                      aCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                              <  1:2750052/100‑1 (MQ=255)
                                                                         gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAt                           <  2:1437722/100‑1 (MQ=255)
                                                                           ggTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAttt                         >  2:3241611/1‑100 (MQ=255)
                                                                            gTAAAAAGTCTGGCTCATTATAGGGAGCCAGCCTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAAAGGAAAACTCTTTTTCACAGTATTTg                        >  1:3588834/1‑100 (MQ=255)
                                                                            aataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                        >  1:2867599/4‑100 (MQ=37)
                                                                            aataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                        >  1:382838/4‑100 (MQ=25)
                                                                               aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt                     <  2:1789089/100‑1 (MQ=255)
                                                                                 aaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt                   >  2:2171250/1‑100 (MQ=38)
                                                                                 aaGTCTGGCTCCCTATAATGAGACAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACCTTATGGGAAAACTCTTTTTTACAATATTTGCAttt                   >  2:1424386/1‑100 (MQ=37)
                                                                                  aGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTTGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg                  >  1:3222810/1‑100 (MQ=37)
                                                                                   gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGa                 >  1:2740338/1‑100 (MQ=38)
                                                                                   gTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGa                 >  1:1381892/1‑100 (MQ=255)
                                                                                   gTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGAATTTACATTTGa                 >  1:3623805/1‑100 (MQ=255)
                                                                                    tCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                <  2:124768/100‑1 (MQ=255)
                                                                                       ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc             >  2:1193265/1‑100 (MQ=255)
                                                                                        gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACATTATTTGCATTTGAGcgcg            <  2:1464376/100‑1 (MQ=255)
                                                                                        gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg            >  1:468946/1‑100 (MQ=255)
                                                                                        gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg            <  1:1937222/100‑1 (MQ=255)
                                                                                         cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa           >  2:2957688/1‑100 (MQ=255)
                                                                                          tCATTATATGGAGCTAGATTTTTATTACCAATCGTAATTAATTGGCTAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGat          <  1:682540/100‑1 (MQ=255)
                                                                                          tCATTATAGTTAGCCAGACTTTTATTACGAACTGCAATTAAGTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc             <  1:3142646/97‑1 (MQ=255)
                                                                                           cATTATAGGGAGCCAGACATTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGTAAAACTCTTTTTCACAGTATTTGCATTTTAGCGCGatat        >  2:147491/1‑100 (MQ=255)
                                                                                            aTTATAGCGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGata         <  2:740034/99‑1 (MQ=255)
                                                                                              tataGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCa      <  1:1454909/100‑1 (MQ=255)
                                                                                                taGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAtt    <  2:3354900/100‑1 (MQ=255)
                                                                                                  gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAtt    >  2:959107/1‑98 (MQ=255)
                                                                                                  gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg  >  1:3931061/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGG  >  NC_000913/4470859‑4471056

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