Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 401,542 G→T 76.8% E53* (GAG→TAG)  iraP → anti‑RssB factor, RpoS stabilzer during Pi starvation; anti‑adapter protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913401,5420GT76.8% 266.5 / 91.9 139E53* (GAG→TAG) iraPanti‑RssB factor, RpoS stabilzer during Pi starvation; anti‑adapter protein
Reads supporting (aligned to +/- strand):  ref base G (16/16);  new base T (54/52);  total (71/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 141 total aligned reads displayed.

AAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTC  >  NC_000913/401443‑401641
                                                                                                   |                                                                                                   
aaaGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTag                                                                                                     >  1:3334972/1‑100 (MQ=255)
    aaCTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAggg                                                                                                 <  1:458245/100‑1 (MQ=255)
     aCTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAgggg                                                                                                <  2:486513/100‑1 (MQ=255)
      cTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAggggg                                                                                               >  1:1938533/1‑100 (MQ=255)
       tgtgtgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGgc                                                                                              <  1:2799461/100‑1 (MQ=255)
        gtgtgCCCAGGTAGAATCTTTGGAGAGAATCGTCACTGCAATGCGTCGCAAGATGGCGCAAAATGACCAACAGCGGTGGATTGATCAGGTATAGGGGgc                                                                                              <  2:3084454/99‑1 (MQ=255)
         tgtgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGgcgc                                                                                            <  1:2798240/100‑1 (MQ=255)
          gtgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCt                                                                                           >  2:1296744/1‑100 (MQ=255)
          gtgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCt                                                                                           >  1:732532/1‑100 (MQ=255)
            gcgcAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTACGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGt                                                                                         >  2:1990547/1‑100 (MQ=255)
             cgcAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGGAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTa                                                                                        <  2:2027133/100‑1 (MQ=255)
             cgcAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTa                                                                                        <  1:3205266/100‑1 (MQ=255)
              gcAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTAc                                                                                       <  1:1582665/100‑1 (MQ=255)
               cAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGCTCAGGTATAGCGGGCGCTGTACg                                                                                      >  2:1100923/1‑100 (MQ=255)
                 ggTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAg                                                                                    <  1:2524508/100‑1 (MQ=255)
                 ggTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAg                                                                                    >  2:2524492/1‑100 (MQ=255)
                 ggTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGCACGAg                                                                                    <  2:1691310/100‑1 (MQ=255)
                  gTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAgg                                                                                   <  1:929854/100‑1 (MQ=255)
                  gTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAgg                                                                                   <  1:2534051/100‑1 (MQ=255)
                  gTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAgg                                                                                   <  1:510288/100‑1 (MQ=255)
                   tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGt                                                                                  <  1:1200174/100‑1 (MQ=255)
                   tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGt                                                                                  <  1:889070/100‑1 (MQ=255)
                    aGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTa                                                                                 >  1:1195238/1‑100 (MQ=255)
                        ggtttGGAGATTATGGTCACTGCAATGCTTTACAATATGGCGTAAAATGACCAACAGGGGTTGGTAGGTCGGGTACATGGGGAGCTGTACGAGGTAAAGc                                                                             >  2:912798/3‑100 (MQ=255)
                          tttGGATCTTATTGTCACGGCAATGCTTCGCATTGTGGCGCATAATGTCCGACAGCGGTTGATTGATGAGGGATATGAGGCACTGTACGAGGTTAAGc                                                                             <  2:3461718/98‑1 (MQ=255)
                          tttGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGccc                                                                           >  1:3432041/1‑100 (MQ=255)
                               agaTTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGc                                                                      >  1:1791014/1‑100 (MQ=255)
                                naTTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGc                                                                      >  1:2904121/2‑99 (MQ=255)
                                  ttATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAg                                                                   <  1:1363602/100‑1 (MQ=255)
                                      cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCAtt                                                               >  2:1424578/1‑100 (MQ=255)
                                      cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCAtt                                                               >  2:1349712/1‑100 (MQ=255)
                                      cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCAtt                                                               >  1:1334506/1‑100 (MQ=255)
                                       gTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCGGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTc                                                              >  2:2083023/1‑100 (MQ=255)
                                        tCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTcc                                                             <  2:2373148/100‑1 (MQ=255)
                                         cACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCt                                                            <  2:1867772/100‑1 (MQ=255)
                                          aCTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTg                                                           <  1:142272/100‑1 (MQ=255)
                                          aCTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTg                                                           <  1:1146098/100‑1 (MQ=255)
                                           cTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATATGGGGAGCTGTACGAGGTAAAGTCCGATGCCAGCATTCCTga                                                          >  2:3068467/1‑100 (MQ=255)
                                             gCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacg                                                        <  1:2299776/100‑1 (MQ=255)
                                             gCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacg                                                        >  1:2215110/1‑100 (MQ=255)
                                             gCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgac                                                         >  1:2766420/1‑99 (MQ=255)
                                                aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGacgacg                                                     <  2:2656686/100‑1 (MQ=255)
                                                 tGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGAcgacga                                                    <  1:1675373/100‑1 (MQ=255)
                                                  gCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGAt                                                   >  1:2137805/1‑100 (MQ=255)
                                                  gCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCAGATGCCAGCATTCCTGACGACGAt                                                   >  2:42919/1‑100 (MQ=255)
                                                   cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                  >  1:2333060/1‑100 (MQ=255)
                                                   cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                  >  1:218225/1‑100 (MQ=255)
                                                    ttCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGGCAGCATTCCTGAAGACGATAc                                                 >  1:2388566/1‑100 (MQ=255)
                                                     tCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCTATGCCAGCATTCCTGACGACGATACg                                                <  2:2685870/100‑1 (MQ=255)
                                                     tCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACg                                                >  2:3182658/1‑100 (MQ=255)
                                                      cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCACGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACgg                                               >  2:1987700/1‑100 (MQ=255)
                                                       gCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGa                                              >  2:2670810/1‑100 (MQ=255)
                                                        cAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAg                                             >  2:2971558/1‑100 (MQ=255)
                                                         aaTATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgc                                            <  1:2615718/100‑1 (MQ=255)
                                                           tatGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctg                                          <  2:3203006/100‑1 (MQ=255)
                                                           tatGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctg                                          <  2:3369126/100‑1 (MQ=255)
                                                             tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctgct                                        >  1:2233927/1‑100 (MQ=255)
                                                              ngcgcAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGctgctg                                       >  2:1143081/2‑100 (MQ=255)
                                                              ggCGCAAAATGACCAACAGCTGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGctgctg                                       <  2:1308646/100‑1 (MQ=255)
                                                              ggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGctgctg                                       >  1:756571/1‑100 (MQ=255)
                                                              ggCGCAAAATGACCAACAGCGGATGATTGATCAGGTAGAGGGGGCGCTGTACGAGTGAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctgc                                         >  2:2121023/1‑98 (MQ=255)
                                                               gcgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCtgctgc                                      <  2:254852/100‑1 (MQ=255)
                                                                cgcAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTgcg                                     >  1:3051166/1‑100 (MQ=255)
                                                                  cAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGcgcg                                   >  1:1920387/1‑100 (MQ=255)
                                                                  cAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGcgcg                                   <  1:2246554/100‑1 (MQ=255)
                                                                  cAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGcgcg                                   <  1:125654/100‑1 (MQ=255)
                                                                  cAAAATGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGcgcg                                   <  1:125653/100‑1 (MQ=255)
                                                                      aTGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTa                               <  2:3381778/100‑1 (MQ=255)
                                                                      aTGACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTa                               >  1:2642995/1‑100 (MQ=255)
                                                                      aTGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTa                               <  2:2890115/100‑1 (MQ=255)
                                                                        gACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACg                             <  1:1098827/100‑1 (MQ=255)
                                                                        gACCAACAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACg                             <  1:3490303/100‑1 (MQ=255)
                                                                          ccAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTa                           >  2:912530/1‑100 (MQ=255)
                                                                            aaCAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                         <  1:2904775/100‑1 (MQ=255)
                                                                            aaCAGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                         <  1:355128/100‑1 (MQ=255)
                                                                            aaCAGCGGTTGAGTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                         >  1:2393624/1‑100 (MQ=255)
                                                                               aGCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaa                      <  1:1211168/100‑1 (MQ=255)
                                                                                gCGGTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                     <  2:656060/100‑1 (MQ=255)
                                                                                gCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                     <  1:3067383/100‑1 (MQ=255)
                                                                                 cGGTTGAGTGATCAGGTAGAGGGGACGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                     <  2:87540/99‑1 (MQ=255)
                                                                                 cGGGTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGt                    <  2:1428425/100‑1 (MQ=255)
                                                                                  ggTTGATTGATCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGtt                   >  2:2354625/1‑100 (MQ=255)
                                                                                   gTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTa                  <  2:2441870/100‑1 (MQ=255)
                                                                                    ttgattgaTCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTAt                 <  2:1562299/100‑1 (MQ=255)
                                                                                      gattgatCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTAtt                >  2:2992477/1‑99 (MQ=255)
                                                                                      gattgatCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTAtt                >  2:1610199/1‑99 (MQ=255)
                                                                                      gattgatCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTg               >  1:922762/1‑100 (MQ=255)
                                                                                      gattgatCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTg               <  1:2260190/100‑1 (MQ=255)
                                                                                       attgatCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGa              >  1:105531/1‑100 (MQ=255)
                                                                                       attgatCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGa              <  1:3031533/100‑1 (MQ=255)
                                                                                        ttgatCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa             <  2:1791006/100‑1 (MQ=255)
                                                                                         tgatTAGGTATAGGGGGCGTTGTACGAGGTAAAGTTTGATGATAGCATTCTTTACGACGATATGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAg            <  2:884661/100‑1 (MQ=255)
                                                                                         tgatCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAg            <  1:1992523/100‑1 (MQ=255)
                                                                                          gatCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGc           >  1:1774663/1‑100 (MQ=255)
                                                                                           atCAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCa          >  2:28314/1‑100 (MQ=255)
                                                                                           agcAGGTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGTGCGATTACGTAAAGAAGTTATTGAAGCa          <  2:3331546/98‑1 (MQ=255)
                                                                                              aGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATcc       >  2:307723/1‑100 (MQ=255)
                                                                                                gTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTc     >  2:708201/1‑100 (MQ=255)
                                                                                                gTATAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTc     >  2:2068209/1‑100 (MQ=255)
                                                                                                   gagGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTc  <  2:2215096/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTC  >  NC_000913/401443‑401641

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