Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,616 C→G 17.2% G145A (GGG→GCG) ‡ 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,6160CG17.2% 401.0 / 69.2 156G145A (GGG→GCG) yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (56/72);  new base G (13/14);  total (69/87)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

Only 100 of 165 total aligned reads displayed.

AGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCT  >  NC_000913/2449517‑2449715
                                                                                                    |                                                                                                    
 gggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGccc                                                                                                     >  1:1010479/1‑100 (MQ=255)
   gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCt                                                                                                   >  2:3621524/1‑100 (MQ=255)
    tGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTg                                                                                                  <  1:2518145/100‑1 (MQ=255)
      aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGAt                                                                                                <  2:3492226/100‑1 (MQ=255)
       gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACCGAAACGGATGGACCAATCGTGTCGTGAC‑TATCAGGCCCCTGAtc                                                                                               >  2:2514648/1‑99 (MQ=255)
       gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATa                                                                                               >  2:789880/1‑100 (MQ=255)
            tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGGGCCTGATA‑GTCAc                                                                                         <  1:436043/100‑1 (MQ=255)
            tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCa                                                                                          <  1:1732054/100‑1 (MQ=255)
            tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGGGACGTATCAGGGGCCTGATagtca                                                                                           <  2:1890675/100‑5 (MQ=255)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATGGTGTAGTGAC‑TATCAGACCCCTGATACGTCACg                                                                                        >  2:354074/1‑100 (MQ=255)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCAc                                                                                         <  1:298679/99‑1 (MQ=38)
                 cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgaca                                                                                     <  1:2024176/100‑4 (MQ=37)
                 cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAACGTGTCGTGACGTAGCAGGAGCCTGATA‑GTCACgaca                                                                                     <  2:148259/100‑4 (MQ=37)
                  ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                                    <  1:666110/100‑5 (MQ=37)
                  ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                                    <  2:1120376/100‑5 (MQ=25)
                   tACTGACCTAATTCTCATTCCATGAAACGTCATAATGTTGAACAGACACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacacg                                                                                   <  1:1499321/100‑6 (MQ=21)
                   tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑ccccc                                                                                  >  2:2567986/1‑100 (MQ=255)
                     cTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGCAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTg                                                                                >  2:150323/1‑100 (MQ=255)
                     cTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTg                                                                                >  1:794600/1‑100 (MQ=255)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGga                                                                              >  1:530225/1‑100 (MQ=255)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGga                                                                              >  1:271029/1‑100 (MQ=255)
                        aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑TCCCCTGgag                                                                             >  2:2766639/1‑100 (MQ=255)
                         cATAATTCTCATTCCATGACACGTCATAATGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGgag                                                                             <  2:720289/99‑1 (MQ=255)
                         cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGgag                                                                             >  1:288171/1‑99 (MQ=255)
                          aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGagag                                                                           <  1:411752/100‑1 (MQ=255)
                           tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAgg                                                                          >  1:2304457/1‑100 (MQ=255)
                            aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAgg                                                                          >  1:2868589/1‑99 (MQ=255)
                            aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGc                                                                         <  1:3502928/100‑1 (MQ=255)
                              ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCtt                                                                       <  2:1204232/100‑1 (MQ=38)
                               tctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTc                                                                      <  2:1664180/100‑1 (MQ=255)
                                ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCa                                                                     >  1:1534484/1‑100 (MQ=255)
                                ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCa                                                                     <  1:789885/100‑1 (MQ=255)
                                 tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCag                                                                    <  2:2954087/100‑1 (MQ=255)
                                  cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCGCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCcga                                                                   >  2:1010930/1‑100 (MQ=255)
                                  cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCaga                                                                   <  2:1198061/100‑1 (MQ=255)
                                  cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCaga                                                                   >  2:2012709/1‑100 (MQ=255)
                                   aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCagag                                                                  <  2:2827653/100‑1 (MQ=255)
                                   aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCagag                                                                  >  2:827780/1‑100 (MQ=255)
                                      ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtca                                                               <  2:848678/100‑1 (MQ=255)
                                      ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtca                                                               <  2:3232649/100‑1 (MQ=255)
                                       cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCagagggcg                                                              >  1:2232281/1‑96 (MQ=255)
                                       cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcat                                                              >  1:2699657/1‑100 (MQ=255)
                                       cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcat                                                              >  1:2169487/1‑100 (MQ=255)
                                        aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcatt                                                             <  1:1520324/100‑1 (MQ=255)
                                        aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcatt                                                             <  1:2278413/100‑1 (MQ=255)
                                        aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcat                                                              >  2:3468288/1‑99 (MQ=255)
                                        aTGAAACGTCATAAGGTCGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATGAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcatt                                                             >  2:2386500/1‑100 (MQ=255)
                                         tgcaaCGTCATAAGGTTGACCAGAAAAGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcatt                                                             <  1:3644323/96‑1 (MQ=255)
                                           aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcattcat                                                          <  1:3329049/100‑1 (MQ=255)
                                           aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcattcat                                                          <  2:1931304/100‑1 (MQ=255)
                                              cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                       <  1:3508339/100‑1 (MQ=255)
                                                 cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                    <  1:3046308/100‑1 (MQ=255)
                                                  aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                   >  1:1896304/1‑100 (MQ=255)
                                                  aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCACCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                    <  2:3262659/99‑1 (MQ=255)
                                                   tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                  >  1:3545184/1‑100 (MQ=255)
                                                   tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                  <  1:2599526/100‑1 (MQ=255)
                                                     aGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                <  2:3113651/100‑1 (MQ=255)
                                                      gatctaaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                               <  2:1940897/95‑1 (MQ=255)
                                                       gTTGAACAGAAACGGATGGACAAATCGTGTCGTGAA‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGGGGCTTCAGAGTCATTCATAATAaccaacca                                              >  1:3208648/1‑100 (MQ=255)
                                                         tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGTAGAGGCTTCAGAGTCATTCATAATAACCAACCAt                                             <  2:583827/99‑1 (MQ=255)
                                                          gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                           >  1:2456017/1‑100 (MQ=255)
                                                          gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                           <  1:438991/100‑1 (MQ=255)
                                                           aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAaccaaccagatc                                          >  1:1940907/1‑96 (MQ=255)
                                                           aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                          >  1:1151813/1‑100 (MQ=255)
                                                               gAAACGGATGGACAAATCGTGTCGTGCC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGaaaa                                       >  2:865302/1‑99 (MQ=255)
                                                               gAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                      >  2:1708398/1‑100 (MQ=255)
                                                               gAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                      <  2:1337530/100‑1 (MQ=255)
                                                                   cGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                  >  2:2708307/1‑100 (MQ=255)
                                                                    ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAGCCAACCATGGACAAGAAAAg                                 <  2:2224440/100‑1 (MQ=255)
                                                                    ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                 <  2:1462876/100‑1 (MQ=255)
                                                                        ggACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                             <  2:1298576/100‑1 (MQ=255)
                                                                         gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                             >  2:3388504/1‑100 (MQ=255)
                                                                         gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                            <  2:1010472/100‑1 (MQ=255)
                                                                          aCAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAt                           <  1:3575/100‑1 (MQ=255)
                                                                             aaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                        >  1:177584/1‑100 (MQ=255)
                                                                               tcgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAt                      <  1:836050/100‑1 (MQ=255)
                                                                                gggggcgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCGCCTGGAGAGGCTTCAGAGTCATACATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                     >  2:52792/6‑100 (MQ=25)
                                                                                 gtgttgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                     <  1:1807694/99‑1 (MQ=255)
                                                                                   gtcgtgAC‑TATCAGGCCCCCGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAAACAACCATGGACAAGAAAAGTAAGATTCAATctct                  <  1:1145241/100‑1 (MQ=255)
                                                                                   ggcgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                  >  1:1029726/3‑100 (MQ=37)
                                                                                   ggcgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                  >  1:2574354/3‑100 (MQ=37)
                                                                                    tcgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                 <  2:444730/100‑1 (MQ=255)
                                                                                    gcgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGCGTCATTCATAATAACCCACCATGGACACGAAAAGTAAGACTCAATCTCTg                 >  2:1533340/2‑100 (MQ=25)
                                                                                     cgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGCGGCTTCAGAGCCATTGATACTCCCCAACCAGGAACACGAAAAGTAAGATTCAAtctc                   >  2:3510698/1‑97 (MQ=16)
                                                                                     cgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTgt                >  1:3070965/1‑99 (MQ=37)
                                                                                     cgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                >  1:1741671/1‑100 (MQ=38)
                                                                                     cgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                <  1:298870/100‑1 (MQ=255)
                                                                                     cgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                >  1:1642467/1‑100 (MQ=255)
                                                                                      gtgacgtATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa               >  1:339017/7‑100 (MQ=38)
                                                                                       tgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag              >  1:3255970/1‑100 (MQ=255)
                                                                                       tgAC‑TATCAGGCCCATGATACGTCACG‑CCCCCTGGAGAGGCATCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag              <  2:3634879/100‑1 (MQ=255)
                                                                                            tATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc          <  1:2346382/100‑1 (MQ=255)
                                                                                             aTCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa         <  1:3621551/100‑1 (MQ=255)
                                                                                             aTCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa         <  2:3310599/100‑1 (MQ=255)
                                                                                               cAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAAt       >  1:316198/1‑100 (MQ=255)
                                                                                                aGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATa      <  2:2201179/100‑1 (MQ=255)
                                                                                                   ccccTGATACGTCACG‑TCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc   <  1:2766657/100‑1 (MQ=255)
                                                                                                   ccccTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc   >  1:294567/1‑100 (MQ=255)
                                                                                                   ccccTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGACTCATTCATAATAAACAACCATGGACAAGAAAAGTAAGATTCAATATCTACACAGACAAGAAc    >  2:3184315/1‑99 (MQ=255)
                                                                                                    cccTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCt  <  1:2383794/100‑1 (MQ=255)
                                                                                                    |                                                                                                    
AGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCT  >  NC_000913/2449517‑2449715

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: ‑