Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,319,373 A→T 5.6% D168E (GAT→GAA rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,319,3730AT 5.6% 273.9 / 3.4 84D168E (GAT→GAArcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base T (2/3);  total (42/42)
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

AGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAA  >  NC_000913/2319233‑2319517
                                                                                                                                            |                                                                                                                                                
aGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTc                                                                                                                                          >  4:64619/1‑149 (MQ=255)
  aTAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCa                                                                                                                                        >  5:378521/1‑149 (MQ=255)
    aaaCTGGCGTCGAAGCGTAAAAATACCCGACGCCAGGGCGCGGACGTGGGGTGAGCCAGTAGAGGTTAGTGCCGGTATGAGGGTGGGGTGCATGTCTAGAGTTTGTGATGCTTTCATGCGGACCTTTCAGCGCGGTATCGCGTGCGACt                                                                                                                                      >  8:434056/1‑149 (MQ=255)
    aaaCTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTACGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTTGTGCATTTCGATATTTATTGATGCGTTAATGCAAAGATTTCAAAGCGGTATCGCGTTCCACt                                                                                                                                      >  6:276380/1‑149 (MQ=255)
     aaCTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTg                                                                                                                                     >  5:125940/1‑149 (MQ=255)
        tGGCGTCAAAGCGTAAAAATACCCGACGCCAGGGCGAGGACATTCGCTGATCCAGTAGAGGTTACTGACGCTATCATATTGTGGAGAATTTAGATAGTTATTGATGAGTTCATGCAATGCTTTCAACGCGGTATCGAGTTCAAATggaa                                                                                                                                  >  4:257504/1‑147 (MQ=255)
          gCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTTTCGCGTTCCACTGGCata                                                                                                                                >  8:430396/1‑149 (MQ=255)
          gCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTTTCGCGTTCCACTGGCata                                                                                                                                >  8:436111/1‑149 (MQ=255)
          gCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTgg                                                                                                                                    <  7:159798/145‑1 (MQ=255)
                        aaaTACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGa                                                                                                                   >  6:174113/1‑148 (MQ=255)
                             cccGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTg                                                                                                              <  4:63284/148‑1 (MQ=255)
                              ccGACGCCAGGGCGCGGACCTTCGCTGATCCAGTCGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCa                                                                                                                         <  2:38083/136‑1 (MQ=255)
                              ccGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCGTTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGc                                                                                                            >  8:196491/1‑149 (MQ=255)
                              ccGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCACCGCGGTATCGCGTTCCACTGGCATATCACGCa                                                                                                                         >  1:38083/1‑136 (MQ=255)
                              ccGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGc                                                                                                            >  5:259339/1‑149 (MQ=255)
                              ccGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGc                                                                                                            >  3:186967/1‑149 (MQ=255)
                                       gggCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCTTATCACGCAGACCGAAATTGGCCATGCagag                                                                                                   >  8:334831/1‑149 (MQ=255)
                                       gggCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCagag                                                                                                   <  7:22474/149‑1 (MQ=255)
                                        ggCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAgg                                                                                                  <  8:181022/149‑1 (MQ=255)
                                         gcgcGGACCTTCTCTGATCCATTAGAGGTTACTGCCTCTTTCATCTTGTGGTGCATTTCGATATTTATTGATTCTTTCATTCAAATCTTTCAACGCTTTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAgg                                                                                                  <  3:257504/148‑1 (MQ=255)
                                          cgcgGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGtt                                                                                                <  6:259339/149‑1 (MQ=255)
                                             ggACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATc                                                                                             <  4:186967/149‑1 (MQ=255)
                                                 cTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCTTATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTg                                                                                         <  7:334831/149‑1 (MQ=255)
                                                    cGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGc                                                                                        <  5:132647/147‑1 (MQ=255)
                                                       tGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGAt                                                                                    >  3:113546/1‑148 (MQ=255)
                                                             aGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAAt                                                                             <  6:378521/149‑1 (MQ=255)
                                                                  agGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCg                                                                        >  1:59528/1‑149 (MQ=255)
                                                                  agGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCg                                                                        <  3:144373/149‑1 (MQ=255)
                                                                    gatACTGCCGCGATCATCTTGTGGTGCATTTCGATATGTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGCCCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGt                                                                      >  6:230258/3‑149 (MQ=255)
                                                                    gTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGt                                                                      >  2:22364/1‑149 (MQ=255)
                                                                     ttACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGtt                                                                     <  6:383553/149‑1 (MQ=255)
                                                                        cTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTa                                                                   <  6:24541/148‑1 (MQ=255)
                                                                         tGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGtt                                                                     >  6:77672/1‑145 (MQ=255)
                                                                           ccGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTTTCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGAt                                                                <  7:436111/148‑1 (MQ=255)
                                                                           ccGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTTTCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGAt                                                                <  7:430396/148‑1 (MQ=255)
                                                                           ccGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGtt                                                                     >  5:422053/1‑143 (MQ=255)
                                                                           ccGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGa                                                                 <  5:174113/147‑1 (MQ=255)
                                                                           ccGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATc                                                               <  4:113546/149‑1 (MQ=255)
                                                                                 tcatcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAg                                                          <  6:422053/148‑1 (MQ=255)
                                                                                 tcatcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGc                                                         <  1:401989/149‑1 (MQ=255)
                                                                                    tcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCg                                                       >  1:196973/1‑148 (MQ=255)
                                                                                           gtgCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTTTCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGaaaa                                                <  1:80794/148‑1 (MQ=255)
                                                                                           gtgCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACg                                                                                                                           >  1:157246/1‑73 (MQ=255)
                                                                                           gtgCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACg                                                                                                                           <  2:157246/73‑1 (MQ=255)
                                                                                           gtgCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGaaaaa                                               >  7:319262/1‑149 (MQ=255)
                                                                                              cATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGATGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTa                                            >  8:155152/1‑149 (MQ=255)
                                                                                                tttCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATc                                          <  3:64619/149‑1 (MQ=255)
                                                                                                tttCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATAGTAAGCCGCAGAAAAATTATc                                          >  2:298360/1‑149 (MQ=255)
                                                                                                 ttCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTaa                                                                                 <  4:336308/109‑1 (MQ=255)
                                                                                                 ttCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTaa                                                                                 >  3:336308/1‑109 (MQ=255)
                                                                                                 ttCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATcg                                         >  8:310774/1‑149 (MQ=255)
                                                                                                   cGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATcgc                                        <  5:252855/148‑1 (MQ=255)
                                                                                                      tatttattGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCgcgc                                      >  6:352844/1‑147 (MQ=255)
                                                                                                         ttattGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAg                                   >  5:6326/1‑147 (MQ=255)
                                                                                                              gATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCgcgc                                      >  2:214692/1‑139 (MQ=255)
                                                                                                              gATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCgcgc                                      <  1:214692/139‑1 (MQ=255)
                                                                                                               aTGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCa                            <  6:6326/148‑1 (MQ=255)
                                                                                                                  cGTTGATGCAAAGCTTTCAACGCGGTATCGCGTTCCATTGGGATATCACGCAGACCGAAATTGGCCATTCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCAGCAGAAAAATTATCGCGCCAGTAGCGCATaaa                        <  7:434056/149‑1 (MQ=255)
                                                                                                                   gTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAAc                       >  6:212664/1‑149 (MQ=255)
                                                                                                                    ttCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCACAAAc                       >  1:161595/1‑148 (MQ=255)
                                                                                                                       attcAAAGCTTTCAACGCGGTATCTCGTTCCACTGGCAAATCACGCAGACATAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCTCATAAACCACg                   <  5:276380/146‑1 (MQ=255)
                                                                                                                       aTGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACg                   >  5:149749/1‑149 (MQ=255)
                                                                                                                        tGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCCCGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGAGTAAAAATACCCGGTTTAGAGCGTAAGCCGCATAAAAATTATCGAGCCAGTAGCGCATCAACCACGc                  >  6:285632/1‑149 (MQ=255)
                                                                                                                         gCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGc                  >  7:298629/1‑148 (MQ=255)
                                                                                                                          cAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGGTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGcc                 >  1:377471/1‑148 (MQ=255)
                                                                                                                           aaaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCaa               >  6:232798/1‑149 (MQ=255)
                                                                                                                           aaaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCaa               <  8:434205/149‑1 (MQ=255)
                                                                                                                           aaaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCaa               >  4:204757/1‑149 (MQ=255)
                                                                                                                            aaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAAt              <  1:22364/149‑1 (MQ=255)
                                                                                                                            aaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAAt              <  2:59528/149‑1 (MQ=255)
                                                                                                                            aaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATAGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAAt              <  1:298360/149‑1 (MQ=255)
                                                                                                                            aaGCTGTCAACGCTGTATCGCGTTCCACTTGCATATCACTCAGACCTAAATTGGCAATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAAt              <  5:285632/149‑1 (MQ=255)
                                                                                                                             aGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAAt              >  2:194897/1‑148 (MQ=255)
                                                                                                                                 ttCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAAt              <  5:390216/144‑1 (MQ=255)
                                                                                                                                 ttCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAAt              >  6:390216/1‑144 (MQ=255)
                                                                                                                                  tCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAGATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTc         >  5:85257/1‑148 (MQ=255)
                                                                                                                                    aaCGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAgcgc                             <  2:310272/126‑1 (MQ=255)
                                                                                                                                    aaCGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAgcgc                             >  1:310272/1‑126 (MQ=255)
                                                                                                                                     aCGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACACGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCaga     <  5:77672/149‑1 (MQ=255)
                                                                                                                                      cgcgGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCagag    <  5:232798/149‑1 (MQ=255)
                                                                                                                                        cgGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGATaa  <  5:212664/149‑1 (MQ=255)
                                                                                                                                        cgGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGaa  <  5:352844/149‑1 (MQ=255)
                                                                                                                                        cgGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGaa  <  7:155152/149‑1 (MQ=255)
                                                                                                                                        cgGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGaa  <  3:204757/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
AGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAA  >  NC_000913/2319233‑2319517

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