Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,139,225 T→C 5.1% I283T (ATT→ACT)  flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,139,2250TC 5.1% 250.5 / 2.9 80I283T (ATT→ACT) flgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base C (2/2);  total (38/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 = 8.74e-01

CACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGT  >  NC_000913/1139078‑1139370
                                                                                                                                                   |                                                                                                                                                 
cACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCAt                                                                                                                                                   >  4:114235/1‑148 (MQ=255)
cACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCAt                                                                                                                                                   >  1:237037/1‑148 (MQ=255)
cACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCAt                                                                                                                                                   >  5:284674/1‑148 (MQ=255)
 aCTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCAt                                                                                                                                                   <  2:237037/147‑1 (MQ=255)
  cTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTc                                                                                                                                                 >  6:214033/1‑148 (MQ=255)
  cTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTc                                                                                                                                                 >  4:288745/1‑148 (MQ=255)
      ttCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGCTCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGAc                                                                                                                                             <  1:8230/148‑1 (MQ=255)
       tCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACa                                                                                                                                            >  3:6241/1‑148 (MQ=255)
          gggAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCTATATTGAGATCCCGGAGATATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTc                                                                                                                                         >  4:52312/1‑148 (MQ=255)
           ggAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTcc                                                                                                                                        >  4:358593/1‑148 (MQ=255)
            gAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCg                                                                                                                                       >  4:347406/1‑148 (MQ=255)
                tACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACCTTCCGTtct                                                                                                                                   >  4:204174/1‑148 (MQ=255)
                  cGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTtctc                                                                                                                                  <  3:288745/147‑1 (MQ=255)
                  cGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTtctc                                                                                                                                  <  7:173102/147‑1 (MQ=255)
                  cGGCGCGGCAACTGGCGGCAGTGCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTtctc                                                                                                                                  <  4:6241/147‑1 (MQ=255)
                    gcgcGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAg                                                                                                                                >  4:109276/1‑147 (MQ=255)
                       cgGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATc                                                                                                                            >  4:257782/1‑148 (MQ=255)
                         gCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTg                                                                                                                          <  2:18866/148‑1 (MQ=255)
                           aaCTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGa                                                                                                                        <  3:114235/148‑1 (MQ=255)
                           aaCTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGa                                                                                                                        <  5:214033/148‑1 (MQ=255)
                             cTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGAcc                                                                                                                      <  1:121622/148‑1 (MQ=255)
                                   gcAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGAcgc                                                                                                                <  6:382684/148‑1 (MQ=255)
                                        tccttccAGCGCTGACCCTTCTCGTAAGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCACGGAGAAAATACTGAATACCGGGTCGCTGGGCGGGAATCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTaa                                                                                                            >  8:55943/1‑147 (MQ=255)
                                         ccttccAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATa                                                                                                          >  7:395508/1‑148 (MQ=255)
                                          cttccAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATAc                                                                                                         >  4:363478/1‑148 (MQ=255)
                                             ccAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGc                                                                                                       <  6:297050/147‑1 (MQ=255)
                                              cAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCt                                                                                                      <  3:109276/147‑1 (MQ=255)
                                                 cgcTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGa                                                                                                  <  3:347406/148‑1 (MQ=255)
                                                   cTGACCCTTCTCGTCCGACTGTCGCTTATGTTGATGGGACGGCAGGCATTATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACa                                                                                                <  5:6445/148‑1 (MQ=255)
                                                   cTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACa                                                                                                <  6:362594/148‑1 (MQ=255)
                                                                    aCTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTg                                                                               <  1:394147/148‑1 (MQ=255)
                                                                       gTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAAc                                                                                              >  4:221929/1‑130 (MQ=255)
                                                                       gTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAAc                                                                                              <  3:221929/130‑1 (MQ=255)
                                                                        tCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGcc                                                                             <  3:204174/146‑1 (MQ=255)
                                                                         cGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGCAGAAATTACTGAATACCGGGTCGCTGGGCGGCCCTCGGCCATACCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTGGGACACCTGGCGCTGGCATTGGCCgcg                                                                          >  4:77757/1‑146 (MQ=255)
                                                                         cGCTTATGTTGATGGGACGGCAGGCAATATTGAGATACCGGAGAAATTAATGAATACCGGGTCGCTGTGCGGCACTCTGACATTCCGTTATCAGGATCTGGAACAGAAGCGTAATACGCTTGGAAAACTGGCGATGGCATTTGCCGAg                                                                          >  4:72808/1‑148 (MQ=255)
                                                                              aTGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTa                                                                                                             >  2:221940/1‑108 (MQ=255)
                                                                              aTGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTa                                                                                                             <  1:221940/108‑1 (MQ=255)
                                                                              aTGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCtt                                                                      <  7:364994/147‑1 (MQ=255)
                                                                                gTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTc                                                                    >  6:268595/1‑147 (MQ=255)
                                                                                 ttGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCaa                                                                  >  2:58681/1‑148 (MQ=255)
                                                                                 ttGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCaa                                                                  <  3:358593/148‑1 (MQ=255)
                                                                                  tGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTgg                                                                                           >  1:403213/1‑122 (MQ=255)
                                                                                  tGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTgg                                                                                           <  2:403213/122‑1 (MQ=255)
                                                                                  tGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCttt                                                                     >  3:291417/1‑144 (MQ=255)
                                                                                  tGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCttt                                                                     <  4:291417/144‑1 (MQ=255)
                                                                                         aCGGCAGGCAATATTGAGATACCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGTCTATCAACACCCaa                                                           >  6:46576/1‑147 (MQ=255)
                                                                                            gcaggcaATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAAcaca                                                       <  6:284674/148‑1 (MQ=255)
                                                                                            gcaggcaATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGAATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGACGAGGCTTTCAACACCCAacac                                                        >  1:61098/1‑147 (MQ=255)
                                                                                                 caATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACCCCCAACACAAAGcc                                                  >  1:44775/1‑148 (MQ=255)
                                                                                                 caATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGcc                                                  <  1:58681/148‑1 (MQ=255)
                                                                                                  aaTATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGAcc                                                                                                                      >  2:294173/1‑79 (MQ=255)
                                                                                                  aaTATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGAcc                                                                                                                      <  1:294173/79‑1 (MQ=255)
                                                                                                    tatTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGa                                               >  5:2646/1‑148 (MQ=255)
                                                                                                     atTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGAt                                              >  4:140728/1‑148 (MQ=255)
                                                                                                     atTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATATGCCGAGGCTTTCAACACCAAACACAAAGCAGGa                                               >  2:191932/1‑147 (MQ=255)
                                                                                                     atTGACATACCGGAGAAAATACAGAAAACCGTGTCTCTGTGCGGCATTCTGACATTCCGTTCTCAGGAACTGGACCAGACGCCAAATACGCTTGGACAACAGGCGCTGGCATTTGCCGAGCCTTTCAACAACCAACACAACGCCGGAt                                              >  6:327324/1‑148 (MQ=255)
                                                                                                      ttGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGTCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGAt                                              <  1:191932/147‑1 (MQ=255)
                                                                                                       tGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGAttt                                            >  7:64159/1‑148 (MQ=255)
                                                                                                         agaTCCCGTAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCTACACCCAACACAAAGCCGGATTTg                                           <  5:46576/147‑1 (MQ=255)
                                                                                                         agaTCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGATGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGa                                          >  1:377362/1‑148 (MQ=255)
                                                                                                          gaTCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGAt                                         <  3:363478/148‑1 (MQ=255)
                                                                                                            tCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAacac                                                        <  7:366662/131‑1 (MQ=255)
                                                                                                            tCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAacac                                                        >  8:366662/1‑131 (MQ=255)
                                                                                                            tCCCGGAGAAATTAATGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAac                                                          <  7:206262/129‑1 (MQ=255)
                                                                                                            tCCCGGAGAAATTAATGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAac                                                          >  8:206262/1‑129 (MQ=255)
                                                                                                                 gagaAATTACTGAATACCGGGTCGCTTGGCGGCACTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAAAACCCAACAAAAAGCCGGATTTGATGCTAACg                                  <  3:72808/148‑1 (MQ=255)
                                                                                                                 gagaAATTACTGAATACCGGGTCGCGGGGCGGCACGCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCGGGCATTGGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACg                                  <  3:77757/148‑1 (MQ=255)
                                                                                                                  agaAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAgg                                                                         <  5:378331/108‑1 (MQ=255)
                                                                                                                  agaAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAgg                                                                         >  6:378331/1‑108 (MQ=255)
                                                                                                                    aaaTTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCg                               <  1:199994/148‑1 (MQ=255)
                                                                                                                     aaTTACTGAATACCGTGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACaa                                                                                               <  5:57867/83‑1 (MQ=255)
                                                                                                                     aaTTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACaa                                                                                               >  6:57867/1‑83 (MQ=255)
                                                                                                                              aTACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTg                      >  1:248478/1‑147 (MQ=255)
                                                                                                                                 ccGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAg                   <  2:44775/147‑1 (MQ=255)
                                                                                                                                  cGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCg                               <  3:27741/134‑1 (MQ=255)
                                                                                                                                  cGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCg                               >  4:27741/1‑134 (MQ=255)
                                                                                                                                    ggTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGAttt               <  2:248478/148‑1 (MQ=255)
                                                                                                                                      tCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCt             <  3:257782/148‑1 (MQ=255)
                                                                                                                                         cTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGTTTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTg          >  3:317817/1‑148 (MQ=255)
                                                                                                                                            ggcggcATTCTGACATACCGTTCTCAGGATCGGGACCAGACGCGGAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTGCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTa       <  3:352305/148‑1 (MQ=255)
                                                                                                                                             gcggcATTCTGACTTTCCGTTCTCAGGATCTGGACCAGACTCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACTCAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTAt      <  7:55943/148‑1 (MQ=255)
                                                                                                                                                  aTTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGAt                                              <  2:123730/103‑1 (MQ=255)
                                                                                                                                                  aTTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGAt                                              >  1:123730/1‑103 (MQ=255)
                                                                                                                                                   ttCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGt  >  3:162003/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACCGGGTCGCTGGGCGGCATTCTGACATTCCGTTCTCAGGATCTGGACCAGACGCGTAATACGCTTGGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGT  >  NC_000913/1139078‑1139370

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: