Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,615 C→G 16.6% G145G (GGG→GGC) ‡ yfcO ← DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,6150CG16.6% 407.5 / 65.0 156G145G (GGG→GGCyfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (56/73);  new base G (13/13);  total (69/87)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 164 total aligned reads displayed.

CAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACC  >  NC_000913/2449516‑2449714
                                                                                                    |                                                                                                   
cAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGc                                                                                                     <  2:487544/100‑1 (MQ=255)
 aGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTTGTGACTATCAGGcc                                                                                                    >  2:2711768/1‑100 (MQ=255)
  gggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGccc                                                                                                   >  1:1010479/1‑100 (MQ=255)
   ggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGcccc                                                                                                  <  1:2582980/100‑1 (MQ=255)
    gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCt                                                                                                 >  2:3621723/1‑100 (MQ=255)
     tGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTg                                                                                                <  2:3095475/100‑1 (MQ=255)
       aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGAt                                                                                              >  1:2264761/1‑100 (MQ=255)
       aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGAt                                                                                              <  2:3492226/100‑1 (MQ=255)
        gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATa                                                                                             >  2:789880/1‑100 (MQ=255)
           gTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGt                                                                                          <  1:3048100/100‑1 (MQ=255)
             tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCa                                                                                        <  1:1732054/100‑1 (MQ=255)
             tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGGGACGTATCAGGGGCCTGATagtca                                                                                         <  2:1890675/100‑5 (MQ=255)
              ttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCAc                                                                                       <  2:2502478/100‑1 (MQ=38)
               tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATGGTGTAGTGACTATCAGACCCCTGATACGTCACg                                                                                      >  2:354074/1‑100 (MQ=255)
               tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCAc                                                                                       <  1:298679/99‑1 (MQ=38)
               tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACg                                                                                      >  1:2979990/1‑100 (MQ=255)
                gccTTACTGCTATAATTCCGAGGCCAAGAAACGTCATAAGGGAGCACAGAACCGGATGGACAAATCGCGACGTGAACAGCAGGCCCCTGAGACGTCa                                                                                        >  2:3565417/3‑97 (MQ=255)
                  cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccc                                                                                   >  1:2720088/1‑100 (MQ=255)
                  cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAACGTGTCGTGACGTAGCAGGAGCCTGATAGTCACgaca                                                                                   <  2:148259/100‑4 (MQ=37)
                   ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacac                                                                                  <  2:1120376/100‑5 (MQ=25)
                   ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacac                                                                                  <  2:114219/100‑5 (MQ=37)
                   ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGcccc                                                                                  >  1:246537/1‑100 (MQ=255)
                    tACTGACCTAATTCTCATTCCATGAAACGTCATAATGTTGAACAGACACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacacg                                                                                 <  1:1499321/100‑6 (MQ=21)
                    tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccccc                                                                                 >  2:2567986/1‑100 (MQ=255)
                      cTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGCAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTg                                                                               >  2:150323/1‑100 (MQ=255)
                      cTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTg                                                                               >  1:794600/1‑100 (MQ=255)
                       tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacacgatt                                                                              <  1:3169092/100‑9 (MQ=21)
                       tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGGCGTGACGTATCAGGGGCCTGATAGTCACgacacgatt                                                                              <  1:28926/100‑9 (MQ=21)
                        gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                             >  2:349383/1‑100 (MQ=255)
                        gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                             >  1:271029/1‑100 (MQ=255)
                          cATAATTCTCATTCCATGACACGTCATAATGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                            <  2:720289/99‑1 (MQ=255)
                            tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                         >  1:2304457/1‑100 (MQ=255)
                             aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                         >  1:2868589/1‑99 (MQ=255)
                             aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                        <  1:3502928/100‑1 (MQ=255)
                               ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCtt                                                                      <  2:1204232/100‑1 (MQ=38)
                                 ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                    <  1:789885/100‑1 (MQ=255)
                                 ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                    >  1:1534484/1‑100 (MQ=255)
                                  tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                   <  1:2630833/100‑1 (MQ=255)
                                   cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCGCCTGATACGTCACGCCCCCTGGAGAGGCTTCcga                                                                  >  2:1010930/1‑100 (MQ=255)
                                   cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga                                                                  >  2:2012709/1‑100 (MQ=255)
                                    aTTCCATGAAACGTCATAAGGTTGAACAGAAACGTATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGGTTCagag                                                                 <  2:1907885/100‑1 (MQ=255)
                                    aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                 <  2:2827653/100‑1 (MQ=255)
                                    aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                 >  1:625474/1‑100 (MQ=255)
                                    aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                 >  2:827780/1‑100 (MQ=255)
                                       ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtca                                                              <  2:848678/100‑1 (MQ=255)
                                        cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagagggcg                                                             >  1:2232281/1‑96 (MQ=255)
                                        cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcat                                                             >  1:2169487/1‑100 (MQ=255)
                                         aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                            <  1:1520324/100‑1 (MQ=255)
                                         aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                            <  2:3201210/100‑1 (MQ=255)
                                         aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcat                                                             >  2:3468288/1‑99 (MQ=255)
                                         aTGAAACGTCATAAGGTCGAACAGAAACGGATGGACAAATCGTGTCGTGACTATGAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                            >  2:2386500/1‑100 (MQ=255)
                                          tgcaaCGTCATAAGGTTGACCAGAAAAGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                            <  1:3644323/96‑1 (MQ=255)
                                              aCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCatan                                                       <  2:1306214/100‑2 (MQ=255)
                                               cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                      >  1:2843062/1‑100 (MQ=255)
                                                   aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                  >  1:1896304/1‑100 (MQ=255)
                                                   aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                  <  1:3029512/100‑1 (MQ=255)
                                                   aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCACCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                   <  2:3262659/99‑1 (MQ=255)
                                                    tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                 >  1:3545184/1‑100 (MQ=255)
                                                    tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                 <  1:3386302/100‑1 (MQ=255)
                                                       gatctaaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                              <  2:1940897/95‑1 (MQ=255)
                                                          tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGTAGAGGCTTCAGAGTCATTCATAATAACCAACCAt                                            <  2:583827/99‑1 (MQ=255)
                                                          tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATCCATAATAACCAACCATg                                           >  2:2913730/1‑100 (MQ=255)
                                                           gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgn                                          <  1:3240306/100‑2 (MQ=255)
                                                           gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                          <  1:438991/100‑1 (MQ=255)
                                                            aaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAaccaaccagatc                                         >  1:1940907/1‑96 (MQ=255)
                                                            aaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                         >  1:1151813/1‑100 (MQ=255)
                                                             aCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGAc                                        <  2:216133/100‑1 (MQ=255)
                                                                gAAACGGATGGACAAATCGTGTCGTGCCTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGaaaa                                      >  2:865302/1‑99 (MQ=255)
                                                                gAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                     >  2:1708398/1‑100 (MQ=255)
                                                                gAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                     >  1:856674/1‑100 (MQ=255)
                                                                    cGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                 >  2:2708307/1‑100 (MQ=255)
                                                                     ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAGCCAACCATGGACAAGAAAAg                                <  2:2224440/100‑1 (MQ=255)
                                                                     ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                <  2:1462876/100‑1 (MQ=255)
                                                                     cttccgatctaTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                <  2:162496/90‑1 (MQ=255)
                                                                         ggACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                            <  2:1298576/100‑1 (MQ=255)
                                                                              aaTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                       >  1:177584/1‑100 (MQ=255)
                                                                                tcgtgtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAt                     <  1:836050/100‑1 (MQ=255)
                                                                                 gggggcgtgACGTATCAGGGGCCTGATAGTCACGCCGCCTGGAGAGGCTTCAGAGTCATACATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                    >  2:52792/6‑100 (MQ=25)
                                                                                 cgtgtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                    >  1:437328/1‑100 (MQ=255)
                                                                                  gtgttgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                    <  1:1807694/99‑1 (MQ=255)
                                                                                   tgtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtctc                  <  1:2359973/100‑1 (MQ=255)
                                                                                    gtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                 <  1:3536882/100‑1 (MQ=255)
                                                                                    gtcgtgACTATCAGGCCCCCGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAAACAACCATGGACAAGAAAAGTAAGATTCAATctct                 <  1:1145241/100‑1 (MQ=255)
                                                                                    ggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                 >  1:1029726/3‑100 (MQ=37)
                                                                                    ggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                 >  1:2574354/3‑100 (MQ=37)
                                                                                     tcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                <  2:444730/100‑1 (MQ=255)
                                                                                     gcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGCGTCATTCATAATAACCCACCATGGACACGAAAAGTAAGACTCAATCTCTg                >  2:1533340/2‑100 (MQ=25)
                                                                                      cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTgt               >  1:3070965/1‑99 (MQ=37)
                                                                                      cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc               >  2:3354571/1‑100 (MQ=38)
                                                                                      cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc               >  1:1741671/1‑100 (MQ=38)
                                                                                      cgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc               >  1:1642467/1‑100 (MQ=255)
                                                                                       gtgacgtATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa              >  1:339017/7‑100 (MQ=38)
                                                                                       gtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa              <  1:1960849/100‑1 (MQ=255)
                                                                                        tgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag             >  1:3255970/1‑100 (MQ=255)
                                                                                        tgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag             <  1:419310/100‑1 (MQ=255)
                                                                                             tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc         <  1:2346382/100‑1 (MQ=255)
                                                                                             tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc         <  1:3579709/100‑1 (MQ=255)
                                                                                              aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa        <  2:3310599/100‑1 (MQ=255)
                                                                                              aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa        <  1:3621750/100‑1 (MQ=255)
                                                                                                    ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGACTCATTCATAATAAACAACCATGGACAAGAAAAGTAAGATTCAATATCTACACAGACAAGAAc   >  2:3184315/1‑99 (MQ=255)
                                                                                                    |                                                                                                   
CAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACC  >  NC_000913/2449516‑2449714

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: