Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,998,039 A→G Y124H (TAC→CAC)  fliY ← putative periplasmic binding transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,998,0390AG99.7% 922.8 / ‑5.1 288Y124H (TAC→CAC) fliYputative periplasmic binding transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (164/123);  total (164/124)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

Only 100 of 300 total aligned reads displayed.

CCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACG  >  NC_000913/1997949‑1998128
                                                                                          |                                                                                         
cccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCg                                                                                         >  1:1131189/1‑93 (MQ=255)
cccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCg                                                                                         <  1:2238774/93‑1 (MQ=255)
cccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGCGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATATCAGAAATGGTGTGCg                                                                                         >  2:311818/1‑93 (MQ=255)
 ccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgg                                                                                        <  1:5048367/93‑1 (MQ=255)
  cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                       >  1:8011914/1‑93 (MQ=255)
   ccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgggg                                                                                      >  2:3576974/1‑93 (MQ=255)
    cACTTTTTTGCCTTTCATATCATCGGCTGGTTTAATGGTGTTTTCGTTACCTTTTTTCACTATCGCCTGAATACCATAAATGGTGCTCggggg                                                                                     >  2:8953238/1‑92 (MQ=255)
        ttgttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTgag                                                                                 <  2:9111781/90‑1 (MQ=255)
         ttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTgaga                                                                                >  1:4651441/1‑93 (MQ=255)
          tttGCCGTTCAGATCATCGGTTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGaa                                                                               <  1:4718286/93‑1 (MQ=255)
            tGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAAt                                                                             >  1:1765947/1‑93 (MQ=255)
            tGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAAt                                                                             <  2:2341160/93‑1 (MQ=255)
            tGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAAt                                                                             >  2:5051803/1‑93 (MQ=255)
              ccTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCg                                                                           <  2:2265489/93‑1 (MQ=255)
               cTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGt                                                                          >  2:2967906/1‑93 (MQ=255)
                 ttCAGATCATCGGGTTTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATTCCAGAAAGGTTGTGAGGGGTTGAGAAATCGTAt                                                                        <  1:2816290/93‑1 (MQ=255)
                 ttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAt                                                                        >  2:1859577/1‑93 (MQ=255)
                  tCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAtt                                                                       >  2:7240312/1‑93 (MQ=255)
                     gATCATCGGCTATTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAttttt                                                                    >  2:5860449/1‑93 (MQ=255)
                      atcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAtttttt                                                                   >  1:5488910/1‑93 (MQ=255)
                       tcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTc                                                                  >  2:2109569/1‑93 (MQ=255)
                       tcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTc                                                                  >  2:8294255/1‑93 (MQ=255)
                       tcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTc                                                                  <  1:3354378/93‑1 (MQ=255)
                       tcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGGGCGGGGTTGAGAAATCGGAttttttt                                                                  >  2:4072008/1‑92 (MQ=255)
                        catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCt                                                                 <  1:4405617/93‑1 (MQ=255)
                        catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCt                                                                 >  1:7274061/1‑93 (MQ=255)
                         atcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCtt                                                                <  1:4847128/93‑1 (MQ=255)
                         atcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCtt                                                                <  2:8573131/93‑1 (MQ=255)
                             gCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgcgc                                                            >  1:5249988/1‑93 (MQ=255)
                             gCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgcgc                                                            >  2:5361037/1‑93 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                         >  2:9233484/1‑93 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                         >  2:8138014/1‑93 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                         >  2:6888286/1‑93 (MQ=255)
                                 ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcat                                                        <  2:6612725/93‑1 (MQ=255)
                                  tttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcatc                                                       >  2:6768640/1‑93 (MQ=255)
                                   ttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcatca                                                      <  2:4664744/93‑1 (MQ=255)
                                    tAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAg                                                     <  1:6670785/93‑1 (MQ=255)
                                    tAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAg                                                     <  2:3034334/93‑1 (MQ=255)
                                    tAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAg                                                     <  1:8543959/93‑1 (MQ=255)
                                     aaTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGa                                                    >  2:3368884/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                  >  2:969988/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                  <  1:8703999/93‑1 (MQ=255)
                                         gTGCCTTCGTTACCTTTTTTCACCAGCGACTGAATACCAGAGATGATGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATg                                                <  1:8921976/93‑1 (MQ=255)
                                         gTGCCTTCGTTACCTTTTTTCACAAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATg                                                <  2:5572742/93‑1 (MQ=255)
                                           gCCTTCGTTACCTTTTTTCCCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGt                                              <  1:9281192/93‑1 (MQ=255)
                                             cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGa                                            >  1:1478938/1‑93 (MQ=255)
                                              ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAc                                           >  1:7588597/1‑93 (MQ=255)
                                              ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAc                                           >  1:9835862/1‑93 (MQ=255)
                                              ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAc                                           <  1:8417464/93‑1 (MQ=255)
                                                cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCt                                         >  2:9939579/1‑93 (MQ=255)
                                                cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCt                                         >  2:1704111/1‑93 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                        >  2:1081899/1‑93 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                        >  2:2423266/1‑93 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                        >  1:2577330/1‑93 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                       <  1:8151056/93‑1 (MQ=255)
                                                    aCCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGTGCTCATCAGAAATGGTGACCTGAtt                                     >  1:5752787/1‑93 (MQ=255)
                                                     ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGCGCGGGGTTGAGAAATCGTATTTTTTCTCGCGCTCATCAGAAATGGTGACCTGATTa                                    >  1:7783252/1‑93 (MQ=255)
                                                      cTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTaa                                   >  1:621110/1‑93 (MQ=255)
                                                       ttttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                  <  2:8496051/93‑1 (MQ=255)
                                                         ttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATca                                <  2:1023766/93‑1 (MQ=255)
                                                          tttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcac                               >  1:7013152/1‑93 (MQ=255)
                                                             cACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcaccac                            <  1:2386484/93‑1 (MQ=255)
                                                             cACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcaccac                            >  2:6524591/1‑93 (MQ=255)
                                                              aCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACcaca                           >  2:3251873/1‑93 (MQ=255)
                                                              aCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACcaca                           <  1:6836798/93‑1 (MQ=255)
                                                              aCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACcaca                           >  1:115372/1‑93 (MQ=255)
                                                               ccAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACAt                          >  2:3783576/1‑93 (MQ=255)
                                                               ccAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACAt                          >  1:1892390/1‑93 (MQ=255)
                                                                cAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATc                         >  2:2601277/1‑93 (MQ=255)
                                                                cAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATc                         <  2:5660514/93‑1 (MQ=255)
                                                                 aGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCa                        >  2:7429806/1‑93 (MQ=255)
                                                                 aGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCa                        <  1:3648872/93‑1 (MQ=255)
                                                                  gcgcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa                       >  2:7743221/1‑93 (MQ=255)
                                                                  gcgcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa                       >  2:1228262/1‑93 (MQ=255)
                                                                  gcgcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa                       <  1:3139251/93‑1 (MQ=255)
                                                                   cgcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAAt                      <  1:3177773/93‑1 (MQ=255)
                                                                   cgcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAAt                      >  1:7933791/1‑93 (MQ=255)
                                                                    gcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATa                     <  2:9804941/93‑1 (MQ=255)
                                                                    gcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATa                     >  1:6783204/1‑93 (MQ=255)
                                                                    gcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATa                     <  2:4077778/93‑1 (MQ=255)
                                                                    gcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATa                     <  2:7838748/93‑1 (MQ=255)
                                                                    gcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATa                     <  2:8011913/93‑1 (MQ=255)
                                                                     ccTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATAc                    >  1:5581336/1‑93 (MQ=255)
                                                                     ccTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATAc                    <  2:4405524/93‑1 (MQ=255)
                                                                      cTGAATACCAGAAATGGTGTGCTGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACg                   <  2:8629552/93‑1 (MQ=255)
                                                                         aaTACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATCCGttt                <  2:6228706/93‑1 (MQ=255)
                                                                         aaTACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGttt                <  1:9201973/93‑1 (MQ=255)
                                                                         aaTACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGttt                <  2:3923792/93‑1 (MQ=255)
                                                                         aaTACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACAACATCAATACGttt                >  1:2881361/1‑93 (MQ=255)
                                                                           tACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTag              >  1:64534/1‑93 (MQ=255)
                                                                           tACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTag              <  1:1640770/93‑1 (MQ=255)
                                                                           tACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGATGACCTGATTAATCACCACATCAATACGTTTag              >  1:9726344/1‑93 (MQ=255)
                                                                            aCCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTaga             <  2:3086177/93‑1 (MQ=255)
                                                                             ccAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTagag            <  2:1350469/93‑1 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc          >  1:10029061/1‑93 (MQ=255)
                                                                                 aaaTGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCa        >  1:1407793/1‑93 (MQ=255)
                                                                                 aaaTGGTGTGCGGGGTCGAGAAATCGTATTTGTTCTTGCTCTCATCAGAAATGGTGACCTGATTAAGCACCATATCAATACGTTTAGAGTCCa        <  2:4898850/93‑1 (MQ=255)
                                                                                     ggTGTGCAGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa    <  1:6018741/93‑1 (MQ=255)
                                                                                      gtgtGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGAc   >  1:6868054/1‑93 (MQ=255)
                                                                                      gtgtGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGAc   >  2:4519523/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACG  >  NC_000913/1997949‑1998128

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