Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 845,227 A→C 19.3% V85G (GTG→GGG)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913845,2270AC19.3% 106.3 / 2.9 66V85G (GTG→GGG) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base A (19/29);  new base C (10/5);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 8.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATACAGCGGCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCAGAAGAGGATGA  >  NC_000913/844989‑845464
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
tcgaaccAGCGGCAGTGCGCACAGGCGAAACAGCCAGTAAAAACCAAACACTAATCCCGCGAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCa                                                                                                                                                                                                                                                   <  2:298441/233‑1 (MQ=255)
       aGCGGCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTTCACCCCGGTTTGCTCTTCGACCAGCGTCGGCGGCCCTAtttt                                                                                                                                                                                                                                           >  1:1308864/1‑239 (MQ=255)
         cGGCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTTCCTTCTGTGGCCTTTTGCAGGCGGGTTTGCTCTTCGACCCGCGCCGGCGGCCCTAtttttt                                                                                                                                                                                                                                         >  1:1069814/1‑238 (MQ=255)
           gCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGCGCCGGGGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCCTCCTCCGGGGACCCTTTGCGGCCCGGTTTTTCCTTCGCCCCGAGTCGGCGGGCCCTTTTTTTgt                                                                                                                                                                                                                                        >  2:1337251/1‑238 (MQ=255)
              gTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCCACGTCGGCGGGCCTATTTTTGGTAcc                                                                                                                                                                                                                                    >  1:633899/1‑239 (MQ=255)
                gcgcACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTAc                                                                                                                                                                                                                                     <  1:708891/236‑1 (MQ=255)
                gcgcAAAGGGAAACCAGCCAGTAAAAACCAAACACTAATACCGGTAACATGGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGGGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTAc                                                                                                                                                                                                                                     <  2:745123/236‑1 (MQ=255)
                  gaacaGGGGAATCAGCCAGTAAAAACCAAACCATAATGCCGCTAACATTGAAAAATGGGTCAGCGGATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGt                                                                                                                                                                                                                                 <  2:191379/236‑1 (MQ=255)
                          aaaTCAGCCAGTAAAAACCAAACCCTAATACCGCTAACATTGAAAAATGGGGCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCCCTATTTTTGGTACCGGTTCAGCAgg                                                                                                                                                                                                                         <  2:1303240/237‑1 (MQ=255)
                                  gagaaaaagCCAAACAGTAATAACACTCAACTGGAAGAAGGGGTCCGTGCATTGGTAAAGCTGTGTGGGTTAAACGGATTAGGGGTGGAGCCCGTGATATTGCGATAAAGTTGCCCGAAGCGGGCTGAAAGGGCTTCACCAGAATTGCGGCTGACTTCGTTGACCTTTTGCACCACGGTTTACGCTTCGACCAGCGTCGGCGGACCTATTTTTGGTAACGGTTCCGCAGGGGGCGt                                                                                                                                                                                                                   <  2:1083236/232‑1 (MQ=255)
                                       aaaaCCAAACCCTAACTCCGCAAACATGGAAAATGGGGTACGCGCATTGCTGAAGGTTTGTGGTATAAACGGATTATGCGGGGAGCCGGGTATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCaa                                                                                                                                                                                                          <  2:869740/239‑1 (MQ=255)
                                       aaaaCCAAACACTAATTCCGCTAACATTGAAAAATGGGTTAGCGCATTGCTGAAGGTTTGTGGATTAAACCGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCaa                                                                                                                                                                                                          <  2:422755/239‑1 (MQ=255)
                                            cAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGGTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTg                                                                                                                                                                                                     <  2:1192135/239‑1 (MQ=255)
                                               acacTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCGTTTGCAGCACGGGTTTGTCTTCGACCAGCGTCGGCGGCACTATTTTTTGTACCGGTTTCGCCAGGGGCGGAgggggcaagaggacg                                                                                                                                                                                                  >  2:139031/1‑228 (MQ=255)
                                                     aTACCGCTAACATGGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAACCGGCTTATGCGGGGAGCCGGTGATAGTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTgg                                                                                                                                                                                            <  1:389702/239‑1 (MQ=255)
                                                           cTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGGCCTTTTGGAGGACGGTTTGCTCTTCGGCCCCCGGCGGCGGCACTATTTTTGGTACCGGTTCACCAGGGGGCGTAGCGGCAAAGGGGCGCaaacgggcggtc                                                                                                                                                                                      >  2:18377/1‑229 (MQ=255)
                                                           cTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCCCTTCGACCAGCGTCGGCGGCCCTATTTTTGGTACCGGTTCCGCAGGGGGGCTAGCGGCCACGGTTGGCCACTGGGCGGTc                                                                                                                                                                                      >  1:362541/1‑239 (MQ=255)
                                                           cTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGGTTGCTCTTCGACCAGCGTCGGCGGCCCTATTTTTGGTACCGGTTCCGCAGGGGGAGTAGCGGCAACGGTGCGCAACTggtcggtc                                                                                                                                                                                      >  1:34381/1‑239 (MQ=255)
                                                            aaactgtGAAAAATGGGTGGGGCGATTGCTGAAGGTTTGGGGATTAAACGGGTTAGGGGGGGAGCCGGGGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGGGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCa                                                                                                                                                                                     <  2:1080112/233‑1 (MQ=255)
                                                             aaaagtGAAAAAGGGGTCAGGGCATTTCTGAAGGTTGGTGGATTAAACGGCTTGTGGGGGGAGCCGGGGATATGTCGATAAAGGTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGAACTTTTGCAGCACGGGTTGCTCTTCGACCAGCGGCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCaa                                                                                                                                                                                    <  2:1352694/234‑1 (MQ=255)
                                                             aaCATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCCTTTGCAGCACGGTTTGCTCTTCGGACCACGCTGGCGGCCCTTATTTTTGACGCGGGTTAGGCCGGGGCGTAGCGGGAAAGGTGCGACCCTGGTCgcggaa                                                                                                                                                                                    >  2:110333/1‑234 (MQ=255)
                                                             aaCATGTAAAAAGTGGGCAGCGCAGTGGTGTAGGGTTGTGGATTAAACGGGTTTGGGGGGGGGGCGGTGAAATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCaa                                                                                                                                                                                    <  2:990284/239‑1 (MQ=255)
                                                                 ttGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCCCGGTTTGCTCTTCGACCAGGGTCTGCGGCCCTATTTTTGGTACCGGTTCAGCAGGGGGGGCAGCGGCCACGGGGGGCAACCGGTGGGTCaaacct                                                                                                                                                                                >  2:420347/1‑235 (MQ=255)
                                                                 ttGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAACGGCAACCGTGGGCAACTGGTCCATCAActct                                                                                                                                                                                >  1:526452/1‑239 (MQ=255)
                                                                      aattggggtcgcgcATTGCTGAAGTTTTGTGGATTAAACGGCTTAGGGGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTAcg                                                                                                                                                                          <  2:673671/229‑1 (MQ=255)
                                                                      aaaTGGGTCGGGGCTTTGTTGAAGGTTTGGGGAGTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGCCTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTGTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTAcg                                                                                                                                                                          <  2:564139/239‑1 (MQ=255)
                                                                          gggTCAGGCCATTGGTGGAGTGTTTGGGATTAAACTGCTTTGGGGGGGGGCGGGTGATATTGCGAAAAAGGTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAg                                                                                                                                                                      <  2:891517/239‑1 (MQ=255)
                                                                           ggTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAg                                                                                                                                                                      <  1:1337251/238‑1 (MQ=255)
                                                                                           aaGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGTCTTCGGTGACCTTTTTAAGCCCGGTTTGCTCTTCGACCGCCGCCGGGGGCGCTTTTTTTTGTTCCGGGTCATCAGGGGGCGTGGCGGGCAACGTGCACACCTGGGCGGTCAACTCTTTACGCGAGGGGGCAATATCCaagc                                                                                                                                                      >  2:565952/1‑235 (MQ=255)
                                                                                           aaGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTAACGGTTCAGCAGGGGGCGTAGGGGCAACGGTGCGCAACTGGTCGATCAACCCTTTACGCGAGGGGGCATTTTCCAGcac                                                                                                                                                     >  1:1360511/1‑239 (MQ=255)
                                                                                                 ttGTGGATTAAACGGGTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGGGGGCAACTGGTCGATCAACTCTTTACCCCAGGGGTCATTATCCAGCAAATcc                                                                                                                                                 >  1:801711/1‑238 (MQ=255)
                                                                                                  ttttgATTAAACCGCTTAGGGGGGGAGCCGGTTAATTTGGGATAAATGTGCCCGAAGCGGGCGGAAAGGGCTTCCCCATAATGGCGGCTGCCTTCGGTGACCTTTGTCAGCCCGGTTTGCTCTTCGACCAGCGTCGGCTGCCCTATTTGTGGTCCCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGc                                                                                                                                              <  2:1376391/235‑1 (MQ=255)
                                                                                                    tgGATTAAACGGTGTAGGGGGGGAGCCGGTGATATTCCGATAAAGTTGCCCGAAGCGGGGGGAAAGGGCTTCACCATAATGGCGGCTGACTCCGGTGACTTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCa                                                                                                                                            <  2:1069814/239‑1 (MQ=255)
                                                                                                         tAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCa                                                                                                                                                                                                                                                   <  2:917256/135‑1 (MQ=255)
                                                                                                         tAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCa                                                                                                                                                                                                                                                   >  1:917256/1‑135 (MQ=255)
                                                                                                         tAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGGAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCCGCACATCCGCCAGGgc                                                                                                                                        >  1:260634/1‑238 (MQ=255)
                                                                                                          aaaCGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCa                                                                                                                                                         <  1:28198/221‑1 (MQ=255)
                                                                                                                          gAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGGCTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGGCCAGCGGCCGCGGCACTATTTTTGGGACCGGGTCAGCAGGGGGCGTAGCGGCGACGGTGGGCCACCGGTCGGTCCACCCTTTACCCGAGGGGGCATTATCCAACACATCCGCCCCGTGCCCCACAGCCCCCttttt                                                                                                                       >  2:1206264/1‑239 (MQ=255)
                                                                                                                                     atTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGGTCAACTCTTTACGCGAGGTGCCATTATCCAGCACATCCGCCAGTGCGCCTTAGGCCGCTTTTTTTTGTTCGata                                                                                                            >  1:655023/1‑238 (MQ=255)
                                                                                                                                         cGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCaa                                                                                                                                                                                    <  1:703717/164‑1 (MQ=255)
                                                                                                                                           aTAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGGTCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGt                                                                                                      >  1:859858/1‑238 (MQ=255)
                                                                                                                                                    cccGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGCTCAGCAAACCCTTACCCAGGGGGGCCTTATCCAGCCCATCCGCCCAGGCGCCCTAAGCCCCTTTTTTTTTTTTGATATCCGGTTCCCGGGGCg                                                                                            >  2:486545/1‑239 (MQ=255)
                                                                                                                                                       gaggcgggcggAAAGGGCTTCACCATATTGGCGGGTGACTTGGGTGCCCTTTGGCAGACGGGTTTGCTCTTCGACCGGCGTGGGCGGCACTATTTTGGGTACCGGTTAAGCAGGGGGCGTAGGGGCAACGGTGCGCAACTGGTCGATCACCTCTTTACGGGGGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGtt                                                                                         <  2:1280161/236‑1 (MQ=255)
                                                                                                                                                          ggggggggAGAGGTTTTCCCCATAATGGGGGGTGGTTTGGGTGACCTTTTGCAGCAGGGTTGCCTCTGCGACCAGCGTCGGGGGCTATTTTTGTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACGGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGtt                                                                                         <  2:260634/234‑1 (MQ=255)
                                                                                                                                                                tgAAAGGGTTTCCCCATAAGGGGGGGGGACCTGGGGGACCTTTTGCAGCGGGGTTTGCTCTCTGACCCGCGTGGGGGGCACTATTTTTGGTACCGGTTCAGCGGGGGGCGTAGGGGCAACGGTGCGCCACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGtt                                                                                 <  2:1015845/237‑1 (MQ=255)
                                                                                                                                                                ggAAGGGGTTTCCCCTTATGGGGGGCTGTCTTGGGTGACCTTTTGCCGCCCGTTTTGCTTTTGGACCAGCGTCGGCGGCACTTTTTTTGGTCCCGTTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCACCTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTg                                                                                <  2:34381/239‑1 (MQ=255)
                                                                                                                                                                        cTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCCTAGGCCGCTTTTTTTTGTTCGGTATCCGGTTCCGGGGCCgggt                                                                                         >  1:776050/1‑220 (MQ=255)
                                                                                                                                                                        cTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAg                                                                                      >  1:274492/1‑225 (MQ=255)
                                                                                                                                                                        cTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCcgt                                                                         <  1:893906/238‑1 (MQ=255)
                                                                                                                                                                        cTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGGGTCATTATCCAGCACATCCGCCAGAGCGCCATAGGCCGCTTTTTTTTGTTCGGTATCCGGTTCCGGGGCCGGTTCAGTCGTTGGGTCcgtt                                                                        >  1:1127342/1‑239 (MQ=255)
                                                                                                                                                                          tcccaaaaaTGGGGGCTGACTTCGGTGCCCTTTTGCAGCCCGGTTTGCTCTTCGACCAGGGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTCGAGTCcgttgt                                                                      <  1:1052195/232‑1 (MQ=255)
                                                                                                                                                                                tAATGGCGGCTGACTTCGGTGACCTTTTGGAGGCCGGTTTGCCCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCcgttgtcgttgt                                                                <  1:18377/239‑1 (MQ=255)
                                                                                                                                                                                tAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCCACGGTGCGCCACTGGTCGGTCAACTCTTTACGCGAGGGGGCATTTTCCAGCACCTCCGCCCGGGCGCCCTAGGCCGCTTTTTTTTTTTCGCTATCCGGTTCCGGGGCCGGTTCGGCCGTTGGGTCcgttttggttgt                                                                >  2:893906/1‑239 (MQ=255)
                                                                                                                                                                                tAATGGCGGCCGGCTTCGGTGAGCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGGCCGCGGCCCTATTTTTGGTACCGGTTCAGCTGGGGGCGTAGCGGCAACGGCGCGCAAATGGTCGATACACTTCTTACGCGAGGGGTCATGATCCAGCACACACGCCCATGCGCCCTAGGGCGCGTTTTTTTGTTCGCAATCCGGTTCCGGGGGCGGGTCAGTGGTTGGGTGcgttgtcgttgg                                                                >  2:1303833/1‑238 (MQ=255)
                                                                                                                                                                                        gcTGACTTCGGTGACCTTTTGCAGCACGGTTTTCTCTTCGACCAGCGTCGGCGGCCCTATTTTTTGTACCGGTTCAGCCGGGGGCGTAGCCGCCACGGTGCGCACCTGGTCGTTCCACTCTTTACTGCGGGGGTCCTTATCCAGCCCCTCCGCCCGTCCGCCCTGGGCCCCTTTTTTTTTTTCCATATTCGGTTCCGGGGCCGGTTTCGGCGGTTAGTCCGTTTTTGTTGTgggttgca                                                        >  2:653988/1‑233 (MQ=255)
                                                                                                                                                                                        gcTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCAGTTCAGCAGGGGGCGGAGCGGCAACGGGGCGCAACTGGTCGATCAACTCTTTACGCGAGGGGTCATTATCCAGCACATCCGCCAGGGCGCCCTAGGCCGCCTTTTTTTGTTCGTATTCCGGTTCCTGGGCCGGTTACGTGTTTGAGTCCCTTGTGTTggtggtgggta                                                        >  2:553922/1‑234 (MQ=255)
                                                                                                                                                                                           tgtttCGGTGACTTTTTGCTGCACGTTTTGCTCGTCGGCCAGGGTCGGCGGCCCTTTTTTTGGTACCGGTTCAGCGGGGGGCGTGGCGGCAACGGTGCGCAACTGGTCGATCAACTTTTTACGCGAGGTGTCATTATCCAGCCCATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGc                                                     <  2:1308864/236‑1 (MQ=255)
                                                                                                                                                                                             cTTCGGTTACCTTTTGCAGCACGGTTTGCTCTTCGGCCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCg                                                   <  1:486545/239‑1 (MQ=255)
                                                                                                                                                                                                       ctttgtcagccCGGTTTGCTTTTCGTCCAGGGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGAt                                          <  1:17455/232‑1 (MQ=255)
                                                                                                                                                                                                           tGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACcgc                                     <  1:420347/239‑1 (MQ=255)
                                                                                                                                                                                                                cccttttttCTCTTCGCCCGGCGGCGGCGGCACTATTGTGGGTCCCGTTTCGGCAGGCGGTGGGGCGGCAACGGTGCGAAAAGGGCCGACCAACCCTGTCCGCGAGGGTCCATTTCCCAGCCCACCCGCCCGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTGCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTggg                                <  2:691641/234‑1 (MQ=255)
                                                                                                                                                                                                                 aCGGTTTGCTCTTCGACCAGCGTCGGCGGCAATATTGTTGGGACCGGGTCAGCAGGGGGCGTAGCGGCAACGGGGCGCAACTGGTCGATCAACTCTTTACGCGAGGGGTCATTGTCCAGCACATCCGCCAGTGAGCCACAGGCCGCATTTTTTTGTTCGATATCCGGTTCCGGGGCCGGGTCAGTCGTTATGTCCGTTGTGGTTGTGGGTGCAACCCCGGGTTGGGATACGGGGTgggg                               >  2:982364/1‑238 (MQ=255)
                                                                                                                                                                                                                     tttGCTTTCCGGCCCGCGTGGGGGGCCGTATTTTTGGTCCCGTTTCGGGAGGGGGCGTGGCGGAACCGGTGCGCACCTGGTGGATCAACTCTTTCCGCGGGGTGTTATTACCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGt                           <  2:1360511/239‑1 (MQ=255)
                                                                                                                                                                                                                             tCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTTTCTTTGTGGTTTTACAGCCGGGTTTGGTAACCCCTGGGGCAGGGGCGccccg                   >  2:800123/1‑237 (MQ=255)
                                                                                                                                                                                                                             tCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTTCCGTTGTGGTTTTCAAGCCGGGTATGGGTACCGCGTGGGCAGGGGGGccccc                   >  2:748295/1‑237 (MQ=255)
                                                                                                                                                                                                                                gACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGAGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGCGTCCGTTGTCGTTGTGGTTTTAACCCCGGGTGTGGATACCGCGCGGGCTGGTGCGCCcagca                 >  1:869151/1‑239 (MQ=255)
                                                                                                                                                                                                                                      cGTGGGCGCCATTTTTTTTGGTCCCGTTTACGCGGGGGGCGTAGCGGCAACGGTGCGCACCTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCaga           <  2:859858/239‑1 (MQ=255)
                                                                                                                                                                                                                                       gTCGGGGGCACTATTTTTGGTCCCGGTTCAGCGGGGGGAGTAGGGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGGGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCaga           <  1:553922/238‑1 (MQ=255)
                                                                                                                                                                                                                                         cggcggCTTTTTTTTTGGTCCCTGTTCGGCAGGGGGCGTAGGGGCAACGGTGCGAAACTGGTCAATCACCTCTTTACGCGAGGTGTCATTATCCAGCACTGCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCagaag         <  2:1127342/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 aCTTTTTTGGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACCCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCAGAAGAGGATGa  <  1:863415/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 aCTGTTTTTGGTGCGGTTTCAGCGGGGGGGGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTAGCCAGCCCATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCAGAAGAGGATGa  <  1:653988/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
GCGATACAGCGGCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAGGTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTCGACCAGCGTCGGCGGCACTATTTTTGGTACCGGTTCAGCAGGGGGCGTAGCGGCAACGGTGCGCAACTGGTCGATCAACTCTTTACGCGAGGTGTCATTATCCAGCACATCCGCCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCAGAAGAGGATGA  >  NC_000913/844989‑845464

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: