Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,056,039 C→A 5.4% D101Y (GAC→TAC) ‡ glnL ← sensory histidine kinase in two‑component regulatory system with GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,056,0390CA 5.4% 427.8 / 7.7 124D101Y (GAC→TAC) ‡glnLsensory histidine kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base C (64/53);  new base A (4/3);  total (68/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 128 total aligned reads displayed.

CCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAAC  >  NC_000913/4055891‑4056187
                                                                                                                                                    |                                                                                                                                                    
    gTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGc                                                                                                                                                  <  4:1040421/149‑1 (MQ=255)
       aaCCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATa                                                                                                                                               >  1:371929/1‑149 (MQ=255)
        aCCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGGACGGCAtga                                                                                                                                              <  1:163577/149‑3 (MQ=255)
          cGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGGACGGCAtgatc                                                                                                                                            <  4:1004545/149‑5 (MQ=255)
             cAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTgg                                                                                                                                         <  4:726780/149‑1 (MQ=255)
                   gATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTc                                                                                                                                   <  3:462312/149‑1 (MQ=255)
                        ttAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGa                                                                                                                              >  3:404622/1‑149 (MQ=255)
                        ttAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGa                                                                                                                              >  3:692915/1‑149 (MQ=255)
                          aaTCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGaaa                                                                                                                            <  4:687159/149‑1 (MQ=255)
                           aTCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAg                                                                                                                           <  3:546554/149‑1 (MQ=255)
                            tctcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAgg                                                                                                                          >  4:552498/1‑149 (MQ=255)
                            tctcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAgg                                                                                                                          >  4:378494/1‑149 (MQ=255)
                             ctcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGa                                                                                                                         <  2:231924/149‑1 (MQ=255)
                              tcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGat                                                                                                                        >  3:970365/1‑149 (MQ=255)
                              tcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGat                                                                                                                        >  1:852272/1‑149 (MQ=255)
                              tcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGat                                                                                                                        >  2:811743/1‑149 (MQ=255)
                              tcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGat                                                                                                                        >  1:972800/1‑149 (MQ=255)
                                  gtgCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGc                                                                                                                    >  1:895341/1‑149 (MQ=255)
                                  gtgCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGc                                                                                                                    <  3:632748/149‑1 (MQ=255)
                                    gCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCCGCATACGCTGGGCCGTCACAGAAAGGATATGCGa                                                                                                                  <  3:1039880/149‑1 (MQ=255)
                                      cAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAgc                                                                                                                <  2:326816/149‑1 (MQ=255)
                                        ggCCGCGCACTGAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAgcgc                                                                                                              >  1:818646/1‑149 (MQ=255)
                                              gcACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCg                                                                                                        <  4:563691/149‑1 (MQ=255)
                                               cACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGTGCCCGTCGa                                                                                                       <  2:741027/149‑1 (MQ=255)
                                               cACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGa                                                                                                       >  1:398647/1‑149 (MQ=255)
                                               cACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGa                                                                                                       <  3:748115/149‑1 (MQ=255)
                                                aCTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAt                                                                                                      <  1:811743/149‑1 (MQ=255)
                                                aCTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAt                                                                                                      <  4:347812/149‑1 (MQ=255)
                                                  tAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGa                                                                                                    >  4:989855/1‑149 (MQ=255)
                                                  tAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGa                                                                                                    >  1:79917/1‑149 (MQ=255)
                                                  tAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGa                                                                                                    >  1:285397/1‑149 (MQ=255)
                                                    aaTCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGAcc                                                                                                  >  2:432783/1‑149 (MQ=255)
                                                      tCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAg                                                                                                >  4:290017/1‑149 (MQ=255)
                                                       cACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGc                                                                                               <  1:444242/149‑1 (MQ=255)
                                                                caACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCg                                                                                      >  1:126085/1‑149 (MQ=255)
                                                                   ccTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTa                                                                                   >  2:1051240/1‑149 (MQ=255)
                                                                     tgctgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATc                                                                                 <  2:483656/149‑1 (MQ=255)
                                                                     tgctgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATc                                                                                 <  3:364541/149‑1 (MQ=255)
                                                                      gctgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCg                                                                                <  2:126085/149‑1 (MQ=255)
                                                                      gctgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCg                                                                                <  2:371929/149‑1 (MQ=255)
                                                                        tgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGt                                                                              <  2:320865/149‑1 (MQ=255)
                                                                         gggCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTa                                                                             >  1:264728/1‑149 (MQ=255)
                                                                          ggCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTaa                                                                            >  1:301247/1‑149 (MQ=255)
                                                                             gTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTc                                                                                       >  3:587950/1‑135 (MQ=255)
                                                                              tGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAAcc                                                                        >  1:715721/1‑149 (MQ=255)
                                                                                cTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCtt                                                                      <  1:500452/149‑1 (MQ=255)
                                                                                     gCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGcccc                                                                 >  2:990839/1‑149 (MQ=255)
                                                                                       tGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGc                                                               >  2:125524/1‑149 (MQ=255)
                                                                                        gTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGcc                                                              <  2:715721/149‑1 (MQ=255)
                                                                                           ccTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTcc                                                           <  3:901676/149‑1 (MQ=255)
                                                                                           ccTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTcc                                                           >  1:1058936/1‑149 (MQ=255)
                                                                                            cTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCa                                                          >  2:664966/1‑149 (MQ=255)
                                                                                              gACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCtt                                                                      <  2:554203/135‑1 (MQ=255)
                                                                                              gACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCtt                                                                      >  1:554203/1‑135 (MQ=255)
                                                                                                 ttAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACtt                                                     >  1:146721/1‑149 (MQ=255)
                                                                                                 ttAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACtt                                                     >  2:980708/1‑149 (MQ=255)
                                                                                                    aGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCt                                                  <  1:126952/149‑1 (MQ=255)
                                                                                                    aGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCt                                                  <  3:989855/149‑1 (MQ=255)
                                                                                                    aGCGGCGCTGGTTATCCATCGGAGCCACCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCAAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCt                                                  <  1:281452/149‑1 (MQ=255)
                                                                                                     gcggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGcc                                                                   >  1:393287/1‑131 (MQ=255)
                                                                                                     gcggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCtt                                                 >  1:1024239/1‑149 (MQ=255)
                                                                                                     gcggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCtt                                                 >  3:978513/1‑149 (MQ=255)
                                                                                                      cggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGGATACGCTGGGCCGTCACAGAAAAGATATGCGAGCCCCCGTCGATGACCAACGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTg                                                <  2:332382/149‑1 (MQ=255)
                                                                                                       ggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGc                                               >  2:812085/1‑149 (MQ=255)
                                                                                                          gcTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATc                                            <  4:978513/149‑1 (MQ=255)
                                                                                                           cTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCa                                           >  3:927213/1‑149 (MQ=255)
                                                                                                            tGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAg                                          <  3:378494/149‑1 (MQ=255)
                                                                                                             ggTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGc                                         >  1:829962/1‑149 (MQ=255)
                                                                                                              gTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCt                                        >  1:642819/1‑149 (MQ=255)
                                                                                                               ttATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTc                                       <  4:589639/149‑1 (MQ=255)
                                                                                                                tATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCg                                      <  2:818646/149‑1 (MQ=255)
                                                                                                                      tCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAttt                                >  4:562218/1‑149 (MQ=255)
                                                                                                                      tCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAttt                                >  3:353619/1‑149 (MQ=255)
                                                                                                                        ggAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAGAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTaa                              >  4:300765/1‑149 (MQ=255)
                                                                                                                        ggAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTaa                              <  1:980708/149‑1 (MQ=255)
                                                                                                                         gAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAAt                             <  2:278009/149‑1 (MQ=255)
                                                                                                                            ccATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGatat                                                                                                                      >  1:238373/1‑57 (MQ=255)
                                                                                                                            ccATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATgag                          <  4:128080/149‑1 (MQ=255)
                                                                                                                              aTCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGaa                        >  4:63309/1‑149 (MQ=255)
                                                                                                                               tctcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGCGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAg                       <  2:464437/149‑1 (MQ=255)
                                                                                                                               tctcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGcc                                                              <  4:289155/110‑1 (MQ=255)
                                                                                                                               tctcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGcc                                                              >  3:289155/1‑110 (MQ=255)
                                                                                                                               tctcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAg                       >  4:990270/1‑149 (MQ=255)
                                                                                                                               tctcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAg                       >  3:621880/1‑149 (MQ=255)
                                                                                                                               tctcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAg                       <  4:353619/149‑1 (MQ=255)
                                                                                                                                 tcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTa                     <  2:972800/149‑1 (MQ=255)
                                                                                                                                 tcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTa                     >  1:830792/1‑149 (MQ=255)
                                                                                                                                 tcCAGCAGGATCAGGCCGTCGGGCATACGCCGGGCCGTCACAGAAAGGATATGGGAGCGCCGTTCTATGACCAGCGTCACTCCTTTACCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTCGCATCAGCTCGATATTTAATGAGAAGTa                     <  2:832592/149‑1 (MQ=255)
                                                                                                                                    agcagGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCt                  >  3:981430/1‑149 (MQ=255)
                                                                                                                                     gctggATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTc                 >  1:940197/4‑149 (MQ=255)
                                                                                                                                     gcagGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTc                 >  1:675435/1‑149 (MQ=255)
                                                                                                                                      ccagcgtaTGCCGGACGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCa                >  4:327749/8‑149 (MQ=255)
                                                                                                                                       agGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAGAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCaa               <  3:300765/149‑1 (MQ=255)
                                                                                                                                         gATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCt                                                  >  2:418265/1‑112 (MQ=255)
                                                                                                                                         gATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa             >  4:1033211/1‑149 (MQ=255)
                                                                                                                                         gATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa             <  3:990270/149‑1 (MQ=255)
                                                                                                                                         gATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGGCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa             >  1:722469/1‑149 (MQ=255)
                                                                                                                                          aTCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAAc                                                                         >  1:37023/1‑88 (MQ=255)
                                                                                                                                          aTCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAAc                                                                         <  2:37023/88‑1 (MQ=255)
                                                                                                                                            cATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGtt          <  2:264728/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAAC  >  NC_000913/4055891‑4056187

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