Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,136,028 A→C 10.2% V155V (GTA→GTC) ‡ flgH → flagellar protein of basal‑body outer‑membrane L ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,136,0280AC10.2% 218.7 / 13.6 75V155V (GTA→GTC) ‡flgHflagellar protein of basal‑body outer‑membrane L ring
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base C (4/4);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGC  >  NC_000913/1135880‑1136171
                                                                                                                                                    |                                                                                                                                               
ttGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTa                                                                                                                                                 >  8:59248/1‑149 (MQ=255)
   aTACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTg                                                                                                                                              <  2:153805/149‑1 (MQ=255)
   aTACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTg                                                                                                                                              <  7:403812/149‑1 (MQ=255)
     aCTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTgg                                                                                                                                             <  4:290058/148‑1 (MQ=255)
       tgtgCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAGAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTc                                                                                                                                           >  8:228193/1‑148 (MQ=255)
         tgCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGTACCTGGTCaa                                                                                                                                         <  8:409296/148‑1 (MQ=255)
         tgCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAAc                                                                                                                                        <  6:163049/149‑1 (MQ=255)
           ccGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGTACCTGGTCAAcc                                                                                                                                       <  5:8975/148‑2 (MQ=255)
             gcgcTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGc                                                                                                                                     >  6:197788/1‑148 (MQ=255)
               gcTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAAc                                                                                                                                  <  3:264310/149‑1 (MQ=255)
                   tttGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTg                                                                                                                               >  3:17710/1‑148 (MQ=255)
                   tttGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGc                                                                                                                              <  5:197788/149‑1 (MQ=255)
                     tGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCa                                                                                                                             <  1:4716/148‑1 (MQ=255)
                      gCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCAt                                                                                                                            <  6:152776/148‑1 (MQ=255)
                         ggggCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAGAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATgtgg                                                                                                                        <  7:228193/149‑1 (MQ=255)
                           ggCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATgtggtg                                                                                                                      >  7:479277/1‑149 (MQ=255)
                            gCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGtggtgg                                                                                                                     <  4:17710/149‑1 (MQ=255)
                            gCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGtggtgg                                                                                                                     <  1:284930/149‑1 (MQ=255)
                              tGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCCTTAGCGGCACGTTGGCGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTggtgggt                                                                                                                   >  7:46990/1‑149 (MQ=255)
                                tttGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGATTGACCAGGTACTGGACAACGGCAACCTGCATGTcgtgggtg                                                                                                                  >  4:9808/1‑148 (MQ=255)
                                     tAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCa                                                                                                                             >  6:434042/1‑132 (MQ=255)
                                     tAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCa                                                                                                                             <  5:434042/132‑1 (MQ=255)
                                      gaCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACC‑‑TAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGaaaaa                                                                                                             >  2:125275/2‑145 (MQ=255)
                                       aCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGaaaaa                                                                                                             >  3:93025/1‑146 (MQ=255)
                                           tCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTg                                                                                                      >  7:468753/1‑149 (MQ=255)
                                              tGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCa                                                                                                   <  5:223376/149‑1 (MQ=255)
                                               gCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCAt                                                                                                  <  8:468405/149‑1 (MQ=255)
                                                  gATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTaa                                                                                               >  8:71153/1‑149 (MQ=255)
                                                       cGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAgg                                                                                          <  4:438741/149‑1 (MQ=255)
                                                        gAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAgg                                                                                          <  7:59248/148‑1 (MQ=255)
                                                         aaGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAggg                                                                                         >  5:406210/1‑148 (MQ=255)
                                                              tCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCg                                                                                    >  4:305208/1‑148 (MQ=255)
                                                                cGGTGGTAACAGGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAAt                                                                                 <  5:391667/149‑1 (MQ=255)
                                                                   tggtAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTa                                                                              >  4:270418/1‑149 (MQ=255)
                                                                      tAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCCTTAGCGGCACGTTGGCGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTc                                                                           <  8:46990/149‑1 (MQ=255)
                                                                        acacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGc                                                                         >  7:477765/1‑149 (MQ=255)
                                                                         cacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCAt                                                                                                                            >  2:375554/1‑97 (MQ=255)
                                                                         cacGTTCAACCGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCAt                                                                                                                            <  1:375554/97‑1 (MQ=255)
                                                                          acGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCtt                                                                       <  1:70596/149‑1 (MQ=255)
                                                                                aaCGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTgg                                                                 >  4:167697/1‑149 (MQ=255)
                                                                                 aCGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGc                                                                >  5:462487/1‑149 (MQ=255)
                                                                                    gAAAGGGCGGGGCCAATGCCAGCAATACC‑‑TAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGt                                                           <  1:125275/149‑1 (MQ=255)
                                                                                      aaGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGt                                                           <  3:270418/149‑1 (MQ=255)
                                                                                      aaGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGt                                                           <  3:9808/149‑1 (MQ=255)
                                                                                              ggCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCAc                                                   <  8:468753/149‑1 (MQ=255)
                                                                                              ggCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCAc                                                   >  2:269450/1‑149 (MQ=255)
                                                                                                  aaTGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCAc                                               >  3:42255/1‑149 (MQ=255)
                                                                                                         gCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGc                                        <  3:305208/149‑1 (MQ=255)
                                                                                                           aaTGCCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCgg                                      >  8:218104/1‑149 (MQ=255)
                                                                                                           aaTACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCgg                                      <  6:406210/149‑1 (MQ=255)
                                                                                                           aaTACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCCTCTCTGGCGTGGTTAATCCACGCACTATCAGCgg                                      >  2:443445/1‑149 (MQ=255)
                                                                                                            aTACCTTTAGCGGCACGTTGACGGTGATGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCgg                                      >  1:239026/1‑148 (MQ=255)
                                                                                                                cTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGt                                                           >  6:172895/1‑123 (MQ=255)
                                                                                                                cTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGt                                                           <  5:172895/123‑1 (MQ=255)
                                                                                                                cTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagc                                  >  6:458871/1‑148 (MQ=255)
                                                                                                                cTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagc                                  >  1:86904/1‑148 (MQ=255)
                                                                                                                cTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagc                                  >  2:34906/1‑148 (MQ=255)
                                                                                                                 tttAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCaa                                >  2:115307/1‑149 (MQ=255)
                                                                                                                 tttAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCaa                                >  1:252347/1‑149 (MQ=255)
                                                                                                                   tAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCaa                                <  3:167697/147‑1 (MQ=255)
                                                                                                                      cGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccg                           <  1:269450/149‑1 (MQ=255)
                                                                                                                        gCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccgta                         >  7:340302/1‑149 (MQ=255)
                                                                                                                                acggtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTa                  <  7:218104/148‑1 (MQ=255)
                                                                                                                                   gttgccgtTGACCAGTACCTGGTCAACCGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagc                                  <  7:442008/123‑1 (MQ=255)
                                                                                                                                   gttgccgtTGACCAGTACCTGGTCAACCGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagc                                  >  8:442008/7‑129 (MQ=255)
                                                                                                                                       ccgtTGACCAGTACCTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTg          >  5:57239/3‑149 (MQ=255)
                                                                                                                                        ggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTgg         >  8:207489/1‑149 (MQ=255)
                                                                                                                                        cgtTGACCAGTACCTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCa              <  4:255198/143‑1 (MQ=255)
                                                                                                                                        cgtTGACCAGTACCTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCa              >  3:255198/2‑144 (MQ=255)
                                                                                                                                        cgtTGACCAGTACCTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTgg         >  5:481647/2‑149 (MQ=255)
                                                                                                                                         gtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGc        <  8:479277/149‑1 (MQ=255)
                                                                                                                                            gACCAGGTACTGGTCAACGACAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGa     >  1:340435/1‑149 (MQ=255)
                                                                                                                                               cAGGTACTGGTCAACGGCACCCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATgc  >  5:169244/1‑149 (MQ=255)
                                                                                                                                               cAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAg             <  8:24617/138‑1 (MQ=255)
                                                                                                                                               cAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAg             >  7:24617/1‑138 (MQ=255)
                                                                                                                                               cAGGTACGGGTCAACGGCACCCTGCATGGGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCCCGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATgc  >  1:223515/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGC  >  NC_000913/1135880‑1136171

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