Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,615 C→G 18.4% 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,6150CG18.4% 406.7 / 75.8 157G145G (GGG→GGC) yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (62/65);  new base G (17/12);  total (79/78)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

AGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACC  >  NC_000913/2449517‑2449714
                                                                                                   |                                                                                                   
aGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGcc                                                                                                    >  2:2678154/1‑100 (MQ=255)
 gggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTAATAAGGTTGAACAGAAACGGATGGACAAATCTTGTCGTGAC‑TATCAGGccc                                                                                                   >  1:959543/1‑100 (MQ=255)
  ggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGcccc                                                                                                  <  1:901430/100‑1 (MQ=255)
    tGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTg                                                                                                <  1:1224351/100‑1 (MQ=255)
     gAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGa                                                                                               <  1:3348116/100‑1 (MQ=255)
     gAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGa                                                                                               <  2:232881/100‑1 (MQ=255)
      aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGAt                                                                                              >  1:2809182/1‑100 (MQ=255)
      aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGAt                                                                                              <  1:1771475/100‑1 (MQ=255)
       gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATa                                                                                             >  1:2203173/1‑100 (MQ=255)
       gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATa                                                                                             <  2:3366295/100‑1 (MQ=255)
       gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATa                                                                                             <  1:774529/100‑1 (MQ=255)
         aGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATa                                                                                             <  1:1497804/98‑1 (MQ=255)
           ttttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTc                                                                                         >  2:1303463/1‑100 (MQ=255)
             ttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCAc                                                                                       <  1:997625/100‑1 (MQ=38)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACg                                                                                      <  1:1735186/100‑1 (MQ=38)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACg                                                                                      <  1:2547498/100‑1 (MQ=38)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACg                                                                                      <  1:3517187/100‑1 (MQ=255)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACg                                                                                      >  2:1239952/1‑100 (MQ=255)
               gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGc                                                                                     <  1:2678171/100‑1 (MQ=255)
               gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGc                                                                                     >  2:2943417/1‑100 (MQ=255)
                 cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgaca                                                                                   <  1:2853232/100‑4 (MQ=37)
                  ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                                  <  1:6712/100‑5 (MQ=25)
                     cTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTg                                                                               >  1:780538/1‑100 (MQ=255)
                      tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                              <  2:2639178/100‑1 (MQ=255)
                      tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                              <  1:1204780/100‑1 (MQ=255)
                      tGACATAATGCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                              <  2:2759161/100‑1 (MQ=255)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                             >  1:1239652/1‑100 (MQ=255)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                             >  1:2313498/1‑100 (MQ=255)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                             >  1:2519122/1‑100 (MQ=255)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                             >  1:1804157/1‑100 (MQ=255)
                        aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                            <  1:1969159/100‑1 (MQ=255)
                         cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                           >  2:591113/1‑100 (MQ=255)
                         cATAATTCTCATGCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                           >  1:2602716/1‑100 (MQ=255)
                          aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                          <  2:1667373/100‑1 (MQ=255)
                          aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                          <  1:1654804/100‑1 (MQ=255)
                           tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                         >  1:2881204/1‑100 (MQ=255)
                            aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                        <  1:281132/100‑1 (MQ=255)
                            aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                        >  2:344106/1‑100 (MQ=255)
                            aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                        <  1:3461728/100‑1 (MQ=255)
                             aTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                       <  2:3947666/100‑1 (MQ=255)
                              ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCtt                                                                      >  1:767145/1‑100 (MQ=255)
                                ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCtt                                                                      >  2:113553/1‑98 (MQ=255)
                                  cATTCCATGAAACGTCATAAGGTGGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                    <  1:3090646/98‑1 (MQ=255)
                                   aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGCGAGGCTTCagag                                                                 >  1:1139014/1‑100 (MQ=255)
                                      ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtca                                                              >  1:2513922/1‑100 (MQ=255)
                                         tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATAAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                           >  2:378735/1‑100 (MQ=255)
                                              cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataa                                                       <  2:416019/99‑1 (MQ=255)
                                               gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataata                                                     <  1:3544301/100‑1 (MQ=255)
                                                tCATAAGGTTGAACAGAAATGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCAtaataa                                                    >  2:2399574/1‑100 (MQ=255)
                                                tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTTACGCCCCCTGGAGAGGCTTCACAGTCATTCataata                                                     >  1:3880160/1‑99 (MQ=255)
                                                 cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                   <  1:790392/100‑1 (MQ=255)
                                                 cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                   <  1:1448214/100‑1 (MQ=255)
                                                 cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATaac                                                   <  2:2443236/100‑1 (MQ=255)
                                                  aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                  <  2:1045733/100‑1 (MQ=255)
                                                   tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                 >  2:2390220/1‑100 (MQ=255)
                                                     agcttGAACAGAACCGGACGTACAAATCGTGTCCTGAC‑TATCAGGCCCCTGATACCTCACGCCCCCTGGAGACGCTTCAGAGTCATTCATAATaac                                                   <  2:2547191/93‑1 (MQ=255)
                                                       gTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAaccaacca                                             >  1:3903357/1‑100 (MQ=255)
                                                        ttGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCGCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAGCCAt                                            <  2:1069352/100‑1 (MQ=255)
                                                         tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                           <  2:3825569/100‑1 (MQ=255)
                                                         tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                           >  2:3061869/1‑100 (MQ=255)
                                                         tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCACCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                           <  2:2161353/100‑1 (MQ=255)
                                                          gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                          <  1:508995/100‑1 (MQ=255)
                                                          gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                          >  1:2108499/1‑100 (MQ=255)
                                                           aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                         <  1:186034/100‑1 (MQ=255)
                                                           aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTGCAGAGTCATTCATAATAACCCACCATGGa                                         >  2:1963663/1‑100 (MQ=255)
                                                           aaCAGAAACGGATGGACAAATCGTGACGTGAC‑TAGCAGGCCCCTGATACGTGACGCCCCCTGGAGAGGCTTCAGCGTCATTAATAATAACTAGCCCTGGa                                         >  1:662334/1‑100 (MQ=255)
                                                             cAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACa                                       <  1:2154373/100‑1 (MQ=255)
                                                               gAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                     >  1:2938214/1‑100 (MQ=255)
                                                                  aCGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaa                                  <  1:1295909/100‑1 (MQ=255)
                                                                    ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                <  1:922304/100‑1 (MQ=255)
                                                                      aTGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTa                              >  1:1602642/1‑100 (MQ=255)
                                                                        ggACAAATCGTGTCGTTAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                            <  2:1534791/100‑1 (MQ=255)
                                                                           cAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                         <  1:1337307/100‑1 (MQ=255)
                                                                            aaaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTc                        <  2:1891849/100‑1 (MQ=255)
                                                                             aaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                       >  1:102483/1‑100 (MQ=255)
                                                                             aaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                       <  2:2582891/100‑1 (MQ=255)
                                                                              aTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCACCCTGGAGAGGCTTCAGAGTAATTCCTAATAACCACCAAGGGACAACAAAAGTAAGATCCaa                      >  2:659508/1‑100 (MQ=255)
                                                                                cgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                    <  1:287470/100‑1 (MQ=255)
                                                                                 ggggcgtgACGTATCAGGGGCCTGATA‑GTCACACCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct                   >  1:3743833/5‑100 (MQ=25)
                                                                                   gtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                 >  2:3517687/1‑100 (MQ=255)
                                                                                   ggcgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                 >  1:3115337/3‑100 (MQ=37)
                                                                                   ggcgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                 >  1:3613679/3‑100 (MQ=37)
                                                                                    tcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                <  2:2313486/100‑1 (MQ=255)
                                                                                     cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc               >  1:3622871/1‑100 (MQ=38)
                                                                                     cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc               >  1:1498306/1‑100 (MQ=38)
                                                                                     cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc               >  1:3479491/1‑100 (MQ=38)
                                                                                      gtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa              <  2:3553864/100‑1 (MQ=255)
                                                                                      gtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa              <  1:1908098/100‑1 (MQ=255)
                                                                                      gtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa              <  1:2390232/100‑1 (MQ=255)
                                                                                       tgacgtATCAGGGGCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa              >  1:929948/6‑100 (MQ=255)
                                                                                       tgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag             <  2:3818553/100‑1 (MQ=255)
                                                                                       tgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag             <  2:911788/100‑1 (MQ=255)
                                                                                         acgtATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag           >  2:3199402/4‑100 (MQ=255)
                                                                                         aC‑TCTCCGGCCCCTGATACGTCACGACCCCTGGAGAGGCCTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag             <  1:3144647/98‑1 (MQ=255)
                                                                                         aC‑TATCAGGACCCTGATACGTCACGCCACCTGGAGAGGCCTCATAGCCATTCATAATAATCAACCATGGACAAGAAAAGTAAGATGCAATCTCTGCagag           <  2:416174/100‑1 (MQ=255)
                                                                                             aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa        >  2:1231097/1‑100 (MQ=255)
                                                                                              tCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACaa       >  1:3470469/1‑100 (MQ=255)
                                                                                                  gCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAc   >  1:3030596/1‑100 (MQ=255)
                                                                                                   ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc  >  1:2349936/1‑100 (MQ=255)
                                                                                                   ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc  >  2:543307/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACC  >  NC_000913/2449517‑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: ‑