Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,395 C→G 5.3% intergenic (+60/+21) pheP → / ← ybdG phenylalanine transporter/mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913603,3950CG 5.3% 679.2 / 20.6 192intergenic (+60/+21)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (88/93);  new base G (3/8);  total (91/101)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

Only 100 of 195 total aligned reads displayed.

GTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATATTCCAGCCACACCACGGTGTTAGTAAACGCGTAGATCTCGAGCGGTAAACCGTTATCACCTGGAGCCAGTTGGCGTACCATTAAGGTCATATCTTTACGAATCCGCGGATGGTTACGCAGATATTCG  >  NC_000913/603130‑603668
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                 
                                                    gttctctcCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                                                                                                                                                                                                      <  2:195428/273‑1 (MQ=255)
                                                     ctctATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                                                                                                                                                                                                                >  2:277425/1‑232 (MQ=255)
                                                         aTCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTa                                                                                                                                                                                                                                                            <  4:277425/232‑1 (MQ=255)
                                                               tCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAc                                                                                                                                                                                                           >  2:176176/1‑275 (MQ=255)
                                                                                    ccTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGCAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAACTGCACCCGCCAGAGAGCGAATCTCATTCCCGGTTGGCGACTGATGAAATCGCAGACCAAAc                                                                                                                                                                                      >  2:210208/1‑275 (MQ=255)
                                                                                                   ttttGCTGCTGAGGTGTACGATGGGTGATATGCGCTTGTCATCGATTCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGTGATTTAAAGCGCTGCGTCGGAAATAAGGCAGTCACGCTACATCCGACAAATCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTGTTTATTACTGCTTAAATGCACCCCGAGGGAGCGAATATCATTGACGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTACGACAGTCGCaa                                                                                                                                                                      <  4:117176/275‑1 (MQ=255)
                                                                                                            tgatttgagcGATGGATGCTATGCGCTTGTCAGTGATCCTGCTTCCGGTGTTGATTGTATTCCTGTTTATGGCATTTAAAACGATGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATAtgtg                                                                                                                                                               <  1:4210/265‑1 (MQ=255)
                                                                                                                     cGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                                                                                                                                                                                                   >  8:124801/1‑149 (MQ=255)
                                                                                                                      gatggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg                                                                                                                                                                                                                                                                                  >  8:878158/1‑149 (MQ=255)
                                                                                                                       atggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa                                                                                                                                                                                                                                                                                 <  5:1286629/149‑1 (MQ=255)
                                                                                                                          gatgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                                                                                                                                                                                              >  8:461916/1‑149 (MQ=255)
                                                                                                                          gatgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                                                                                                                                                                                              <  5:1059991/149‑1 (MQ=255)
                                                                                                                            tgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgaga                                                                                                                                                                                                                                                                            <  5:776142/149‑1 (MQ=255)
                                                                                                                            tgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgaga                                                                                                                                                                                                                                                                            >  7:407568/1‑149 (MQ=255)
                                                                                                                             gatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGagag                                                                                                                                                                                                                                                                           >  8:732649/1‑149 (MQ=255)
                                                                                                                              atatGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                                                                                                                                                                                                          <  6:1102104/149‑1 (MQ=255)
                                                                                                                              atatGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                                                                                                                                                                                                          <  7:778523/149‑1 (MQ=255)
                                                                                                                                     cTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttt                                                                                                                                                                                                                                                                   <  6:504455/149‑1 (MQ=255)
                                                                                                                                       tGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                                                                                                                                                                                                                 <  8:451433/149‑1 (MQ=255)
                                                                                                                                       tGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                                                                                                                                                                                                                 <  5:526720/149‑1 (MQ=255)
                                                                                                                                        gTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                                                                                                                                                                                                                <  8:391473/149‑1 (MQ=255)
                                                                                                                                             cGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTAc                                                                                                                                                                                                                                                           <  7:1334137/149‑1 (MQ=255)
                                                                                                                                                tCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGc                                                                                                                                                                                                                                                        <  8:394525/149‑1 (MQ=255)
                                                                                                                                                      tgcCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                                                                                                                                                                                                  <  8:331036/149‑1 (MQ=255)
                                                                                                                                                      tgcCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                                                                                                                                                                                                  <  5:1206320/149‑1 (MQ=255)
                                                                                                                                                              tgGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGc                                                                                                                                                                                                                                          <  7:384524/149‑1 (MQ=255)
                                                                                                                                                                gATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCa                                                                                                                                                                                                                                        >  5:631972/1‑149 (MQ=255)
                                                                                                                                                                 aTTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCag                                                                                                                                                                                                                                       >  6:752647/1‑149 (MQ=255)
                                                                                                                                                                  ttGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCaga                                                                                                                                                                                                                                      >  7:727275/1‑149 (MQ=255)
                                                                                                                                                                      aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGTGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGc                                                                                                                                                                                                                                  >  5:51423/1‑149 (MQ=255)
                                                                                                                                                                           tGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                                                                                                                                                                                                             >  5:1437710/1‑149 (MQ=255)
                                                                                                                                                                           tGTTAATGGCATTTAAAACGCTGCGTCGGAAATAAGGCACTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                                                                                                                                                                                                             <  8:1229350/149‑1 (MQ=255)
                                                                                                                                                                             tttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                                                                                                                                                                                                >  8:635197/1‑96 (MQ=255)
                                                                                                                                                                             tttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                                                                                                                                                                                                           >  6:1037836/1‑149 (MQ=255)
                                                                                                                                                                               tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATAAATCAGCCTGAATGCTTTCATATTCCAGCCACACCACGGGTTTAGTAAACGCGTAGAttc                                                                                           >  4:215743/1‑273 (MQ=255)
                                                                                                                                                                                   gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                                                                                                                                                                                                     <  7:466264/149‑1 (MQ=255)
                                                                                                                                                                                   gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                                                                                                                                                                                                     >  8:309696/1‑149 (MQ=255)
                                                                                                                                                                                     aTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGcc                                                                                                                                                                                                                                         >  7:314876/1‑127 (MQ=255)
                                                                                                                                                                                      tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGt                                                                                                                                                                                                                  >  8:139717/1‑149 (MQ=255)
                                                                                                                                                                                      tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGt                                                                                                                                                                                                                  >  8:421953/1‑149 (MQ=255)
                                                                                                                                                                                       ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                                                                                                                                                                                                 <  8:581001/149‑1 (MQ=255)
                                                                                                                                                                                       ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                                                                                                                                                                                                 <  7:922182/149‑1 (MQ=255)
                                                                                                                                                                                       ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                                                                                                                                                                                                 <  6:1437710/149‑1 (MQ=255)
                                                                                                                                                                                       ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                                                                                                                                                                                                 <  6:631972/149‑1 (MQ=255)
                                                                                                                                                                                          aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                                                                                                                                                                                                   >  6:558275/1‑112 (MQ=255)
                                                                                                                                                                                          aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                                                                                                                                                                                                              <  7:211364/149‑1 (MQ=255)
                                                                                                                                                                                                 gCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATATTCCAGCCACACCACGGTGTTAGTAAACGCGTAGATCTCGAGCGGTAAACCGTTAt                                                                         >  2:236981/1‑275 (MQ=255)
                                                                                                                                                                                                    tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                                                                                                                                                                                                    >  6:1174384/1‑149 (MQ=255)
                                                                                                                                                                                                    tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                                                                                                                                                                                                    <  7:958088/149‑1 (MQ=255)
                                                                                                                                                                                                    tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCGCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                                                                                                                                                                                                    <  8:1434135/149‑1 (MQ=255)
                                                                                                                                                                                                     cGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAg                                                                                                                                                                                                   <  7:568144/149‑1 (MQ=255)
                                                                                                                                                                                                         aaTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                                                                                                                                                                                               >  7:129783/1‑149 (MQ=255)
                                                                                                                                                                                                          aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                                                                                                                                                                                              <  8:1238518/149‑1 (MQ=255)
                                                                                                                                                                                                            aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                                                                                                                                                                                                <  5:982350/65‑1 (MQ=255)
                                                                                                                                                                                                            aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                                                                                                                                                                                                >  6:982350/1‑65 (MQ=255)
                                                                                                                                                                                                            aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa                                                                                                                                                                                            <  6:906400/149‑1 (MQ=255)
                                                                                                                                                                                                                cATTCATGCTACATCCGAAAAAACGAAGTCAACCTTCCGAAACCACTCTCATCCTGTGGATGAGGGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATATTCCAGCCACACCACGGTGTTAGTAAACGCGTAGATCTCGAGCGGTAAACCGTTATCACCTGGGGCCAGt                                                           <  4:183680/274‑1 (MQ=255)
                                                                                                                                                                                                                 aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCCGACCaaa                                                                                                                                                                                       <  6:806072/149‑1 (MQ=255)
                                                                                                                                                                                                                 aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa                                                                                                                                                                                       >  7:363385/1‑149 (MQ=255)
                                                                                                                                                                                                                       gCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc                                                                                                                                                                                 <  7:271298/149‑1 (MQ=255)
                                                                                                                                                                                                                       gCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc                                                                                                                                                                                 <  7:444023/149‑1 (MQ=255)
                                                                                                                                                                                                                       gCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc                                                                                                                                                                                 <  6:458563/149‑1 (MQ=255)
                                                                                                                                                                                                                            aTCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACaa                                                                                                                                                                            >  5:1313230/1‑149 (MQ=255)
                                                                                                                                                                                                                              ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATc                                                                                                                                                                          <  5:767739/149‑1 (MQ=255)
                                                                                                                                                                                                                              ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATc                                                                                                                                                                          >  5:597698/1‑149 (MQ=255)
                                                                                                                                                                                                                                gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAtt                                                                                                                                                                                                                        <  5:691813/101‑1 (MQ=255)
                                                                                                                                                                                                                                gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAtt                                                                                                                                                                                                                        >  6:691813/1‑101 (MQ=255)
                                                                                                                                                                                                                                 aCAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAatat                                                                                                                                                                                                                            >  7:1077262/1‑96 (MQ=255)
                                                                                                                                                                                                                                 aCAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATATTCCCGCCCCACACCGTTGTTAGTAAACGCGTCGATCTCACGCGGTAAACCGTTATCCCCTGGAGCCAGTTGGCGTACCATTAAGGTc                                         >  4:4967/1‑275 (MQ=255)
                                                                                                                                                                                                                                 aCAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCCCCAGAGAGCGAatat                                                                                                                                                                                                                            <  8:1077262/96‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCaaa                                                                                                                                                                     >  7:648639/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCaaa                                                                                                                                                                     >  7:74686/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCaaa                                                                                                                                                                     >  7:539456/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCaaa                                                                                                                                                                     >  8:167325/1‑149 (MQ=255)
                                                                                                                                                                                                                                    aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATACACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAat                                                                                                                                                                    >  8:1027465/1‑149 (MQ=255)
                                                                                                                                                                                                                                     aaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAata                                                                                                                                                                   >  5:259561/1‑149 (MQ=255)
                                                                                                                                                                                                                                     aaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATATTCCAGCCACACCACGGTGTTAGTAAACGCGTAGATCTCGAGCGGTAAACCGTTATCACCTGGAGCCAGTTGGCGTACCATTAAGGTCat                                       <  2:4967/273‑1 (MQ=255)
                                                                                                                                                                                                                                      aCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAatat                                                                                                                                                                  <  7:88000/149‑1 (MQ=255)
                                                                                                                                                                                                                                      aCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAatat                                                                                                                                                                  >  7:1149118/1‑149 (MQ=255)
                                                                                                                                                                                                                                          tGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGa                                                                                                                                                              >  5:458038/1‑149 (MQ=255)
                                                                                                                                                                                                                                            tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATc                                                                                                                                                            <  8:129783/149‑1 (MQ=255)
                                                                                                                                                                                                                                             cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCg                                                                                                                                                           >  6:746920/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGa                                                                                                                                                          <  5:820822/149‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGa                                                                                                                                                          >  5:33280/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGa                                                                                                                                                          >  8:75934/1‑149 (MQ=255)
                                                                                                                                                                                                                                                ccATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCa                                                                                                                                                                       <  7:616879/134‑1 (MQ=255)
                                                                                                                                                                                                                                                 cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAata                                                                                                                                                       >  6:292922/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAata                                                                                                                                                       <  5:1014054/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   tCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATAtgt                                                                                                                                                                <  5:740364/138‑1 (MQ=255)
                                                                                                                                                                                                                                                     cGAAACCGCTCTCATCCATTGGATGAGAGCGGATTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTACCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATc                                                                                                                                                   >  8:956458/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        aaCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGc                                                                                                                                                >  6:1276546/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTg                                                                                                                                             >  5:431110/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTg                                                                                                                                             >  7:739504/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCTAATATATCAGCCTGa                                                                                                                                            >  6:1468409/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              tctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAAt                                                                                                                                          >  8:649435/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTc                                                                                                                                    <  4:320634/152‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCtt                                                                                                                                      >  8:427359/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCttt                                                                                                                                     >  6:709166/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        tCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTg                                                                                                                                             >  2:366947/1‑136 (MQ=255)
                                                                                                                                                                                                                                                                        tCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGACTGAATGCTTTCATATTCCAGCCACACCACGGTGTTAGTAAACGCGTAGATCTCGAGCGGTAAACCGTTATCACCTGGAGCCAGTTGGCGTACCATTAAGGTCATCTCTTTACGAATCCGCGGATGGTTACGCAGATATTCg  >  2:112542/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                 
GTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATATTCCAGCCACACCACGGTGTTAGTAAACGCGTAGATCTCGAGCGGTAAACCGTTATCACCTGGAGCCAGTTGGCGTACCATTAAGGTCATATCTTTACGAATCCGCGGATGGTTACGCAGATATTCG  >  NC_000913/603130‑603668

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: