Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,106,370 G→C 5.4% 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.4% 451.4 / 12.1 144intergenic (+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/84);  new base C (5/3);  total (57/87)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

Only 100 of 144 total aligned reads displayed.

AAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCG  >  NC_000913/4106222‑4106516
                                                                                                                                                    |                                                                                                                                                    
aaaCATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGc                                                                                                                                                     >  1:626410/1‑150 (MQ=255)
 aaCATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCt                                                                                                                                                    >  1:349010/1‑150 (MQ=255)
     tCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCACCACATGCAGTgg                                                                                                                                                >  1:1133833/1‑150 (MQ=255)
     tCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGATCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAATCTGTgg                                                                                                                                                >  2:1523796/1‑150 (MQ=255)
      caacaaCGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTggg                                                                                                                                               <  1:1293719/150‑1 (MQ=255)
      caacaaCGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTggg                                                                                                                                               <  1:169043/150‑1 (MQ=255)
       aacaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCATGTCTTCCCCCACATGCTGTgggg                                                                                                                                              <  1:413994/150‑1 (MQ=255)
         caacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTggggga                                                                                                                                            <  2:459912/150‑2 (MQ=255)
         caacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGCAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑t                                                                                                                                          >  1:50687/1‑150 (MQ=255)
           acGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTgggggaag                                                                                                                                          <  2:197554/150‑4 (MQ=255)
            cGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑tttt                                                                                                                                       >  1:797173/1‑150 (MQ=255)
            cGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑tttt                                                                                                                                       <  1:365498/150‑1 (MQ=255)
                 ggAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑tttttgttt                                                                                                                                  >  1:1164682/1‑149 (MQ=255)
                   aGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATc                                                                                                                                <  1:1018307/150‑1 (MQ=255)
                   aGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATc                                                                                                                                <  1:1642339/150‑1 (MQ=255)
                    gCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATcc                                                                                                                               >  2:1434053/1‑150 (MQ=255)
                    gCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATcc                                                                                                                               <  2:667515/150‑1 (MQ=255)
                      cattGCGTCACGTCACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGGGGGGG‑‑TTTTTTTTATCCta                                                                                                                             >  2:1726522/3‑149 (MQ=255)
                      aGTTGCTTGACGTGACGCAATGGCAAAAAAGGTGCTCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTc                                                                                                                             <  2:324475/150‑1 (MQ=255)
                         ggCGTGACGGGACGCAATGGCAAACAAGGTCATCGGTGACGCTAGTGAGTAGTAGCAACTCCCGTTCCCAGTAGTACCCCCTGTTTTCCTTGCCCTAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTGTATCCTCAAt                                                                                                                          <  1:117054/149‑1 (MQ=255)
                          gCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAAtt                                                                                                                         >  2:744846/1‑150 (MQ=255)
                              gacgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGcc                                                                                                                     <  2:171520/150‑1 (MQ=255)
                                cgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctg                                                                                                                   <  1:1209320/150‑1 (MQ=255)
                                cgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctg                                                                                                                   >  2:133693/1‑150 (MQ=255)
                                cgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctg                                                                                                                   <  2:1060324/150‑1 (MQ=255)
                                cgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctg                                                                                                                   <  2:381016/150‑1 (MQ=255)
                                 gtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctgc                                                                                                                  >  1:1003225/1‑150 (MQ=255)
                                  tgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctgct                                                                                                                 <  1:223042/150‑1 (MQ=255)
                                  tgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctgct                                                                                                                 <  2:1517680/150‑1 (MQ=255)
                                    acgCAATGGCAAATAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCtgctgc                                                                                                               <  2:421427/150‑1 (MQ=255)
                                    aataaaTGGAAAAAAAGTCATAGGTGAACTCAATGACTTCTAGGAACAGACCTTCCCAGTACTTAACCCTTTTTACCATGCCAATGTCACCAACCATGGCTTGCCCCACCAGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTgctgct                                                                                                              <  2:1159618/146‑1 (MQ=255)
                                       cAATGGCAAAAAAGTTCATCGTTGAAGCTTTTGAGGAGTAGCAACTCACGTTCCCTGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTa                                                                                                            <  1:252972/150‑1 (MQ=255)
                                        aaTGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTaa                                                                                                           >  1:1035922/1‑150 (MQ=255)
                                           ggCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGc                                                                                                        <  1:810468/150‑1 (MQ=255)
                                             cAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                                                                                      <  1:1319338/150‑1 (MQ=255)
                                             cAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                                                                                      >  1:6159/1‑150 (MQ=255)
                                              aaaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAtt                                                                                                     <  2:1341570/150‑1 (MQ=255)
                                              aaaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAtt                                                                                                     <  2:563929/150‑1 (MQ=255)
                                               aaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                                                                                    <  2:851456/150‑1 (MQ=255)
                                                aaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGACATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGc                                                                                                   <  1:1242172/150‑1 (MQ=255)
                                                 aaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCa                                                                                                  <  1:200346/150‑1 (MQ=255)
                                                    gTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAgat                                                                                               <  1:1224163/150‑1 (MQ=255)
                                                     ttCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAgatg                                                                                              <  2:1010434/150‑1 (MQ=255)
                                                             tGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCtt                                                                                      >  2:1402394/1‑150 (MQ=255)
                                                              gAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTc                                                                                     >  2:513547/1‑150 (MQ=255)
                                                               aaGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAgatgat                                                                                            >  1:1342278/1‑142 (MQ=255)
                                                               aaGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTcc                                                                                    <  1:1626371/150‑1 (MQ=255)
                                                               aaGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGGTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCGGCTTAATGCATTGCAGATGATTTGCTTcc                                                                                    <  1:1452214/150‑1 (MQ=255)
                                                                tgCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCg                                                                                   <  1:1388490/149‑1 (MQ=255)
                                                                aGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAATCCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCg                                                                                   >  1:77904/1‑150 (MQ=255)
                                                                 gCTATTGAGTAGTAGCAACTCACGTTGCCAGTAGTAAACCCTGTTTTCCTTGCCATAGTCACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGt                                                                                  <  1:1392916/150‑1 (MQ=255)
                                                                  cTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGtt                                                                                 <  2:1704439/150‑1 (MQ=255)
                                                                   tATTGAGTAGCAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTta                                                                                <  2:50687/150‑1 (MQ=255)
                                                                    aTTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATGCTCAATTTGCCTTCTGCTTAATGCATTGCAGATGATTTGCTTACGtttt                                                                               >  1:1309140/1‑148 (MQ=255)
                                                                       gAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACt                                                                            <  2:1071514/150‑1 (MQ=255)
                                                                         gtagtaGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAg                                                                          <  2:396706/150‑1 (MQ=255)
                                                                          tagtagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGc                                                                         <  1:18091/150‑1 (MQ=255)
                                                                           agtagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCg                                                                        >  1:763471/1‑150 (MQ=255)
                                                                           agtagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TGTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCg                                                                        >  2:446303/1‑150 (MQ=255)
                                                                               gCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTc                                                                      <  1:955566/150‑1 (MQ=255)
                                                                               gCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAg                                                                    <  1:744846/150‑1 (MQ=255)
                                                                                cAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGt                                                                   <  1:1458787/150‑1 (MQ=255)
                                                                                 aaCTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGGCTGTGGGGG‑‑TTTTTTTTATACTAAATTTGCCTGGTGCTTAAAGGATTGCAGATGGTTTGGGTGCGGTTTACTAGCGTCAGtt                                                                  >  1:1543952/1‑150 (MQ=255)
                                                                                   cTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGa                                                                >  2:965951/1‑150 (MQ=255)
                                                                                    tCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGAt                                                               >  1:960280/1‑150 (MQ=255)
                                                                                    tCACGTTCCCAGTAGTAAACCCTGTGTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTAACCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGAt                                                               <  1:826709/150‑1 (MQ=255)
                                                                                     cACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATa                                                              >  2:921691/1‑150 (MQ=255)
                                                                                       cGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTGTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTGGATAGc                                                            >  2:676842/1‑150 (MQ=255)
                                                                                          tCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                         <  1:992706/150‑1 (MQ=255)
                                                                                          tCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                         <  2:506933/150‑1 (MQ=255)
                                                                                          tCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                         <  2:946635/150‑1 (MQ=255)
                                                                                          tCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                         <  2:1254067/150‑1 (MQ=255)
                                                                                           cccAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGa                                                        <  1:1327200/150‑1 (MQ=255)
                                                                                            ccAGTATTAAAGCCTTTATTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAg                                                       <  1:113209/150‑1 (MQ=255)
                                                                                                tagtaAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTATTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTAtt                                                   <  1:438722/150‑1 (MQ=255)
                                                                                                 atgaaaCCCTGTTTTCCTTGCCATAGCCACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTAttt                                                  <  2:771226/147‑1 (MQ=255)
                                                                                                  gtaAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTTATAGCGGTAGTATTTa                                                 >  1:798139/1‑150 (MQ=255)
                                                                                                    aaaCCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATg                                               <  1:390163/150‑1 (MQ=255)
                                                                                                        ccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTACTGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatc                                           <  2:1477344/150‑1 (MQ=255)
                                                                                                            ttttCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatcaatc                                       <  2:1543122/150‑1 (MQ=255)
                                                                                                            tgttCCTTGCCATAGACAGCATCCCTGTCTTCCCACACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatcaatc                                       <  2:36480/148‑1 (MQ=255)
                                                                                                             tttCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCt                                      <  1:1330233/150‑1 (MQ=255)
                                                                                                              ttCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCtt                                     >  1:18169/1‑150 (MQ=255)
                                                                                                                 cTTGCCATAGACACCATCCCAGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATg                                  <  2:760466/150‑1 (MQ=255)
                                                                                                                      cATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACgg                             >  1:845361/1‑150 (MQ=255)
                                                                                                                      cATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACgg                             >  2:1639650/1‑150 (MQ=255)
                                                                                                                      cATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATACTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACgg                             >  1:356293/1‑150 (MQ=255)
                                                                                                                       aTAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGc                            <  1:392069/150‑1 (MQ=255)
                                                                                                                        tAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCt                           <  2:1646635/150‑1 (MQ=255)
                                                                                                                          gaaacaaTCCCTGTCTTCCGACACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTgg                         <  2:1706033/144‑1 (MQ=255)
                                                                                                                               cATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGt                    <  2:218185/150‑1 (MQ=255)
                                                                                                                                        cTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGgcggcg           >  2:797766/1‑150 (MQ=255)
                                                                                                                                        cTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGgcggcg           >  2:1004309/1‑150 (MQ=255)
                                                                                                                                        cTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGgcggcg           >  2:1344433/1‑150 (MQ=255)
                                                                                                                                        cTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCATAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGgcggcg           >  2:465150/1‑150 (MQ=255)
                                                                                                                                         aacccccACAGCATGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGa          >  2:1229883/3‑150 (MQ=255)
                                                                                                                                          acccccACAGCATGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGAt         >  2:889989/2‑150 (MQ=255)
                                                                                                                                           cccccACATGCAGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGAtt        <  1:866277/150‑1 (MQ=255)
                                                                                                                                           cccccACAGCATGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGAtt        >  1:513256/1‑150 (MQ=255)
                                                                                                                                                acaTGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtgctgc   >  1:1405148/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCG  >  NC_000913/4106222‑4106516

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