Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,106,369 T→G 6.6% 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 6.6% 394.1 / 12.8 144intergenic (+49/‑100)cpxP/fieFinhibitor of the cpx response; periplasmic adaptor protein/ferrous iron and zinc transporter
Reads supporting (aligned to +/- strand):  ref base T (49/83);  new base G (7/3);  total (57/87)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 145 total aligned reads displayed.

AGAAACATCAACAACGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCG  >  NC_000913/4106220‑4106516
                                                                                                                                                          |                                                                                                                                                     
agaAACATCAACAACGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGATGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGAAGTAAACCCTGTTTTCCTAGCCATAGACACCATACCTGTCTTCACCCACAt                                                                                                                                                       >  1:36480/1‑150 (MQ=255)
  aaaCATCAACAACGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGc                                                                                                                                                     >  1:626410/1‑150 (MQ=255)
   aaCATCAACAACGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCt                                                                                                                                                    >  1:349010/1‑150 (MQ=255)
      aTCAACAACGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCtgtg                                                                                                                                                 >  2:427507/1‑150 (MQ=255)
       tCAACAACGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCACCACATGCAGTgg                                                                                                                                                >  1:1133833/1‑150 (MQ=255)
       tCAACAACGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGATCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAATCTGTgg                                                                                                                                                >  2:1523796/1‑150 (MQ=255)
        caacaaCGAATGGAGCAGTTGCGTGACGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCt                                                                                                                                                    >  2:18091/1‑150 (MQ=255)
        caacaaCGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTggg                                                                                                                                               <  1:886720/150‑1 (MQ=255)
        caacaaCGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTggg                                                                                                                                               <  1:1293719/150‑1 (MQ=255)
        caacaaCGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTggg                                                                                                                                               <  1:169043/150‑1 (MQ=255)
           caacGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCCTCCCTGTCTTCCCCCACATGGGGTgggggg                                                                                                                                            >  1:1543122/1‑149 (MQ=255)
           caacGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTggggga                                                                                                                                            <  2:459912/150‑2 (MQ=255)
              cGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑tttt                                                                                                                                       <  1:365498/150‑1 (MQ=255)
                  tGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑tttttttt                                                                                                                                   <  2:171316/150‑1 (MQ=255)
                   ggAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑tttttgttt                                                                                                                                  >  1:1164682/1‑149 (MQ=255)
                     aGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATc                                                                                                                                <  1:1642339/150‑1 (MQ=255)
                     aGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATc                                                                                                                                <  1:1018307/150‑1 (MQ=255)
                      gCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATcc                                                                                                                               <  2:667515/150‑1 (MQ=255)
                        cattG‑‑‑‑‑CGTCACGTCACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGGGGGGG‑‑TTTTTTTTATCCta                                                                                                                             >  2:1726522/3‑149 (MQ=255)
                        aGTTG‑‑‑‑‑CTTGACGTGACGCAATGGCAAAAAAGGTGCTCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTc                                                                                                                             <  2:324475/150‑1 (MQ=255)
                           gg‑‑‑‑‑CGTGACGGGACGCAATGGCAAACAAGGTCATCGGTGACGCTAGTGAGTAGTAGCAACTCCCGTTCCCAGTAGTACCCCCTGTTTTCCTTGCCCTAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTGTATCCTCAAt                                                                                                                          <  1:117054/149‑1 (MQ=255)
                                     gacgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGcc                                                                                                                     <  1:795800/150‑1 (MQ=255)
                                     gacgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGcc                                                                                                                     <  2:1409686/150‑1 (MQ=255)
                                       cgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctg                                                                                                                   <  2:381016/150‑1 (MQ=255)
                                       cgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctg                                                                                                                   >  2:133693/1‑150 (MQ=255)
                                       cgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctg                                                                                                                   <  1:1209320/150‑1 (MQ=255)
                                        gtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctgc                                                                                                                  >  1:1003225/1‑150 (MQ=255)
                                         tgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCctgct                                                                                                                 <  1:223042/150‑1 (MQ=255)
                                           acgCAATGGCAAATAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCtgctgc                                                                                                               <  2:421427/150‑1 (MQ=255)
                                             gCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCtt                                                                                                             <  2:744535/150‑1 (MQ=255)
                                               aaTGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTaa                                                                                                           >  1:1035922/1‑150 (MQ=255)
                                                    cAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                                                                                      >  2:1129434/1‑150 (MQ=255)
                                                    cAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                                                                                      <  1:1319338/150‑1 (MQ=255)
                                                     aaaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAtt                                                                                                     <  2:563929/150‑1 (MQ=255)
                                                     aaaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAtt                                                                                                     <  2:1341570/150‑1 (MQ=255)
                                                      aaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                                                                                    <  2:851456/150‑1 (MQ=255)
                                                       aaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGc                                                                                                   <  1:404453/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‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAgatgat                                                                                            <  2:1342278/142‑1 (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)
                                                                          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)
                                                                                      gCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTc                                                                      <  1:955566/150‑1 (MQ=255)
                                                                                      gCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAg                                                                    <  1:744846/150‑1 (MQ=255)
                                                                                        aaCTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGGCTGTGGGGG‑‑TTTTTTTTATACTAAATTTGCCTGGTGCTTAAAGGATTGCAGATGGTTTGGGTGCGGTTTACTAGCGTCAGtt                                                                  >  1:1543952/1‑150 (MQ=255)
                                                                                         aCTCACGTTCCCAGTAGTAAACCCTTTTTTCCTTGCCATAGACACCATCACTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGACTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTg                                                                 <  1:836042/150‑1 (MQ=255)
                                                                                          cTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGa                                                                >  2:965951/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:1254067/150‑1 (MQ=255)
                                                                                                 tCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                         <  2:946635/150‑1 (MQ=255)
                                                                                                 tCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                         <  2:506933/150‑1 (MQ=255)
                                                                                                  cccAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGa                                                        <  2:1498234/150‑1 (MQ=255)
                                                                                                  cccAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGa                                                        <  1:1327200/150‑1 (MQ=255)
                                                                                                  acctgtagtaAACCCTGTTTTCCTTGGCATAGACACAATCACTTTCTTCCCCCTCATGCTGTGTGGG‑‑TTTTTTTTATCCTCAATTTGGCTGATGCTTAATGCATTGAAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGa                                                        <  1:450849/146‑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)
                                                                                                        agtaAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTAttt                                                  >  1:218185/1‑150 (MQ=255)
                                                                                                         gtaAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTTATAGCGGTAGTATTTa                                                 >  1:798139/1‑150 (MQ=255)
                                                                                                               ccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTACTGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatc                                           <  2:1477344/150‑1 (MQ=255)
                                                                                                               ccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatc                                           <  2:1045564/150‑1 (MQ=255)
                                                                                                                  gTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatcaat                                        >  1:574794/1‑150 (MQ=255)
                                                                                                                  gTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCATCAATTTGACTGCTGCTTAATGGATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTTTGaataaat                                        >  2:169141/1‑150 (MQ=255)
                                                                                                                   ttttCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatcaatc                                       <  2:1543122/150‑1 (MQ=255)
                                                                                                                   tgttCCTTGCCATAGACAGCATCCCTGTCTTCCCACACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatcaatc                                       <  2:36480/148‑1 (MQ=255)
                                                                                                                    tttCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCt                                      <  2:1003225/150‑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                             >  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)
                                                                                                                                   cacCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTc                       >  1:673618/1‑150 (MQ=255)
                                                                                                                                      cATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGt                    <  2:218185/150‑1 (MQ=255)
                                                                                                                                            aaaaaacccccACAGCATGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGggcg              >  1:295041/7‑150 (MQ=255)
                                                                                                                                               cTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGgcggcg           >  2:797766/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)
                                                                                                                                                     ccACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGCAGTATTAATGAATCAATCTAATGTACGGCTGGACAGTCGGGCTGCGATtgat     >  1:1244612/1‑148 (MQ=255)
                                                                                                                                                        caTGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCg  <  1:800483/150‑1 (MQ=255)
                                                                                                                                                        caTGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCg  >  1:1410251/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                     
AGAAACATCAACAACGAATGGAGCAGTTG‑‑‑‑‑CGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGG‑‑TTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCG  >  NC_000913/4106220‑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: