Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 380,021 (G)10→8 100% intergenic (‑140/+47) frmR ← / ← yaiO regulator protein that represses frmRAB operon/outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,0130G.95.6% 587.6 / 16.5 136intergenic (‑132/+56)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base . (59/71);  total (63/73)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01
*NC_000913380,0140G.75.0% 306.9 / 149.2 136intergenic (‑133/+55)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (18/14);  new base . (44/58);  total (63/73)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 168 total aligned reads displayed.

TAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAA  >  NC_000913/379914‑380113
                                                                                                   ||                                                                                                   
taggggggtatgcattgacatatagaatacccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAg                                                                                                      >  1:1265451‑M2/46‑100 (MQ=255)
  ggggggtatgcattgacatatagaatacccccctatagtatatAACACTGGAGAATAAAATTTCTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                                    >  2:1702924‑M2/44‑100 (MQ=255)
  ggggggtatgcattgacatatagaatacccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                                    >  1:3533546‑M2/44‑100 (MQ=255)
  ggggggtatgcattgacatatagaatacccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                                    >  1:1465965‑M2/44‑100 (MQ=255)
  ggggggtatgcattgacatatagaatacccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                                    >  2:3511882‑M2/44‑100 (MQ=255)
  GGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                                    <  2:3481097‑M1/100‑1 (MQ=255)
  GGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                                    <  2:3160127‑M1/100‑1 (MQ=255)
  GGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                                    <  1:1323738‑M1/100‑1 (MQ=255)
  GGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTAACCGGTGAGTGTGGTCGGAAAACAAAGAGGAAAg                                                                                                      <  1:3555550‑M1/98‑1 (MQ=255)
   gggggtatgcattgacatatagaatacccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                                   >  1:1234525‑M2/43‑100 (MQ=255)
   gggggtatgcattgacatatagaatacccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                                   >  2:750113‑M2/43‑100 (MQ=255)
    ggggtatgcattgacatatagaatacccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAAGAAAcggg                                                                                                   >  2:1397502‑M2/42‑99 (MQ=255)
    ggggtatgcattgacatatagaatacccccctatagtatatAACACTGGAGAATAAAATTTATCAGGTGAATGTAGTCGGAAAACAAAGAGGAAAggggg                                                                                                  >  2:2538108‑M2/42‑100 (MQ=255)
    GGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                                                                                  <  2:652782‑M1/100‑1 (MQ=255)
     GGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggg                                                                                                 <  2:2425403‑M1/100‑1 (MQ=255)
      GGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGCCGGAAAACAAAGAGGAAAggggggg                                                                                                <  2:2545362‑M1/100‑1 (MQ=255)
        TATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAaggggggggc                                                                                              >  1:1104434‑M1/1‑99 (MQ=255)
          TTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAaagggggggggct                                                                                            >  1:2506900‑M1/1‑98 (MQ=255)
          TTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTGATCCGGTGAATGTGGTCGGAAAACAAAGAGGAaaggggggggct                                                                                             >  2:1541310‑M1/1‑97 (MQ=255)
           tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggctaa                                                                                           >  2:3544631/1‑96 (MQ=255)
                atGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCgg                                                                                     >  2:2372436/1‑100 (MQ=255)
                atGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGAGGCTAATCgg                                                                                     >  2:1483960/1‑99 (MQ=255)
                 tGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAACTTTATCCGGTGAATGTGGTCGGAAAACAACGAGGAA‑‑GGGGGGGGGCTAATCGGCa                                                                                   >  2:1216006/1‑100 (MQ=255)
                 tGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGc                                                                                    <  2:26974/100‑1 (MQ=255)
                  ntCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAg                                                                                  >  1:2689213/2‑100 (MQ=255)
                  gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCa                                                                                   >  1:2440150/1‑100 (MQ=255)
                  gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAg                                                                                  >  2:3748122/1‑100 (MQ=255)
                  gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAg                                                                                  >  2:1351954/1‑100 (MQ=255)
                  gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGAGGGGGCTAATCGGCAg                                                                                  >  2:2233001/1‑100 (MQ=255)
                   tCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAg                                                                                  <  1:3010057/100‑1 (MQ=255)
                     aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAggg                                                                                <  1:3385747/99‑1 (MQ=255)
                     aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGa                                                                               >  1:1404609/1‑100 (MQ=255)
                       tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAg                                                                             <  2:21139/100‑1 (MQ=255)
                         cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGgc                                                                           <  1:2662500/100‑1 (MQ=255)
                          atacccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGgcc                                                                          <  1:601348‑M2/81‑1 (MQ=255)
                          atacccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGgc                                                                           <  2:1308225‑M2/80‑1 (MQ=255)
                            acccccctatagtatatAACACTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGgcc                                                                          <  1:2962906‑M2/81‑1 (MQ=255)
                            aacaccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑AGGGGGGGCTAATCGGCAGGGAAGgccgc                                                                        <  2:649838‑M2/83‑1 (MQ=255)
                             cccccctatagtatatAACACTGGAGCATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccgc                                                                        >  2:3526062‑M2/17‑100 (MQ=255)
                             cccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccgc                                                                        >  1:197576‑M2/17‑100 (MQ=255)
                             cccccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGgccgcc                                                                       <  1:504228‑M2/84‑1 (MQ=255)
                               cccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                                     <  1:3090559‑M2/86‑1 (MQ=255)
                               cccctatagtatatAACACTGGAGAATAAAAGTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                                     <  2:1896483‑M2/86‑1 (MQ=255)
                                ccctatagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                                     <  1:2485229‑M2/87‑1 (MQ=255)
                                     tagtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGTGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAg                                                               >  1:1059448‑M2/9‑100 (MQ=255)
                                       gtatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGCCCGCCCCGGATAGc                                                              >  1:1761907‑M2/7‑100 (MQ=255)
                                        tatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCgg                                                            <  1:822045‑M2/95‑1 (MQ=255)
                                        tatatAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGCAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCgg                                                            <  2:1169174‑M2/95‑1 (MQ=255)
                                           atAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGc                                                          <  1:2299618‑M2/98‑1 (MQ=255)
                                           atAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCg                                                         <  1:2923488‑M2/98‑1 (MQ=255)
                                           atAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCg                                                         <  2:1869264‑M2/98‑1 (MQ=255)
                                            taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg                                                        >  2:2877831/1‑100 (MQ=255)
                                            taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGggcggc                                                       >  2:242906/1‑100 (MQ=255)
                                             aacgcgggAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggc                                                          <  2:1627832/92‑1 (MQ=255)
                                             aaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCa                                                      <  1:3160151/100‑1 (MQ=255)
                                              acacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg                                                        >  2:2186711/1‑100 (MQ=255)
                                                  tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCTCGGATAGCGGGCGGCAGAAg                                                  <  2:2333120/100‑1 (MQ=255)
                                                     agaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGa                                                <  2:2170401/100‑1 (MQ=255)
                                                     agaATAAAATTTAGCGGGGGAAGTGGGGCGGAAAAAAAAAGGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAg                                                  <  2:2935542/99‑1 (MQ=255)
                                                      gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGaa                                               >  1:1192051/1‑100 (MQ=255)
                                                       aaTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATc                                             >  1:3716294/1‑100 (MQ=255)
                                                       aaTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATc                                             <  1:94031/100‑1 (MQ=255)
                                                        aTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCa                                            >  2:500664/1‑100 (MQ=255)
                                                         tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAg                                           <  1:3664465/100‑1 (MQ=255)
                                                         tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAg                                           <  1:2216462/100‑1 (MQ=255)
                                                         tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAg                                           >  1:3181797/1‑100 (MQ=255)
                                                         tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATAAg                                           >  1:1690362/1‑100 (MQ=255)
                                                          aaaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                          <  1:2630753/100‑1 (MQ=255)
                                                           aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa                                         <  1:2820229/100‑1 (MQ=255)
                                                           aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa                                         <  2:422147/100‑1 (MQ=255)
                                                            aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGTCTAATCGGCAGGGAAGGCCGTCCCGGATAGCGGGCGGCAGAAGGAATCAGact                                        >  2:1967966/1‑98 (MQ=255)
                                                            aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCACCGGATCGCGGGCGGCCGAAGGAATCAGAAt                                        >  2:3388788/1‑100 (MQ=255)
                                                             aTTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                       <  1:2661311/100‑1 (MQ=255)
                                                                   ccGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAgg                                 <  1:2233015/100‑1 (MQ=255)
                                                                   ccGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAgg                                 <  2:1147118/100‑1 (MQ=255)
                                                                    cGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGt                                >  1:3010080/1‑100 (MQ=255)
                                                                    cGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGt                                >  1:1643350/1‑100 (MQ=255)
                                                                    cGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGt                                >  1:3776955/1‑100 (MQ=255)
                                                                       tGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                               >  2:2328722/1‑100 (MQ=255)
                                                                        tggattgttCGGAAATCAAGGAGGAATAGGGGGGGGTTAATCGTCTGAGAAGGCGGCCGCTGATAGCGCGCGGCAGAAGGAATCAGAATTTCCAgg                                 <  2:744083/88‑1 (MQ=255)
                                                                        gAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                            <  2:1727813/100‑1 (MQ=255)
                                                                             tgGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc                       >  1:2091026/1‑100 (MQ=255)
                                                                             tgGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc                       <  1:871864/100‑1 (MQ=255)
                                                                              ggTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                     >  2:1034808/1‑100 (MQ=255)
                                                                              ggTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                     >  2:1098736/1‑100 (MQ=255)
                                                                              ggTAGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCTGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                      <  1:3628106/100‑1 (MQ=255)
                                                                               gTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                     >  1:2923712/1‑100 (MQ=255)
                                                                                tCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                    >  1:723685/1‑100 (MQ=255)
                                                                                 cGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                   <  1:790173/100‑1 (MQ=255)
                                                                                 cGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                   >  1:3534448/1‑100 (MQ=255)
                                                                                 cGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                   <  1:3636748/100‑1 (MQ=255)
                                                                                  ggAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                   >  2:2107262/1‑100 (MQ=255)
                                                                                   gAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                 >  1:2495021/1‑100 (MQ=255)
                                                                                     aaaCAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt               <  2:1761894/100‑1 (MQ=255)
                                                                                      aaCAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg              >  2:3389506/1‑100 (MQ=255)
                                                                                       aCAAAGAGGAAA‑‑‑GGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa            <  1:3809686/100‑1 (MQ=255)
                                                                                        cAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCTGGGAAGGCCGCCCCGGATCGCGGGCGGCCGAAGACAGCAGAAATTCCAGGTCAGACGGGCTCCAAGTTGCa            >  2:2075430/1‑100 (MQ=255)
                                                                                        cAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa            <  2:2276421/100‑1 (MQ=255)
                                                                                                  aagggggggggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGTCTGCAACTTGCAGCCCGtct    <  2:1510788/98‑3 (MQ=255)
                                                                                                    gggggggggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaa  >  2:245770/1‑100 (MQ=255)
                                                                                                   ||                                                                                                   
TAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAA  >  NC_000913/379914‑380113

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: