Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,106,370 G→C 5.5% intergenic (+50/‑99) cpxP → / → fieF inhibitor of the cpx response; periplasmic adaptor protein/ferrous iron and zinc transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,106,3700GC 5.5% 329.7 / 10.2 109intergenic (+50/‑99)cpxP/fieFinhibitor of the cpx response; periplasmic adaptor protein/ferrous iron and zinc transporter
Reads supporting (aligned to +/- strand):  ref base G (51/52);  new base C (3/3);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

Only 100 of 109 total aligned reads displayed.

CAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTG  >  NC_000913/4106228‑4106514
                                                                                                                                              |                                                                                                                                                
ctccaacGAAAGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTggg                                                                                                                                           <  2:1157272/147‑1 (MQ=255)
 aacaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAATCCCTGTTTTCCTTGCCATAGAAACCATCCCTGTCTTCCCCACCATGCTGTgggg                                                                                                                                          >  2:144379/1‑150 (MQ=255)
  acaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGCAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTTCCATAGACACGATCCCTGTCTTCCCCCACATGCTGgggggg                                                                                                                                         >  2:124817/1‑150 (MQ=255)
   caacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCGGTGGGGGt                                                                                                                                        >  2:1542206/1‑150 (MQ=255)
    aacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtt                                                                                                                                       <  2:842397/150‑1 (MQ=255)
     acGACTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTgggggaag                                                                                                                                      <  1:965234/150‑4 (MQ=255)
     acGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCACCTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttt                                                                                                                                      >  2:1644668/1‑150 (MQ=255)
     acGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTgggggaag                                                                                                                                      <  1:1447305/150‑4 (MQ=255)
      cGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGAAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGGCGGGGGGGGgttt                                                                                                                                     >  1:319859/1‑150 (MQ=255)
       gAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttttt                                                                                                                                    >  1:269866/1‑150 (MQ=255)
        aaTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttttt                                                                                                                                   <  2:43749/150‑1 (MQ=255)
        aaTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttttt                                                                                                                                   <  1:1599995/150‑1 (MQ=255)
         aTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttttttt                                                                                                                                  <  1:1021320/150‑1 (MQ=255)
         aTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttttttt                                                                                                                                  <  2:35684/150‑1 (MQ=255)
         aTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTgggggaagacag                                                                                                                                  <  1:43758/150‑8 (MQ=255)
           ggAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTa                                                                                                                                <  2:331276/150‑1 (MQ=255)
           ggAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGGTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTa                                                                                                                                <  1:1304691/150‑1 (MQ=255)
            gAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTAt                                                                                                                               <  1:1073476/150‑1 (MQ=255)
             aGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATc                                                                                                                              >  2:761037/1‑150 (MQ=255)
              gCATTTGCGTGACGTGACGCATTGGCAAAAAATTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTTATCCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATcc                                                                                                                             <  1:268581/150‑1 (MQ=255)
                aGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCCTCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGGTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACCTGCTGTGGGGGTTTTTTTTATCCTc                                                                                                                           <  1:1681528/150‑1 (MQ=255)
                      gtgacgtgacGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTg                                                                                                                     >  1:573944/1‑150 (MQ=255)
                             gacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctgctg                                                                                                              >  1:1528931/1‑150 (MQ=255)
                              acgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCtgctgc                                                                                                             <  1:1608805/150‑1 (MQ=255)
                              acgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCtgctgc                                                                                                             <  1:1601193/150‑1 (MQ=255)
                              acgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCtgctgc                                                                                                             <  1:1203904/150‑1 (MQ=255)
                              acgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCtgctgc                                                                                                             <  1:285710/150‑1 (MQ=255)
                              acgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCtgctgc                                                                                                             <  2:802031/150‑1 (MQ=255)
                                gCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCtt                                                                                                           <  2:1530527/150‑1 (MQ=255)
                                   aTGGCAAAAAATTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAAt                                                                                                        <  2:520791/150‑1 (MQ=255)
                                    aggCAAATAATATCAGCGTTGAAGCTATTGAGTACTCGCCACTCCCGCTCCCACTACTAAACCCTGTATTCCATCCCAACGACACCATCCCTGTCTTCCCCCCCCTGCTGTGGTGGTTTTTTTTATCCTCAATTGGCCTGCTGCTTAATg                                                                                                       <  2:324739/149‑1 (MQ=255)
                                       cAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                                                                                    >  2:793130/1‑150 (MQ=255)
                                       cAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                                                                                    >  2:809318/1‑150 (MQ=255)
                                       cAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGATGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                                                                                    >  2:666175/1‑150 (MQ=255)
                                         aaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                                                                                  <  1:51863/150‑1 (MQ=255)
                                         aaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                                                                                  >  1:652588/1‑150 (MQ=255)
                                         aaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACAATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAAAGCTTTg                                                                                                  >  1:357386/1‑150 (MQ=255)
                                          aaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCTACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGc                                                                                                 <  1:958514/150‑1 (MQ=255)
                                            aaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAg                                                                                               <  1:902271/150‑1 (MQ=255)
                                              tttCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAgat                                                                                             <  1:1282050/149‑1 (MQ=255)
                                                    ggTTGAAGCAATTGAGAAGTAGCTACTCACGTTGCCAGTAGTAAACCGTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCGTATTGCATTGCAGATGATTTg                                                                                       <  2:568258/149‑1 (MQ=255)
                                                      ttGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCt                                                                                     <  1:77764/150‑1 (MQ=255)
                                                       tGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGGTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCtt                                                                                    >  2:1656844/1‑150 (MQ=255)
                                                        gAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTACCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGAAGATTATTTGCCTc                                                                                   >  2:1435778/1‑150 (MQ=255)
                                                           gCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGt                                                                                >  1:308206/1‑150 (MQ=255)
                                                           gCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGt                                                                                <  1:1092995/150‑1 (MQ=255)
                                                             tATTGAGTAGTAGCAACTCTCGTTCCCAGTAGTAACCCCTGATTTGCTTGCCATAGACACCATCCCTGTCTTCCACCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGGAGATGATTTGCTTCCGttt                                                                              >  2:936673/1‑149 (MQ=255)
                                                             tATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTta                                                                              <  2:964638/150‑1 (MQ=255)
                                                             tATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTta                                                                              <  1:1709330/150‑1 (MQ=255)
                                                             tATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTta                                                                              <  1:1583889/150‑1 (MQ=255)
                                                             tATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGGTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGttt                                                                              >  2:592953/1‑149 (MQ=255)
                                                              aTTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACCTGCGGTGGGGGGTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATTATTTGCTTCCGTtat                                                                             >  1:1386069/1‑150 (MQ=255)
                                                              aTTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTtat                                                                             >  2:1072391/1‑150 (MQ=255)
                                                              aTTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTtat                                                                             <  1:733649/150‑1 (MQ=255)
                                                               ttGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGTTTTGCTTCCGTtatt                                                                            >  2:1330346/1‑149 (MQ=255)
                                                               ttGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTtata                                                                            >  1:1141106/1‑150 (MQ=255)
                                                                 gAGTAGTAGCAACTCCCGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGGGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACt                                                                          >  1:826394/1‑150 (MQ=255)
                                                                   gtagtaGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCATGTCTTCCCCCACATGCTGTGGGGGTTTATTTTATCCTCAATTTGCCTGAGGGTTAAAGCATTGCAGATGATTTGCTTCAGTTATACTAg                                                                        >  2:1587869/1‑150 (MQ=255)
                                                                     agtagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCg                                                                      <  1:1487667/150‑1 (MQ=255)
                                                                     agtagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGGTTTTTTTATGCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCg                                                                      >  1:1754328/1‑150 (MQ=255)
                                                                     agtagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGGCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCg                                                                      <  1:1480166/150‑1 (MQ=255)
                                                                      gtagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTTATGCATTGCAGATGATTTGCTTCCGTTATACTAGCgg                                                                     >  1:1730777/1‑149 (MQ=255)
                                                                        agCAACTCACGTTCCCAGTAGTAAACCCTGGTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCGGTGGGGGTTTTTTTTATCCTCAATTTGCGTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCa                                                                   >  2:226635/1‑150 (MQ=255)
                                                                           aaCTCACGTTCCCAGTATTAAACCCTGTTTTCCGTGCCATAGCCACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGtt                                                                <  2:1401206/150‑1 (MQ=255)
                                                                           aaCTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCCTCCCTGTCTTCCCCCACATGCTTTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCTGATGATTTGCTTCCGTTATACTAGCGTCAGtt                                                                <  1:1841978/150‑1 (MQ=255)
                                                                             cTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGa                                                              <  1:1707506/150‑1 (MQ=255)
                                                                               cACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATACCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATa                                                            <  2:653819/150‑1 (MQ=255)
                                                                                   ttCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCgg                                                        <  2:1261365/150‑1 (MQ=255)
                                                                                    tCCCAGTAGTAACCCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                       <  2:269866/150‑1 (MQ=255)
                                                                                    tCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                       <  2:1655092/150‑1 (MQ=255)
                                                                                     cccAGTAGTAAACCCTGTTTTACTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGa                                                      <  1:144379/150‑1 (MQ=255)
                                                                                      ccAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCCCATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGGTGCTTAATGCATTGCAGATGATTTGCTTCCGGTTTAATAGGGTCAGTTGATAGCGGGAg                                                     >  2:1183885/1‑150 (MQ=255)
                                                                                             taAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTAt                                              >  2:970872/1‑150 (MQ=255)
                                                                                                aCCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCACTTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaa                                           <  2:1636840/150‑1 (MQ=255)
                                                                                                 cccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTTATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaat                                          >  2:1462872/1‑150 (MQ=255)
                                                                                                 cccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaat                                          >  2:1507960/1‑150 (MQ=255)
                                                                                                 cccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTACCCCACCAGCTGTGGGGGTTTTTTTTATCCTCAAATTGCATGAATCTAAATGCATTGCAGATGATTTGCGTACGATATACTAGCGTAAATTGATAGCCGGAGTATTTATGaat                                          >  1:1372876/1‑150 (MQ=255)
                                                                                                   cTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatca                                        >  2:9436/1‑150 (MQ=255)
                                                                                                    tGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTTTGaatcaa                                       >  1:531886/1‑150 (MQ=255)
                                                                                                    tGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatcaa                                       <  1:433054/150‑1 (MQ=255)
                                                                                                       tttCCTTGCCATCCACACCTTCCCTGTGTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCt                                    <  1:1301557/150‑1 (MQ=255)
                                                                                                        ttCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCtt                                   <  2:43302/150‑1 (MQ=255)
                                                                                                         tCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTa                                  >  2:469433/1‑150 (MQ=255)
                                                                                                          ccTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTAt                                 >  2:686158/1‑150 (MQ=255)
                                                                                                              gCCATAGCCTCCATCACTGTCTTCCCCCACATGCTGTGGGGGTTTTGTTTATCCTCAATTTGCCTGCTGCTTAATGAATTGCAGATGATTTGCTTCCGTTATACTAGCGGCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGAc                             >  2:74883/1‑150 (MQ=255)
                                                                                                              gCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGAc                             <  2:1375089/150‑1 (MQ=255)
                                                                                                                    gccacCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTgg                       >  2:1723585/3‑150 (MQ=255)
                                                                                                                      cacGATCCCTGTCTTCCCGCACCTGCTGTGGGGGGTTTTTTTATCCTCAATTTTCCTGCTGCTTAATTCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGGTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTc                     <  1:774448/150‑1 (MQ=255)
                                                                                                                      cacCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTc                     <  1:1458391/150‑1 (MQ=255)
                                                                                                                                 tCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGggcggc          >  2:50963/1‑150 (MQ=255)
                                                                                                                                  cTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCTGgctgcg         >  1:36995/1‑150 (MQ=255)
                                                                                                                                   ttCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGGTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGa        <  1:1286295/150‑1 (MQ=255)
                                                                                                                                    tCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGAt       >  1:515694/1‑150 (MQ=255)
                                                                                                                                    acccccACAGCATGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGAt       >  1:852004/2‑150 (MQ=255)
                                                                                                                                    acccccACAGCATGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGAt       >  1:1871589/2‑150 (MQ=255)
                                                                                                                                      ccccACAGCATGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGATACAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtg     >  2:596568/1‑150 (MQ=255)
                                                                                                                                       cccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtgc    >  1:1095894/1‑150 (MQ=255)
                                                                                                                                       cccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtgc    >  1:525007/1‑150 (MQ=255)
                                                                                                                                        ccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtgct   >  1:850206/1‑150 (MQ=255)
                                                                                                                                         cacaTGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtgctg  >  2:1712888/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTG  >  NC_000913/4106228‑4106514

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