Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 779,659 C→A 5.1% P21H (CCC→CAC) ‡ ybgF → periplasmic TolA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913779,6590CA 5.1% 546.0 / 12.2 156P21H (CCC→CAC) ‡ybgFperiplasmic TolA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (85/63);  new base A (2/6);  total (87/69)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

GGTAAAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATT  >  NC_000913/779511‑779807
                                                                                                                                                    |                                                                                                                                                    
 gTAAAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGccccc                                                                                                                                                     >  2:130241/1‑149 (MQ=255)
    aaGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTgg                                                                                                                                                  <  1:235437/149‑1 (MQ=255)
    aaGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTgg                                                                                                                                                  >  2:195166/1‑149 (MQ=255)
     aGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTggg                                                                                                                                                 <  2:580337/149‑1 (MQ=255)
      gAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGc                                                                                                                                                >  2:13193/1‑149 (MQ=255)
          aaCCTGCAGTACGGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCt                                                                                                                                            <  1:951171/149‑1 (MQ=255)
           aCCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCtt                                                                                                                                           <  4:502424/149‑1 (MQ=255)
           aCCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCta                                                                                                                                           <  2:456776/149‑2 (MQ=255)
            ccTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCtat                                                                                                                                          <  2:1070386/149‑3 (MQ=255)
            ccTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCtat                                                                                                                                          <  4:370561/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)
                     cTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGc                                                                                                                                 <  4:604574/149‑1 (MQ=255)
                            aTGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATc                                                                                                                          >  3:189417/1‑149 (MQ=255)
                            aTGACGAAGCGACATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATc                                                                                                                          <  3:788200/149‑1 (MQ=255)
                              gACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAAt                                                                                                                           >  3:557345/1‑146 (MQ=255)
                              gACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAAt                                                                                                                           <  4:557345/146‑1 (MQ=255)
                              gACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCag                                                                                                                        >  2:996993/1‑149 (MQ=255)
                                 gAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAAt                                                                                                                           >  4:79703/1‑143 (MQ=255)
                                        cATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTGAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGc                                                                                                              <  3:644055/149‑1 (MQ=255)
                                          tACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAAt                                                                                                                           >  3:11776/1‑134 (MQ=255)
                                          tACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAAt                                                                                                                           <  4:11776/134‑1 (MQ=255)
                                            cTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAg                                                                                                          >  2:904905/1‑149 (MQ=255)
                                             tCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAgg                                                                                                         >  2:666995/1‑149 (MQ=255)
                                             tCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAgg                                                                                                         >  2:502044/1‑149 (MQ=255)
                                               cAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCt                                                                                                       >  1:1071470/1‑149 (MQ=255)
                                                aaaaaCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCtata                                                                                                                                         >  3:469408/1‑111 (MQ=255)
                                                    aCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCagta                                                                                                                      <  3:814456/129‑1 (MQ=255)
                                                      cgtcgtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGa                                                                                                >  4:1029176/1‑149 (MQ=255)
                                                       gtcgtGCGGTACTGGTTTACTAAGAGAATTTCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGaa                                                                                               <  1:730945/149‑1 (MQ=255)
                                                       gtcgtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGaa                                                                                               >  3:981103/1‑149 (MQ=255)
                                                       gtcgtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGaa                                                                                               <  3:326089/149‑1 (MQ=255)
                                                       gtcgtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGaa                                                                                               <  2:395044/149‑1 (MQ=255)
                                                            gCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACcg                                                                                          >  4:1061085/1‑149 (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)
                                                                     gTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATc                                                                                                                          <  3:125798/108‑1 (MQ=255)
                                                                     gTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCa                                                                                 >  2:85864/1‑149 (MQ=255)
                                                                          cTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTg                                                                            >  2:518395/1‑149 (MQ=255)
                                                                           tAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTTTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGa                                                                           <  3:679873/149‑1 (MQ=255)
                                                                           tAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGa                                                                           >  3:354707/1‑149 (MQ=255)
                                                                           tAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGa                                                                           <  3:1069570/149‑1 (MQ=255)
                                                                             agagAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                         >  4:677581/1‑149 (MQ=255)
                                                                             agagAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                         >  3:445463/1‑149 (MQ=255)
                                                                             agagAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                         >  1:833438/1‑149 (MQ=255)
                                                                              gagaATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCg                                                                        >  1:396504/1‑149 (MQ=255)
                                                                                  aTTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAgg                                                                                                         <  3:56200/112‑1 (MQ=255)
                                                                                   ttGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTAttt                                                                   >  2:922536/1‑149 (MQ=255)
                                                                                       aTGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTaa                                                               >  4:672901/1‑149 (MQ=255)
                                                                                       aTGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTaa                                                               <  1:518395/149‑1 (MQ=255)
                                                                                       aTGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTaa                                                               <  3:311532/149‑1 (MQ=255)
                                                                                              gTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTcac                                                        >  4:432440/1‑149 (MQ=255)
                                                                                               tAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTcaca                                                       >  4:972527/1‑149 (MQ=255)
                                                                                                 aCTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTAt                                                                     <  2:559746/133‑1 (MQ=255)
                                                                                                 aCTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTAt                                                                     >  1:559746/1‑133 (MQ=255)
                                                                                                  cTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCt                                                           >  3:144999/1‑142 (MQ=255)
                                                                                                  cTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCAcagcc                                                    <  4:445463/149‑1 (MQ=255)
                                                                                                   ttCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCAcagcca                                                   >  3:390813/1‑149 (MQ=255)
                                                                                                   ttCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCAcagcca                                                   >  1:736671/1‑149 (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:1071613/149‑1 (MQ=255)
                                                                                                           tCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                           <  4:22489/149‑1 (MQ=255)
                                                                                                            cAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAc                                          <  3:852555/149‑1 (MQ=255)
                                                                                                             aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCtttt                                             >  4:357898/1‑145 (MQ=255)
                                                                                                             aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                           >  1:189499/1‑147 (MQ=255)
                                                                                                             aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                           <  1:623448/147‑1 (MQ=255)
                                                                                                             aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                           <  2:189499/147‑1 (MQ=255)
                                                                                                             aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAcc                                         >  2:82943/1‑149 (MQ=255)
                                                                                                             aaCTATTGACTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTATTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAcc                                         >  3:964393/1‑149 (MQ=255)
                                                                                                              aCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAccc                                        >  4:964938/1‑149 (MQ=255)
                                                                                                               cTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCa                                       >  2:523684/1‑149 (MQ=255)
                                                                                                                tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGGTTTTAACCCaa                                      <  1:577069/149‑1 (MQ=255)
                                                                                                                tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                      >  1:155715/1‑149 (MQ=255)
                                                                                                                tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                      >  2:807879/1‑149 (MQ=255)
                                                                                                                tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                      >  1:695031/1‑149 (MQ=255)
                                                                                                                 aTTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAAc                                     >  1:862509/1‑149 (MQ=255)
                                                                                                                  ttGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACt                                    >  2:687801/1‑149 (MQ=255)
                                                                                                                    gAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTcc                                  <  2:979828/149‑1 (MQ=255)
                                                                                                                     aGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCa                                 <  3:990335/149‑1 (MQ=255)
                                                                                                                         tGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGcaa                             >  2:312851/1‑149 (MQ=255)
                                                                                                                            cGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGcaacaa                          <  3:485273/149‑1 (MQ=255)
                                                                                                                             gTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCaacaac                         >  1:44990/1‑149 (MQ=255)
                                                                                                                             gTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCaacaac                         >  4:466661/1‑149 (MQ=255)
                                                                                                                             gTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCaacaac                         >  2:162118/1‑149 (MQ=255)
                                                                                                                               tACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACtt                       >  3:235944/1‑149 (MQ=255)
                                                                                                                                 cTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTc                     >  4:213822/1‑149 (MQ=255)
                                                                                                                                    gttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGa                  <  4:189417/149‑1 (MQ=255)
                                                                                                                                    gttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGa                  >  1:1035454/1‑149 (MQ=255)
                                                                                                                                    gttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCAAACTCCAGCAACAACTTTCTGa                  <  1:130241/149‑1 (MQ=255)
                                                                                                                                     ttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGAt                 <  1:677579/149‑1 (MQ=255)
                                                                                                                                           aaGCGGCCCAGGGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaatcaa           >  2:1081514/2‑149 (MQ=255)
                                                                                                                                             gCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaatcaatc         >  4:800084/1‑149 (MQ=255)
                                                                                                                                                   cccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGatat   >  3:84846/1‑149 (MQ=255)
                                                                                                                                                    ccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATAtt  >  2:1024631/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTAAAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATT  >  NC_000913/779511‑779807

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: