Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,106,369 T→G 5.2% intergenic (+49/‑100) 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,3690TG 5.2% 298.8 / 6.3 100intergenic (+49/‑100)cpxP/fieFinhibitor of the cpx response; periplasmic adaptor protein/ferrous iron and zinc transporter
Reads supporting (aligned to +/- strand):  ref base T (39/55);  new base G (4/2);  total (43/57)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGA  >  NC_000913/4106222‑4106517
                                                                                                                                                   |                                                                                                                                                    
aaaCATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGc                                                                                                                                                    <  2:247880/150‑1 (MQ=255)
 aaCATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCt                                                                                                                                                   >  2:355432/1‑150 (MQ=255)
   cATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCtgt                                                                                                                                                 >  2:1348683/1‑150 (MQ=255)
    aTCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCtgtg                                                                                                                                                <  2:553860/150‑1 (MQ=255)
     tCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgg                                                                                                                                               <  2:1256885/150‑1 (MQ=255)
     tCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgg                                                                                                                                               <  2:1256924/150‑1 (MQ=255)
     tCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgg                                                                                                                                               <  2:1256879/150‑1 (MQ=255)
     tCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgg                                                                                                                                               >  2:641815/1‑150 (MQ=255)
       aacaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAACAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTACACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgggg                                                                                                                                             <  1:1004295/150‑1 (MQ=255)
       aacaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgggg                                                                                                                                             <  1:1077275/150‑1 (MQ=255)
       aacaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgggg                                                                                                                                             <  1:1077276/150‑1 (MQ=255)
       aacaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgggg                                                                                                                                             <  2:825704/150‑1 (MQ=255)
       aacaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgggg                                                                                                                                             <  1:279356/150‑1 (MQ=255)
        acaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATAGAGTACAAGCAACTCACGTTCCCAGAAGTCAACCCAGATTTCCTAGCCATAGACACCAACCCTGTCTTCCCCCCCAGGCTGTggggg                                                                                                                                            >  1:13198/1‑150 (MQ=255)
         caacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGt                                                                                                                                           >  2:997406/1‑150 (MQ=255)
         caacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGt                                                                                                                                           >  2:997412/1‑150 (MQ=255)
          aacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTgggggaa                                                                                                                                          <  1:539866/150‑3 (MQ=255)
          aacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCGGTCTTCCCCCACATGCTGTGGGGGtt                                                                                                                                          >  2:707189/1‑150 (MQ=255)
           acGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttt                                                                                                                                         >  1:587035/1‑150 (MQ=255)
              aaTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttttt                                                                                                                                      <  2:1140499/150‑1 (MQ=255)
              aaTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttttt                                                                                                                                      <  2:1140520/150‑1 (MQ=255)
              aaTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttttt                                                                                                                                      <  2:869472/150‑1 (MQ=255)
               aTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttttttt                                                                                                                                     >  2:558219/1‑150 (MQ=255)
               aTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttttttt                                                                                                                                     >  2:194891/1‑150 (MQ=255)
               aTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttttttt                                                                                                                                     <  1:698193/150‑1 (MQ=255)
               aTGGAGCAGTTGCGAGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCCTAGACACCATCCCTCTCTGCCCCCACATGCGGTGGGGgtttgttt                                                                                                                                     >  2:940569/1‑150 (MQ=255)
                tGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttttttt                                                                                                                                    >  1:905956/1‑150 (MQ=255)
                  gatcAGTTGCGTGACGTGACGCAATGGCAAAAAATTTCATCGTTGAAGCTATTGAGTAGTTGCTACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATATACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTAt                                                                                                                                  <  1:1144483/147‑1 (MQ=255)
                  gAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTAt                                                                                                                                  <  2:584181/150‑1 (MQ=255)
                    gCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGTAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATcc                                                                                                                                >  2:598514/1‑150 (MQ=255)
                    gCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATcc                                                                                                                                >  2:598528/1‑150 (MQ=255)
                       gTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCa                                                                                                                             >  1:1387802/1‑150 (MQ=255)
                        ttGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCaa                                                                                                                            >  1:314568/1‑150 (MQ=255)
                        ttGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCaa                                                                                                                            >  1:1114002/1‑150 (MQ=255)
                           cgtgacgtgaCGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTAAAttt                                                                                                                         >  2:1267848/1‑150 (MQ=255)
                               acgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCct                                                                                                                     <  1:381416/150‑1 (MQ=255)
                                       cAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTa                                                                                                             <  1:563609/150‑1 (MQ=255)
                                       cAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTa                                                                                                             <  1:684689/150‑1 (MQ=255)
                                        aaTGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTaa                                                                                                            <  2:689362/150‑1 (MQ=255)
                                        aaTGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTaa                                                                                                            <  2:475605/150‑1 (MQ=255)
                                          tGGCAAAAAAGTTCATCGTTGGAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATg                                                                                                          <  1:598514/150‑1 (MQ=255)
                                          tGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATg                                                                                                          <  1:598528/150‑1 (MQ=255)
                                             cAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                                                                                       >  1:1192206/1‑150 (MQ=255)
                                               aaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                                                                                     <  1:950454/150‑1 (MQ=255)
                                                   aGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAga                                                                                                 <  2:96788/150‑1 (MQ=255)
                                                   aGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAga                                                                                                 <  2:21164/150‑1 (MQ=255)
                                                        aTCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGAtt                                                                                            <  2:587035/150‑1 (MQ=255)
                                                        aTCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGACATAGACACCATCCCTGTCTTCCCCCACATTCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGAtt                                                                                            <  1:1188025/150‑1 (MQ=255)
                                                         tCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCAGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGAttt                                                                                           <  2:1403911/150‑1 (MQ=255)
                                                              gAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTc                                                                                      >  1:760541/1‑150 (MQ=255)
                                                              gAAGATATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCCCATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTc                                                                                      <  1:82527/150‑1 (MQ=255)
                                                              gAAGATATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACCTGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTc                                                                                      <  1:82526/150‑1 (MQ=255)
                                                               aaGCTATTGAGTAGTAGCAACTCCCGTTCCCAGGAGTAACCCCTGTTTTCCTTGCCATAGACCCCATCCCGCTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTcc                                                                                     <  2:314856/150‑1 (MQ=255)
                                                               aaGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTcc                                                                                     <  2:1401089/150‑1 (MQ=255)
                                                                    aTTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTtat                                                                                >  2:1352132/1‑150 (MQ=255)
                                                                       gAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACt                                                                             >  2:1311413/1‑150 (MQ=255)
                                                                       gAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACt                                                                             >  1:401565/1‑150 (MQ=255)
                                                                       gAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACt                                                                             >  1:751404/1‑150 (MQ=255)
                                                                       gAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACt                                                                             >  1:751403/1‑150 (MQ=255)
                                                                        agtagtAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTa                                                                            <  2:764342/150‑1 (MQ=255)
                                                                        agtagtAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTa                                                                            <  1:385331/150‑1 (MQ=255)
                                                                             tagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTc                                                                       >  1:1287447/1‑150 (MQ=255)
                                                                             tagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTc                                                                       >  2:456570/1‑150 (MQ=255)
                                                                             tagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTc                                                                       >  1:561641/1‑150 (MQ=255)
                                                                             tagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTc                                                                       >  1:408162/1‑150 (MQ=255)
                                                                                cAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGt                                                                    <  2:401997/150‑1 (MQ=255)
                                                                                     cACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATa                                                               <  2:172637/150‑1 (MQ=255)
                                                                                     cACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATa                                                               <  2:172658/150‑1 (MQ=255)
                                                                                        gTTCCAAGTAGTAAACCCTGTTTTCCTTTCCATAGACACCATCCCTTTCTTCCCCCACATGCTGTGGTGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATATCg                                                            <  2:704946/150‑1 (MQ=255)
                                                                                         ttCCCAGTAGTAACCCCTGTTTTCCGTGCCATCGACACCCTCCCTGTCGTCCCCCACAGGCTGTGGGGGTTTTTTTTCTCCTCAATGTGCCTGCTGCTGAATGCATTGCAGATGATGTGCGTCCGTTATACTAGCGTCAGTTGATAGCgg                                                           <  2:502429/150‑1 (MQ=255)
                                                                                           cccATTAGTAAACCCTGGTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCCCCTGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGa                                                         <  1:166955/150‑1 (MQ=255)
                                                                                            ccAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAg                                                        >  2:1216677/1‑150 (MQ=255)
                                                                                            ccAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAg                                                        <  1:91687/150‑1 (MQ=255)
                                                                                             cAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGTTAGCGGGAGt                                                       >  2:1179294/1‑150 (MQ=255)
                                                                                                   taAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTAt                                                 <  2:958581/150‑1 (MQ=255)
                                                                                                   taAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGGGGGGGTTTTTTTTAGCCTCAAGTGGCCTGCTGCGTAATGCATTGCCGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGGATTGAt                                                 >  2:690844/1‑150 (MQ=255)
                                                                                                     aaCCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGa                                               >  2:784115/1‑150 (MQ=255)
                                                                                                       cccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaat                                             <  1:1074130/150‑1 (MQ=255)
                                                                                                        ccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatc                                            <  2:1020866/150‑1 (MQ=255)
                                                                                                        ccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatc                                            <  2:1020860/150‑1 (MQ=255)
                                                                                                            ttttCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatcaatc                                        <  2:905956/150‑1 (MQ=255)
                                                                                                              ttCCGTGCCATAGACACCATCCCTGTCTTCCCCCACCTGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCtt                                      <  2:130662/150‑1 (MQ=255)
                                                                                                                ccTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTAt                                    <  1:562412/150‑1 (MQ=255)
                                                                                                                 cTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATg                                   >  2:324952/1‑150 (MQ=255)
                                                                                                                 cTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATg                                   >  2:1170036/1‑150 (MQ=255)
                                                                                                                 cTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTATTATCCTCAATTTGCCTGCTGCTTAATGCCTTGCAGCTGATTTGCTTCCGGTCTACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATg                                   >  2:324886/1‑150 (MQ=255)
                                                                                                                  ttGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTTCAGATGATTTGGTTCCGGTATACTAGCGTCAGTTGATAGCGGGGGTATTTATGAATCAATCTTATgg                                  >  1:179935/1‑150 (MQ=255)
                                                                                                                  ttGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATgg                                  >  2:104268/1‑150 (MQ=255)
                                                                                                                  ttGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGGTATACTAGCGTCAGTTGATAGCGGGGGTATTTATGAATCAATCTTATgg                                  >  1:179908/1‑150 (MQ=255)
                                                                                                                  ttGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTGAATGCATTGCAGATGATTTGCGTCCGTTATACTAGCGGCAGGTGATAGCGGGAGTATTTATGAATCAATCGTATgg                                  >  1:1255481/1‑150 (MQ=255)
                                                                                                                          gACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTgg                          <  2:117704/150‑1 (MQ=255)
                                                                                                                          gACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTgg                          <  1:215806/150‑1 (MQ=255)
                                                                                                                           acacCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGt                         <  1:1086851/150‑1 (MQ=255)
                                                                                                                               cATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGt                     <  1:245807/150‑1 (MQ=255)
                                                                                                                                      aaaaacccccACAGCATGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGggcgg              >  1:1185517/6‑150 (MQ=255)
                                                                                                                                      aaaaacccccACAGCATGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCGGGTCAGTCGggcgg              >  1:1185523/6‑150 (MQ=255)
                                                                                                                                         aacccccACAGCATGTGGGGGTTTTTTTTATACTCAATTTGCCTGCTGCTTAAAGCATTGCAGATGGTTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTTTTTATGAATCAATATTATGGACGGCTGGTCAGTCGGGCGGCGa           >  1:365872/3‑150 (MQ=255)
                                                                                                                                              ccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtgct      <  2:146247/150‑1 (MQ=255)
                                                                                                                                              ccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtgct      <  1:1116088/150‑1 (MQ=255)
                                                                                                                                                  aTGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGa  <  2:401414/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGA  >  NC_000913/4106222‑4106517

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: