Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,957 G→A 16.2% intergenic (‑44/+29) 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,9570GA16.2% 383.0 / 55.3 160intergenic (‑44/+29)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (70/64);  new base A (10/16);  total (80/80)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

Only 100 of 168 total aligned reads displayed.

TGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGG  >  NC_000913/4470858‑4471056
                                                                                                     |                                                                                                   
 gTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTCTTACCGCGGTA‑‑AAAAGTCTGGCTCATTATAgg                                                                                                    >  2:726922/1‑100 (MQ=255)
   aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGa                                                                                                  >  2:1494024/1‑100 (MQ=255)
    cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAg                                                                                                 >  1:69225/1‑100 (MQ=255)
    cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAg                                                                                                 >  2:423152/1‑100 (MQ=255)
    cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTGCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAg                                                                                                 >  2:2333247/1‑100 (MQ=255)
      tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGcc                                                                                               >  1:1525329/1‑100 (MQ=255)
         ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                            >  2:1562010/1‑100 (MQ=255)
          gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  2:1796317/100‑1 (MQ=38)
          gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  2:2629747/100‑1 (MQ=38)
          gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  1:2266665/100‑1 (MQ=38)
          gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  2:122247/100‑1 (MQ=38)
          gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  2:1850630/100‑1 (MQ=38)
          gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  2:2089529/100‑1 (MQ=38)
          gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  1:1123630/100‑1 (MQ=38)
          gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                           <  1:3063705/100‑1 (MQ=38)
           aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                          <  1:1046231/100‑1 (MQ=38)
           aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                          <  2:2860909/100‑1 (MQ=38)
           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:893726/1‑100 (MQ=255)
            ccGATAGCTGCCGGTGCATTTTCCGTCGAGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                         <  2:701284/100‑1 (MQ=255)
              gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                       >  1:1491033/1‑100 (MQ=255)
               aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                       >  1:1722835/1‑99 (MQ=255)
               aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                      <  1:1755429/100‑1 (MQ=255)
                tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                     >  1:3258462/1‑100 (MQ=255)
                    tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                                 >  1:2560156/1‑100 (MQ=255)
                     gCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTGCCa                                                                                >  1:105677/1‑100 (MQ=255)
                       cGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                              <  2:2982082/100‑1 (MQ=255)
                       cGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                              <  2:843112/100‑1 (MQ=255)
                         gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAATAAAAGTCTGGCTCCCTATAATGAGCCAGACtttttaccgctgt                                                                              <  1:700050/100‑10 (MQ=17)
                         gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                            <  2:1602399/100‑1 (MQ=255)
                          tacATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                           <  1:1912124/98‑1 (MQ=255)
                           gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCa                                                                          >  2:3323309/1‑100 (MQ=255)
                            cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaa                                                                         >  1:2617804/1‑100 (MQ=255)
                               tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatta                                                                      >  1:1377026/1‑100 (MQ=255)
                                ttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTCTAGGGAGCCAGACTTTTATTACCAACCGCaattaa                                                                     >  2:2742581/1‑100 (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)
                                        cgATAGTTTTGCTCATGATTTCTCCTTTATTACTGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                             >  2:2385801/1‑100 (MQ=255)
                                           tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                          <  1:2708912/100‑1 (MQ=255)
                                             gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                        >  1:3416401/1‑100 (MQ=255)
                                              ttttGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcaccac                                                       <  1:3323334/100‑1 (MQ=255)
                                                ttGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                     >  1:839693/1‑100 (MQ=255)
                                                  gCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTa                                                   >  1:408131/1‑100 (MQ=255)
                                                   cTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                                  <  2:1459506/100‑1 (MQ=255)
                                                   cTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                                  >  2:446632/1‑100 (MQ=255)
                                                    tCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                 <  1:102250/100‑1 (MQ=255)
                                                    tCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                 <  1:1060948/100‑1 (MQ=255)
                                                     cATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                >  2:388825/1‑100 (MQ=255)
                                                      aTGATTTCTCCTTTATTACAGCTGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATggg                                               <  2:2461245/100‑1 (MQ=255)
                                                       tGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGa                                              >  1:2056863/1‑100 (MQ=255)
                                                        gATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                             >  2:3138400/1‑100 (MQ=255)
                                                        gATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                             <  2:3515351/100‑1 (MQ=255)
                                                        gATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                             <  1:1562018/100‑1 (MQ=255)
                                                         aTTTCTCCTTTATTACAGCGATA‑‑AAAAGTCTGGCTCATTATAGGGAAACAGACTTTTATTACCAACCGCAATTAATTGGCCAGAAACACATTATagg                                               <  2:1878175/97‑1 (MQ=255)
                                                          tttCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTTGCCAGCACCACATTATGGGaaaa                                           >  2:1200510/1‑100 (MQ=255)
                                                           ttCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAc                                          <  1:1650949/100‑1 (MQ=255)
                                                             ctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                        >  2:1311088/1‑100 (MQ=255)
                                                               ccTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGCGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActct                                       >  2:2639594/1‑99 (MQ=255)
                                                               ccTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTTGCCAGCACCACATTATGGGAAAACTCtt                                      >  1:1442089/1‑100 (MQ=255)
                                                               ccTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                      >  1:3576712/1‑100 (MQ=255)
                                                               ccTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                      >  2:530241/1‑100 (MQ=255)
                                                               ccTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                      <  2:827823/100‑1 (MQ=255)
                                                                cTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                     <  1:2894660/100‑1 (MQ=255)
                                                                   tattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAcagatcgga                                  >  1:2639276/1‑92 (MQ=255)
                                                                   tattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                  >  2:1785390/1‑100 (MQ=255)
                                                                   tattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                  >  1:2012023/1‑100 (MQ=255)
                                                                   tattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                  <  2:3572828/100‑1 (MQ=255)
                                                                    attaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                 >  1:2331166/1‑100 (MQ=255)
                                                                     ttaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                <  2:1026852/100‑1 (MQ=255)
                                                                       aCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                              <  2:36438/100‑1 (MQ=255)
                                                                          gCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGTCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTa                            <  2:1941524/99‑1 (MQ=255)
                                                                             gtaataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                          >  1:2981655/6‑100 (MQ=25)
                                                                             gTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                        >  1:1420844/1‑100 (MQ=255)
                                                                             gTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                        <  1:519954/100‑1 (MQ=255)
                                                                               a‑‑aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCa                      >  2:283765/1‑100 (MQ=255)
                                                                                 taaaaGTCTGGCTCCCTATAATGAGCCATGCTTTTATTAAGAACCGCAATTAATTGGTCAGGAACACATTATGGAAAAACACTTTTTTATAGTATTTgg                       >  2:2899543/2‑98 (MQ=16)
                                                                                  aaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt                     <  2:893719/100‑1 (MQ=255)
                                                                                   aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt                    >  1:981776/1‑100 (MQ=37)
                                                                                   aaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt                    >  2:3154951/1‑100 (MQ=255)
                                                                                    aaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt                   >  1:1985052/1‑100 (MQ=38)
                                                                                    aaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt                   >  1:224420/1‑100 (MQ=38)
                                                                                     aGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg                  <  1:2319317/100‑1 (MQ=255)
                                                                                     aGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg                  <  1:388827/100‑1 (MQ=255)
                                                                                     aGTCTGGCTCATGATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt                    <  2:3033535/98‑1 (MQ=255)
                                                                                       tCTGGCTCATTATAGGGAGCCAGAGTTTTATTACCAACCGCAATTAATGGGCCAGCACCACATTAGGGGAAAACTCGTTTTCACAGTATTTGCATTTGAg                <  2:2914120/100‑1 (MQ=255)
                                                                                        cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc               >  1:3173378/1‑100 (MQ=255)
                                                                                        cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc               >  1:2899157/1‑100 (MQ=255)
                                                                                          ngCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc             >  1:3044509/2‑100 (MQ=255)
                                                                                           gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATGATGGGAAACCGCTTTTTCACAGTATTTGCATTTAAgcgc             >  2:1578505/1‑99 (MQ=255)
                                                                                              cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGata         <  2:2432314/100‑1 (MQ=255)
                                                                                               aTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGTCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat        <  1:1578578/100‑1 (MQ=255)
                                                                                                ttATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc       >  2:3230063/1‑100 (MQ=255)
                                                                                                 tataGGGAGCCAGACTTTTATTACCAACGCCAATTAATTGGCCAGCACAAAGTTATGGAATAACTCTCTTTCAACGTATTTGCATTTCACCGCGata         >  2:1771620/1‑97 (MQ=255)
                                                                                                   taGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAtt    >  1:1007539/1‑100 (MQ=255)
                                                                                                   taGGGAGCCAGACTTTTATGACAAAACGAAATTAATTGGCCAGCACCACATTATGGGAAAAGTCTTTTGCATAGTATTTGCATTTGAGCGGCATAGCAtt    <  2:543967/100‑1 (MQ=255)
                                                                                                    aGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTg   <  1:933488/100‑1 (MQ=255)
                                                                                                     gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg  <  1:1998176/100‑1 (MQ=255)
                                                                                                     |                                                                                                   
TGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGG  >  NC_000913/4470858‑4471056

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: