Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,951 A→C 17.1% intergenic (‑38/+35) 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,9510AC17.1% 363.8 / 66.0 156intergenic (‑38/+35)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base A (69/59);  new base C (10/17);  total (79/77)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 161 total aligned reads displayed.

ACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATAT  >  NC_000913/4470852‑4471050
                                                                                                     |                                                                                                   
aCGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCa                                                                                                     <  1:886750/100‑1 (MQ=255)
 cGCCCTGTACGTAAGGACCGATACCTGCCGGTGCACTATCCAAAACGATAGTTTTGCTCATGATTTCTCCTTTAAAAAAAAAAAA‑‑AAAAGTCTGGCTCa                                                                                                     <  2:1006425/99‑1 (MQ=255)
    ccTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATtat                                                                                                 >  1:2660867/1‑100 (MQ=255)
       gTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCTCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAgg                                                                                              <  1:1037937/100‑1 (MQ=255)
       gTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTCTTACCGCGGTA‑‑AAAAGTCTGGCTCATTATAgg                                                                                              >  2:726922/1‑100 (MQ=255)
         aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGa                                                                                            >  2:1494024/1‑100 (MQ=255)
           gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGc                                                                                          <  2:296864/100‑1 (MQ=255)
            tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGcc                                                                                         >  1:1525329/1‑100 (MQ=255)
               ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                      >  1:834127/1‑100 (MQ=255)
               ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                      >  2:1562010/1‑100 (MQ=255)
                gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                     <  2:122247/100‑1 (MQ=38)
                gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                     <  1:3063705/100‑1 (MQ=38)
                gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                     <  2:542332/100‑1 (MQ=38)
                gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                     <  2:205620/100‑1 (MQ=38)
                gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                     <  2:2089529/100‑1 (MQ=38)
                gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                     <  1:3159849/100‑1 (MQ=255)
                 aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTC‑‑AAAAGTCTGGCTCCTTATAGGGCGCCAGa                                                                                      >  1:3039313/1‑98 (MQ=255)
                 aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                    >  1:1347155/1‑100 (MQ=255)
                  ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                   <  1:462579/100‑1 (MQ=38)
                  ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                   >  2:858489/1‑100 (MQ=255)
                  ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                   <  1:2010799/100‑1 (MQ=255)
                    gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                 >  1:1491033/1‑100 (MQ=255)
                    gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                 >  1:640433/1‑100 (MQ=255)
                     aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                 <  2:2889052/99‑1 (MQ=255)
                     aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                <  1:1755429/100‑1 (MQ=255)
                      tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                               >  1:3258462/1‑100 (MQ=255)
                        gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                             >  1:2366602/1‑100 (MQ=255)
                         cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACtttttaccg                                                                            <  1:1675761/100‑6 (MQ=25)
                         cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAc                                                                            <  1:1528628/100‑1 (MQ=255)
                         cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAc                                                                            <  1:478685/100‑1 (MQ=255)
                         cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAc                                                                            >  1:3179375/1‑100 (MQ=255)
                          tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                           >  1:2560156/1‑100 (MQ=255)
                           gCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTGCCa                                                                          >  1:105677/1‑100 (MQ=255)
                             cGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                        >  1:405548/1‑100 (MQ=255)
                              ggTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                       >  2:411603/1‑100 (MQ=255)
                              ggTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                       <  1:1050070/100‑1 (MQ=255)
                              ggTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCAGTATAGGGTGCCAGACTTTTATTACCAAcc                                                                       >  1:1426812/1‑100 (MQ=255)
                               gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAATAAAAGTCTGGCTCCCTATAATGAGCCAGACtttttaccgctgt                                                                        <  1:700050/100‑10 (MQ=17)
                               gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                      <  2:1602399/100‑1 (MQ=255)
                                 gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCa                                                                    >  2:3323309/1‑100 (MQ=255)
                                  cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaa                                                                   >  1:2617804/1‑100 (MQ=255)
                                      ttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTCTAGGGAGCCAGACTTTTATTACCAACCGCaattaa                                                               >  2:2742581/1‑100 (MQ=255)
                                         cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCTGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                             <  2:2690474/99‑1 (MQ=255)
                                         cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTg                                                            >  1:2160833/1‑100 (MQ=255)
                                          gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                           <  1:3581983/100‑1 (MQ=255)
                                          gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAACGCAATTAATTgg                                                           >  2:2390718/1‑100 (MQ=255)
                                             gcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCa                                                        <  2:1872784/100‑1 (MQ=255)
                                              cgATAGTTTTGCTCATGATTTCTCCTTTATTACTGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                       >  2:2385801/1‑100 (MQ=255)
                                                aTAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGca                                                     <  1:1647840/100‑1 (MQ=255)
                                                 tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                    <  1:2708912/100‑1 (MQ=255)
                                                   gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                  >  1:3416401/1‑100 (MQ=255)
                                                    ttttGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcaccac                                                 <  2:2588731/100‑1 (MQ=255)
                                                    ttttGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcaccac                                                 <  2:726684/100‑1 (MQ=255)
                                                    ttttGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcaccac                                                 <  1:3323334/100‑1 (MQ=255)
                                                         cTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                            >  2:446632/1‑100 (MQ=255)
                                                         cTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                            <  2:1459506/100‑1 (MQ=255)
                                                          tCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                           <  1:102250/100‑1 (MQ=255)
                                                           cATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                          >  2:388825/1‑100 (MQ=255)
                                                            aTGATTTCTCCTTTATTACAGCTGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATggg                                         <  2:2461245/100‑1 (MQ=255)
                                                              gATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                       >  2:3138419/1‑100 (MQ=255)
                                                              gATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                       >  2:3138400/1‑100 (MQ=255)
                                                              gATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                       <  2:3515351/100‑1 (MQ=255)
                                                               aTTTCTCCTTTATTACAGCGATA‑‑AAAAGTCTGGCTCATTATAGGGAAACAGACTTTTATTACCAACCGCAATTAATTGGCCAGAAACACATTATagg                                         <  2:1878175/97‑1 (MQ=255)
                                                                tttCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTTGCCAGCACCACATTATGGGaaaa                                     >  2:1200510/1‑100 (MQ=255)
                                                                 tccgatcttattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAc                                    <  2:2639258/93‑1 (MQ=255)
                                                                  tctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAct                                   >  2:85734/1‑100 (MQ=255)
                                                                   ctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                  >  2:1311088/1‑100 (MQ=255)
                                                                    tcCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActct                                 >  1:2041774/1‑100 (MQ=255)
                                                                     ccTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                <  2:827823/100‑1 (MQ=255)
                                                                     ccTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                >  2:530241/1‑100 (MQ=255)
                                                                      cTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                               <  1:2894660/100‑1 (MQ=255)
                                                                         tattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                            <  2:3572828/100‑1 (MQ=255)
                                                                         tattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                            >  2:1785390/1‑100 (MQ=255)
                                                                          attaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                           >  1:2270080/1‑100 (MQ=255)
                                                                           ttaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                          <  2:1026852/100‑1 (MQ=255)
                                                                             aCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                        <  2:36438/100‑1 (MQ=255)
                                                                                gCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAt                     >  2:1571603/1‑100 (MQ=255)
                                                                                   gTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                  >  1:1420844/1‑100 (MQ=255)
                                                                                     aataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                  >  1:2184131/4‑100 (MQ=25)
                                                                                     a‑‑aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCa                >  2:283765/1‑100 (MQ=255)
                                                                                     a‑‑aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCa                >  1:2255771/1‑100 (MQ=255)
                                                                                      ttaaaaGTCTGGCTCTTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCTCCACATTATGGGAAAACTCTTTTTCACAGTATTTg                  <  1:926863/97‑1 (MQ=255)
                                                                                      ataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                    >  2:2410328/3‑97 (MQ=37)
                                                                                        aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt               <  2:893719/100‑1 (MQ=255)
                                                                                        aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt               >  1:1219681/1‑100 (MQ=255)
                                                                                         aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTTTCCAGCACCACATTATTGGAAAACTCTGTGTCACAGTATTTGCAtt              >  2:300138/1‑100 (MQ=25)
                                                                                         aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt              >  2:180004/1‑100 (MQ=37)
                                                                                         aaaGTCTGGCTCACTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACAACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt              >  1:346527/1‑100 (MQ=38)
                                                                                          aaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt             >  1:1985052/1‑100 (MQ=38)
                                                                                           aGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg            <  1:388827/100‑1 (MQ=255)
                                                                                           aGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg            <  1:2319317/100‑1 (MQ=255)
                                                                                            gTCTTGCTCCCTATAGTGAGCCAGACTTTTATTACCAACCTCAATTAATTTGCCAGCTCCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg            <  1:886390/99‑1 (MQ=255)
                                                                                            gTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAACACAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTGGCATTTg            <  2:1780175/99‑1 (MQ=255)
                                                                                             tCTGGCTCATTATAGGGAGCCAGAGTTTTATTACCAACCGCAATTAATGGGCCAGCACCACATTAGGGGAAAACTCGTTTTCACAGTATTTGCATTTGAg          <  2:2914120/100‑1 (MQ=255)
                                                                                              cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc         >  1:2899157/1‑100 (MQ=255)
                                                                                                 gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg      <  1:1079838/100‑1 (MQ=255)
                                                                                                 gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAATCTTTTTCACAGTATTTGCATTTGAGcgcg      >  2:3207069/1‑100 (MQ=255)
                                                                                                   tCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGat    <  1:1785398/100‑1 (MQ=255)
                                                                                                    cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGat    <  2:2842212/99‑1 (MQ=255)
                                                                                                     aTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGTCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat  <  1:1578578/100‑1 (MQ=255)
                                                                                                     |                                                                                                   
ACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATAT  >  NC_000913/4470852‑4471050

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: