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.

AGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCT  >  NC_000913/2449517‑2449715
                                                                                                    |                                                                                                    
 gggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGccc                                                                                                     >  1:1010479/1‑100 (MQ=255)
   gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCt                                                                                                   >  2:3621524/1‑100 (MQ=255)
    tGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTg                                                                                                  <  1:2518145/100‑1 (MQ=255)
      aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGAt                                                                                                <  2:3492226/100‑1 (MQ=255)
       gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACCGAAACGGATGGACCAATCGTGTCGTGACTATCAGGCCCCTGAtc                                                                                               >  2:2514648/1‑99 (MQ=255)
       gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATa                                                                                               >  2:789880/1‑100 (MQ=255)
            tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGGGCCTGATAGTCAc                                                                                         <  1:436043/100‑1 (MQ=255)
            tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCa                                                                                          <  1:1732054/100‑1 (MQ=255)
            tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGGGACGTATCAGGGGCCTGATagtca                                                                                           <  2:1890675/100‑5 (MQ=255)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATGGTGTAGTGACTATCAGACCCCTGATACGTCACg                                                                                        >  2:354074/1‑100 (MQ=255)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCAc                                                                                         <  1:298679/99‑1 (MQ=38)
                 cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgaca                                                                                     <  1:2024176/100‑4 (MQ=37)
                 cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAACGTGTCGTGACGTAGCAGGAGCCTGATAGTCACgaca                                                                                     <  2:148259/100‑4 (MQ=37)
                  ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacac                                                                                    <  1:666110/100‑5 (MQ=37)
                  ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacac                                                                                    <  2:1120376/100‑5 (MQ=25)
                   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)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                              >  1:530225/1‑100 (MQ=255)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                              >  1:271029/1‑100 (MQ=255)
                        aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGTCCCCTGgag                                                                             >  2:2766639/1‑100 (MQ=255)
                         cATAATTCTCATTCCATGACACGTCATAATGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                             <  2:720289/99‑1 (MQ=255)
                         cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                             >  1:288171/1‑99 (MQ=255)
                          aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                           <  1:411752/100‑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)
                               tctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTc                                                                      <  2:1664180/100‑1 (MQ=255)
                                ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                     >  1:1534484/1‑100 (MQ=255)
                                ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                     <  1:789885/100‑1 (MQ=255)
                                 tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                    <  2:2954087/100‑1 (MQ=255)
                                  cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCGCCTGATACGTCACGCCCCCTGGAGAGGCTTCcga                                                                   >  2:1010930/1‑100 (MQ=255)
                                  cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga                                                                   <  2:1198061/100‑1 (MQ=255)
                                  cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga                                                                   >  2:2012709/1‑100 (MQ=255)
                                   aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                  <  2:2827653/100‑1 (MQ=255)
                                   aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                  >  2:827780/1‑100 (MQ=255)
                                      ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtca                                                               <  2:848678/100‑1 (MQ=255)
                                      ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtca                                                               <  2:3232649/100‑1 (MQ=255)
                                       cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagagggcg                                                              >  1:2232281/1‑96 (MQ=255)
                                       cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcat                                                              >  1:2699657/1‑100 (MQ=255)
                                       cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcat                                                              >  1:2169487/1‑100 (MQ=255)
                                        aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                             <  1:1520324/100‑1 (MQ=255)
                                        aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                             <  1:2278413/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)
                                           aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                          <  1:3329049/100‑1 (MQ=255)
                                           aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                          <  2:1931304/100‑1 (MQ=255)
                                              cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                       <  1:3508339/100‑1 (MQ=255)
                                                 cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                    <  1:3046308/100‑1 (MQ=255)
                                                  aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                   >  1:1896304/1‑100 (MQ=255)
                                                  aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCACCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                    <  2:3262659/99‑1 (MQ=255)
                                                   tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                  >  1:3545184/1‑100 (MQ=255)
                                                   tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                  <  1:2599526/100‑1 (MQ=255)
                                                     aGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                <  2:3113651/100‑1 (MQ=255)
                                                      gatctaaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                               <  2:1940897/95‑1 (MQ=255)
                                                       gTTGAACAGAAACGGATGGACAAATCGTGTCGTGAATATCAGGCCCCTGATACGTCACGCCCCCTGGAGGGGCTTCAGAGTCATTCATAATAaccaacca                                              >  1:3208648/1‑100 (MQ=255)
                                                         tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGTAGAGGCTTCAGAGTCATTCATAATAACCAACCAt                                             <  2:583827/99‑1 (MQ=255)
                                                          gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                           >  1:2456017/1‑100 (MQ=255)
                                                          gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                           <  1:438991/100‑1 (MQ=255)
                                                           aaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAaccaaccagatc                                          >  1:1940907/1‑96 (MQ=255)
                                                           aaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                          >  1:1151813/1‑100 (MQ=255)
                                                               gAAACGGATGGACAAATCGTGTCGTGCCTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGaaaa                                       >  2:865302/1‑99 (MQ=255)
                                                               gAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                      >  2:1708398/1‑100 (MQ=255)
                                                               gAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                      <  2:1337530/100‑1 (MQ=255)
                                                                   cGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                  >  2:2708307/1‑100 (MQ=255)
                                                                    ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAGCCAACCATGGACAAGAAAAg                                 <  2:2224440/100‑1 (MQ=255)
                                                                    ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                 <  2:1462876/100‑1 (MQ=255)
                                                                        ggACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                             <  2:1298576/100‑1 (MQ=255)
                                                                         gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                             >  2:3388504/1‑100 (MQ=255)
                                                                         gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                            <  2:1010472/100‑1 (MQ=255)
                                                                          aCAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAt                           <  1:3575/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)
                                                                                 gtgttgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                     <  1:1807694/99‑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)
                                                                                     cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGCGGCTTCAGAGCCATTGATACTCCCCAACCAGGAACACGAAAAGTAAGATTCAAtctc                   >  2:3510698/1‑97 (MQ=16)
                                                                                     cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTgt                >  1:3070965/1‑99 (MQ=37)
                                                                                     cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                >  1:1741671/1‑100 (MQ=38)
                                                                                     cgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                <  1:298870/100‑1 (MQ=255)
                                                                                     cgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                >  1:1642467/1‑100 (MQ=255)
                                                                                      gtgacgtATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa               >  1:339017/7‑100 (MQ=38)
                                                                                       tgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag              >  1:3255970/1‑100 (MQ=255)
                                                                                       tgACTATCAGGCCCATGATACGTCACGCCCCCTGGAGAGGCATCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag              <  2:3634879/100‑1 (MQ=255)
                                                                                            tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc          <  1:2346382/100‑1 (MQ=255)
                                                                                             aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa         <  1:3621551/100‑1 (MQ=255)
                                                                                             aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa         <  2:3310599/100‑1 (MQ=255)
                                                                                               cAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAAt       >  1:316198/1‑100 (MQ=255)
                                                                                                aGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATa      <  2:2201179/100‑1 (MQ=255)
                                                                                                   ccccTGATACGTCACGTCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc   <  1:2766657/100‑1 (MQ=255)
                                                                                                   ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc   >  1:294567/1‑100 (MQ=255)
                                                                                                   ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGACTCATTCATAATAAACAACCATGGACAAGAAAAGTAAGATTCAATATCTACACAGACAAGAAc    >  2:3184315/1‑99 (MQ=255)
                                                                                                    cccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCt  <  1:2383794/100‑1 (MQ=255)
                                                                                                    |                                                                                                    
AGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCT  >  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: