Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 779,661 T→G 5.0% W22G (TGG→GGG)  ybgF → periplasmic TolA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913779,6610TG 5.0% 550.4 / 18.1 160W22G (TGG→GGG) ybgFperiplasmic TolA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (87/64);  new base G (2/6);  total (89/71)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 161 total aligned reads displayed.

AAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGA  >  NC_000913/779515‑779809
                                                                                                                                                  |                                                                                                                                                    
aaGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTgg                                                                                                                                                    >  2:195166/1‑149 (MQ=255)
aaGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTgg                                                                                                                                                    <  1:235437/149‑1 (MQ=255)
aaGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTgg                                                                                                                                                    >  1:589711/1‑149 (MQ=255)
 aGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTggg                                                                                                                                                   <  2:580337/149‑1 (MQ=255)
 aGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTggg                                                                                                                                                   >  2:933787/1‑149 (MQ=255)
  gAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGc                                                                                                                                                  >  2:13193/1‑149 (MQ=255)
  gAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGc                                                                                                                                                  <  1:696248/149‑1 (MQ=255)
       aCCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCta                                                                                                                                             <  1:691038/149‑2 (MQ=255)
       aCCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCta                                                                                                                                             <  2:456776/149‑2 (MQ=255)
        ccTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCtat                                                                                                                                            <  4:370561/149‑3 (MQ=255)
        ccTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCtat                                                                                                                                            <  2:1070386/149‑3 (MQ=255)
            cAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCt                                                                                                                                        <  1:373097/149‑1 (MQ=255)
              gTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAAGTGCATGAGCAGTAACTTCAGACAGCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCa                                                                                                                                      <  4:743152/149‑1 (MQ=255)
                 cTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGc                                                                                                                                   >  4:824628/1‑149 (MQ=255)
                      tCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCaa                                                                                                                              >  3:676461/1‑149 (MQ=255)
                      tCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACGAAGAGAATTGCATGAGTAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGACCCCCTGTGCCGCTTTTGCTCAGGCACCaa                                                                                                                              >  2:835350/1‑149 (MQ=255)
                        aTGGCGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGCTTACTAAGAGAATTGCATGAGCAATAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATc                                                                                                                            >  2:677579/1‑149 (MQ=255)
                        aTGACGAAGCGACATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATc                                                                                                                            <  3:788200/149‑1 (MQ=255)
                          gACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAAt                                                                                                                             <  4:557345/146‑1 (MQ=255)
                          gACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCag                                                                                                                          >  2:823766/1‑149 (MQ=255)
                             gAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAAt                                                                                                                             <  3:79703/143‑1 (MQ=255)
                             gAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAAt                                                                                                                             >  4:79703/1‑143 (MQ=255)
                                   gCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCgg                                                                                                                 <  2:188847/149‑1 (MQ=255)
                                    cATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTGAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGc                                                                                                                <  3:644055/149‑1 (MQ=255)
                                      tACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAAt                                                                                                                             >  3:11776/1‑134 (MQ=255)
                                         tCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAgg                                                                                                           >  4:640770/1‑149 (MQ=255)
                                         tCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAgg                                                                                                           >  2:553219/1‑149 (MQ=255)
                                           cAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCt                                                                                                         >  1:1071470/1‑149 (MQ=255)
                                            aaaaaCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTc                                                                                                        >  3:438070/1‑149 (MQ=255)
                                            aaaaaCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCtata                                                                                                                                           <  4:469408/114‑4 (MQ=255)
                                                  cgtcgtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGa                                                                                                  >  4:1029176/1‑149 (MQ=255)
                                                   gtcgtGCGGTACTGGTTTACTAAGAGAATTTCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGaa                                                                                                 <  1:730945/149‑1 (MQ=255)
                                                   gtcgtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGaa                                                                                                 <  3:326089/149‑1 (MQ=255)
                                                         cGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACcgc                                                                                           <  3:1014853/149‑1 (MQ=255)
                                                           gTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCTGCTCAGGCTCGGTCGAAGACCGCGt                                                                                         >  1:79117/1‑149 (MQ=255)
                                                           gTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGt                                                                                         <  4:916997/149‑1 (MQ=255)
                                                            tACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTc                                                                                        >  1:312472/1‑149 (MQ=255)
                                                             aCTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCa                                                                                       <  1:891371/149‑1 (MQ=255)
                                                                 gTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATc                                                                                                                            >  4:125798/1‑108 (MQ=255)
                                                                 gTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCa                                                                                   >  2:85864/1‑149 (MQ=255)
                                                                      cTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTg                                                                              >  2:518395/1‑149 (MQ=255)
                                                                       tAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGa                                                                             >  3:354707/1‑149 (MQ=255)
                                                                       tAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGa                                                                             <  3:1069570/149‑1 (MQ=255)
                                                                         agagAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                           <  3:264091/149‑1 (MQ=255)
                                                                         agagAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                           >  1:833438/1‑149 (MQ=255)
                                                                         agagAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                           >  4:677581/1‑149 (MQ=255)
                                                                          gagaATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCg                                                                          >  1:396504/1‑149 (MQ=255)
                                                                          gagaATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCg                                                                          >  1:479844/1‑149 (MQ=255)
                                                                              aTTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAgg                                                                                                           <  3:56200/112‑1 (MQ=255)
                                                                              aTTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAgg                                                                                                           >  4:56200/1‑112 (MQ=255)
                                                                              aTTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTAtt                                                                      >  1:750220/1‑149 (MQ=255)
                                                                               ttGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTAttt                                                                     >  1:547179/1‑149 (MQ=255)
                                                                                 gCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCt                                                                   <  2:479844/149‑1 (MQ=255)
                                                                                           tAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTcaca                                                         >  4:972527/1‑149 (MQ=255)
                                                                                             aCTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTAt                                                                       >  1:559746/1‑133 (MQ=255)
                                                                                              cTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCt                                                             <  4:144999/142‑1 (MQ=255)
                                                                                              cTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCt                                                             >  3:144999/1‑142 (MQ=255)
                                                                                               ttCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCAcagcca                                                     >  2:221913/1‑149 (MQ=255)
                                                                                                tCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTaa                                                                 >  3:136738/1‑136 (MQ=255)
                                                                                                tCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTaa                                                                 <  4:136738/136‑1 (MQ=255)
                                                                                                 cAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCAcagccagc                                                   <  1:922536/149‑1 (MQ=255)
                                                                                                      aTCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTa                                              >  3:417999/1‑149 (MQ=255)
                                                                                                       tCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                             <  4:22489/149‑1 (MQ=255)
                                                                                                       tCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                             <  1:865509/149‑1 (MQ=255)
                                                                                                       tCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                             <  4:1071613/149‑1 (MQ=255)
                                                                                                        cAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAc                                            <  3:852555/149‑1 (MQ=255)
                                                                                                         aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCtttt                                               <  3:357898/145‑1 (MQ=255)
                                                                                                         aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                             >  1:189499/1‑147 (MQ=255)
                                                                                                         aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                             <  1:623448/147‑1 (MQ=255)
                                                                                                         aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                             >  2:623448/1‑147 (MQ=255)
                                                                                                         aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAcc                                           >  4:103444/1‑149 (MQ=255)
                                                                                                         aaCTATTGACTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTATTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAcc                                           >  3:964393/1‑149 (MQ=255)
                                                                                                          aCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAccc                                          >  4:964938/1‑149 (MQ=255)
                                                                                                            tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                        >  1:155715/1‑149 (MQ=255)
                                                                                                            tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                        >  2:807879/1‑149 (MQ=255)
                                                                                                            tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                        <  3:496562/149‑1 (MQ=255)
                                                                                                            tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                        >  2:514750/1‑149 (MQ=255)
                                                                                                             aTTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAAc                                       >  3:454080/1‑149 (MQ=255)
                                                                                                                gAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTcc                                    <  2:979828/149‑1 (MQ=255)
                                                                                                                 aGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCa                                   <  3:990335/149‑1 (MQ=255)
                                                                                                                        cGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGcaacaa                            <  3:485273/149‑1 (MQ=255)
                                                                                                                         gTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCaacaac                           >  2:162118/1‑149 (MQ=255)
                                                                                                                         gTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCaacaac                           >  4:466661/1‑149 (MQ=255)
                                                                                                                           tACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACtt                         >  3:235944/1‑149 (MQ=255)
                                                                                                                             cTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTc                       >  4:213822/1‑149 (MQ=255)
                                                                                                                               ggttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTg                     <  3:972527/149‑1 (MQ=255)
                                                                                                                                gttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGa                    >  1:1035454/1‑149 (MQ=255)
                                                                                                                                gttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCAAACTCCAGCAACAACTTTCTGa                    <  1:130241/149‑1 (MQ=255)
                                                                                                                                   ggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaa                 >  4:1077900/1‑149 (MQ=255)
                                                                                                                                      ataGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaatca              >  4:414258/1‑149 (MQ=255)
                                                                                                                                       aaGCGGCCCAGGGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaatcaa             >  2:1081514/2‑149 (MQ=255)
                                                                                                                                         gCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaatcaatc           <  4:438070/149‑1 (MQ=255)
                                                                                                                                         gCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaatcaatc           <  4:249945/149‑1 (MQ=255)
                                                                                                                                         gCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaatcaatc           >  4:800084/1‑149 (MQ=255)
                                                                                                                                            gCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGa        <  2:918310/149‑1 (MQ=255)
                                                                                                                                              ccccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGAGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGata      <  1:996993/149‑1 (MQ=255)
                                                                                                                                               cccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGatat     <  2:319559/149‑1 (MQ=255)
                                                                                                                                               cccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGatat     >  3:84846/1‑149 (MQ=255)
                                                                                                                                                ccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATAtt    >  2:1024631/1‑149 (MQ=255)
                                                                                                                                                  tGGGCCGCTTTTGCTCAGGCACCAATCAGGAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGa  <  3:815376/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGA  >  NC_000913/779515‑779809

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