Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,106,370 G→C 5.8% 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.8% 586.4 / 21.8 172intergenic (+50/‑99)cpxP/fieFinhibitor of the cpx response; periplasmic adaptor protein/ferrous iron and zinc transporter
Reads supporting (aligned to +/- strand):  ref base G (80/81);  new base C (2/8);  total (82/90)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

Only 100 of 172 total aligned reads displayed.

ACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGA  >  NC_000913/4106224‑4106517
                                                                                                                                                  |                                                                                                                                                     
    caacaaCGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTgg                                                                                                                                                <  1:406519/149‑1 (MQ=255)
       caacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTggggg                                                                                                                                             <  1:452339/149‑1 (MQ=255)
       caacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTgggg                                                                                                                                              <  3:145757/148‑1 (MQ=255)
        aacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGt                                                                                                                                           >  1:821442/1‑149 (MQ=255)
        aacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTggggga                                                                                                                                            <  1:242010/149‑2 (MQ=255)
        aacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTggggg                                                                                                                                             <  4:542718/148‑1 (MQ=255)
         acGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTggggga                                                                                                                                            <  4:479867/148‑2 (MQ=255)
          cGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtt                                                                                                                                          <  3:275154/148‑1 (MQ=255)
          cGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACAGCATGTgggggaag                                                                                                                                          <  1:965507/149‑4 (MQ=255)
           gAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttt                                                                                                                                        <  2:750223/149‑1 (MQ=255)
            aaTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttttt                                                                                                                                       <  1:175266/149‑1 (MQ=255)
            aaTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttt                                                                                                                                        <  3:526675/148‑1 (MQ=255)
             aTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGATTTCCTTGCCATAGACACCATCCCTGTCTT‑‑CCCACATGCTGTGGGGGttttttt                                                                                                                                     >  4:264236/1‑148 (MQ=255)
              tGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttttt                                                                                                                                      >  3:211879/1‑148 (MQ=255)
                   cAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTAt                                                                                                                                  >  4:757031/1‑147 (MQ=255)
                      ttGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTAt                                                                                                                                  <  3:842976/144‑1 (MQ=255)
                      ttGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTAt                                                                                                                                  >  4:842976/1‑144 (MQ=255)
                       tGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCaa                                                                                                                            <  2:26993/149‑1 (MQ=255)
                       tGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCa                                                                                                                             <  4:641514/148‑1 (MQ=255)
                            gacgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGc                                                                                                                       <  1:546602/149‑1 (MQ=255)
                              cgtgacgGAATGGCAAAAAAGTTCATCCTTGAAGCTATTGACTAGTAGCAACCCACGTTCCAAGTAGTATACGCTCTTATCCTTTCTATAGACACCATCCCAGTCTTCCCCCACATGCTGTCGGGGTTTTTATTATCCTCAATTTGCct                                                                                                                     <  2:773393/149‑1 (MQ=255)
                              cgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGcc                                                                                                                      >  3:379278/1‑148 (MQ=255)
                                 gacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctgct                                                                                                                  <  2:953315/149‑1 (MQ=255)
                                     cAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTgctgct                                                                                                               >  4:473825/1‑148 (MQ=255)
                                      aaTGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTa                                                                                                             >  2:823816/1‑149 (MQ=255)
                                        tGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAAt                                                                                                           <  1:627984/149‑1 (MQ=255)
                                         gccAAAAAAGTTCATCGTTAAAGCTATTAATAATTGGCACCTCCCTTTCCCAGTAGTAACCCCGGTTTTCCTTGCAATAGACCCCTCCCCTGCCTCCCCCCCAATGCTGGGGGGGTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGc                                                                                                         <  2:821442/147‑1 (MQ=255)
                                            aaaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                                                                                       >  2:518225/1‑149 (MQ=255)
                                               aaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                                                                                     <  3:578225/148‑1 (MQ=255)
                                               aaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                                                                                     >  4:118647/1‑148 (MQ=255)
                                               aaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGc                                                                                                    <  2:1031777/149‑1 (MQ=255)
                                                   ttCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCTTCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTCCCCGCTTCTTAATGCATCGCCgat                                                                                                >  2:842300/1‑149 (MQ=255)
                                                      aTCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAgatgat                                                                                             >  1:361742/1‑149 (MQ=255)
                                                        cGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGAttt                                                                                           <  1:248891/149‑1 (MQ=255)
                                                        cGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGAtt                                                                                            >  3:598597/1‑148 (MQ=255)
                                                            gAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCtt                                                                                       >  1:800913/1‑149 (MQ=255)
                                                            gAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCtt                                                                                       >  2:571636/1‑149 (MQ=255)
                                                                cTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCg                                                                                    <  4:211879/148‑1 (MQ=255)
                                                                cTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCg                                                                                    <  3:149188/148‑1 (MQ=255)
                                                                 tATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGtt                                                                                  >  2:372478/1‑149 (MQ=255)
                                                                  aTTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTta                                                                                 >  1:607196/1‑149 (MQ=255)
                                                                   ttGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttt                                                                                                                                        <  1:116419/93‑1 (MQ=255)
                                                                   ttGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttt                                                                                                                                        >  2:116419/1‑93 (MQ=255)
                                                                    tGAGTAGTAGCAACTCACGTTGCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCTTCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCGCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCATtat                                                                                >  4:660566/1‑148 (MQ=255)
                                                                    tGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTtat                                                                                <  4:597428/148‑1 (MQ=255)
                                                                     gAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATAc                                                                              <  2:326581/149‑1 (MQ=255)
                                                                       gtagtaGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACt                                                                             <  3:118647/148‑1 (MQ=255)
                                                                       gtagtaGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTa                                                                            >  2:97998/1‑149 (MQ=255)
                                                                        tagtagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGc                                                                                                                       >  1:966251/1‑105 (MQ=255)
                                                                        tagtagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAg                                                                           <  1:823816/149‑1 (MQ=255)
                                                                           tagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAtt                                                                                                      >  1:1038419/1‑119 (MQ=255)
                                                                           tagCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAtt                                                                                                      <  2:1038419/119‑1 (MQ=255)
                                                                             gCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTc                                                                       >  3:466033/1‑148 (MQ=255)
                                                                               aaCTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGt                                                                    >  2:1790/1‑149 (MQ=255)
                                                                                aCTCGCGTTCCCCGTTGTAAAACCTTTTTTTCTTGCCATAGACACCCACCCTGTCTTCCCCCACCTGATGTGGGGGTTTTTTTTTCCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGtt                                                                   <  1:842300/149‑1 (MQ=255)
                                                                                    aCGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTGGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATa                                                               <  2:123156/149‑1 (MQ=255)
                                                                                       ttCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGc                                                             >  4:218011/1‑148 (MQ=255)
                                                                                       ttCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGc                                                             >  3:843391/1‑148 (MQ=255)
                                                                                        tCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCg                                                            >  4:763667/1‑148 (MQ=255)
                                                                                          ccAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAAt                                                                                                           >  1:56310/1‑99 (MQ=255)
                                                                                          ccAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGa                                                         <  1:22931/149‑1 (MQ=255)
                                                                                             gtagtaAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGt                                                       >  3:1055700/1‑148 (MQ=255)
                                                                                                gtaAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTAttt                                                   >  2:974858/1‑149 (MQ=255)
                                                                                                 taAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTa                                                  <  2:46592/149‑1 (MQ=255)
                                                                                                  aaaCCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTa                                                  <  4:394204/148‑1 (MQ=255)
                                                                                                   aaCCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTAt                                                 <  2:480639/149‑1 (MQ=255)
                                                                                                      ccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTTCTGCTTAATGCATTGCAGATGATTTGCTTCTGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaa                                              <  3:405527/148‑1 (MQ=255)
                                                                                                      ccTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaat                                             >  1:717907/1‑149 (MQ=255)
                                                                                                       cTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatc                                            >  1:335201/1‑149 (MQ=255)
                                                                                                        tGTTTTCCTTGCCATAGACACCATCCCTGTCTT‑‑CCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatcaa                                          <  3:264236/148‑1 (MQ=255)
                                                                                                          ttttCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTAt                                                                                                                                  <  4:232770/60‑1 (MQ=255)
                                                                                                          ttttCCGTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGaatcaa                                          >  4:958652/1‑148 (MQ=255)
                                                                                                               cTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTAt                                    <  1:472151/149‑1 (MQ=255)
                                                                                                                ttGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTAt                                    <  4:147521/148‑1 (MQ=255)
                                                                                                                ttGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATg                                   >  1:52294/1‑149 (MQ=255)
                                                                                                                ttGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATg                                   >  1:442594/1‑149 (MQ=255)
                                                                                                                  gCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGa                                 <  2:204644/149‑1 (MQ=255)
                                                                                                                   ccATAGACGCCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGACTGCTGCTTAATGCATTGCAGATTATTTGCTTCCGTTATACTCGCGTCAGCTGATAGCGGGAGTATTTATGAATCAATCTTATGGa                                 <  3:660566/148‑1 (MQ=255)
                                                                                                                   ccATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGa                                 <  4:13299/148‑1 (MQ=255)
                                                                                                                        gACACCATCCCTGTCGTCCCCCTCATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCATGCTGCTTAATGCATTGCAGATGATTTGCTTCCGCTATACGAGCGTCAGTTGATAGCGGGAGTATTGATGATTCAATCTTATGGACGGCTg                           <  2:135473/149‑1 (MQ=255)
                                                                                                                            ccATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTc                        >  3:809895/1‑148 (MQ=255)
                                                                                                                            ccATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTc                        >  3:155009/1‑148 (MQ=255)
                                                                                                                            ccATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTCCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTc                        >  4:792152/1‑148 (MQ=255)
                                                                                                                              aTCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                                                                                     <  4:1033801/69‑1 (MQ=255)
                                                                                                                              aTCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                                                                                     >  3:1033801/1‑69 (MQ=255)
                                                                                                                               tCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGt                     >  3:425476/1‑148 (MQ=255)
                                                                                                                                cccTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCg                   >  1:713564/1‑149 (MQ=255)
                                                                                                                                     tCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGggcgg              >  2:125029/1‑149 (MQ=255)
                                                                                                                                        tCCCCCACATGCTATGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGTAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGCTAGCGGAAGTATTTATAAATCAATCTTATGGACGCCTGGTCAGTCGGGCGGCGa           >  2:409356/1‑149 (MQ=255)
                                                                                                                                         cccccACAGCATGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGa                                               <  3:922379/112‑1 (MQ=255)
                                                                                                                                           cccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtg        <  2:212904/149‑1 (MQ=255)
                                                                                                                                           cccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtg        <  1:125029/149‑1 (MQ=255)
                                                                                                                                            ccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCtt                                      <  1:1000971/118‑1 (MQ=255)
                                                                                                                                            ccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCtt                                      >  2:1000971/1‑118 (MQ=255)
                                                                                                                                            ccACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtgc       >  2:663405/1‑149 (MQ=255)
                                                                                                                                             cacaTGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATtgc       >  3:854879/1‑148 (MQ=255)
                                                                                                                                              acaTGCTGTGGGGGTTTTTTTTATCCCCAATTTGCCTGCTGCTTAATGCATTGCAGATGACTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATAGACGGCTGGTCAGTCGGGCGGCCAGtgctg     <  2:733330/149‑1 (MQ=255)
                                                                                                                                               caTGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCATATGGACGGCTGGTCAGTCGGTCGGCGATtgctg     >  3:190575/1‑148 (MQ=255)
                                                                                                                                                 tGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                          >  4:1002842/1‑93 (MQ=255)
                                                                                                                                                 tGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCggg                                                          <  3:1002842/93‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTATGAATCAATCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGA  >  NC_000913/4106224‑4106517

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