Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 4,111,596 G→A 100% L179L (CTG→TTG)  yiiQ ← DUF1454 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,111,5960GA94.6% 175.4 / ‑4.1 70L179L (CTG→TTG) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (37/24);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCG  >  NC_000913_3_pae_tpiA/4111454‑4111737
                                                                                                                                              |                                                                                                                                             
ggAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTTATTTTTTCAAACCTTCGAGCGCATAAGATAGCACCACCGTAATCGGTTCAACAGCGAGGGTCAACCCCtt                                                                                                                                         >  2:733177/1‑149 (MQ=255)
    gATGGCCGCCCCGCAGGGCAGCAGGTCTGCGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTNAACAGCGAAGGTCAACCCCTTTTCg                                                                                                                                     >  1:615184/1‑149 (MQ=255)
         ccgccCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTGATTTATTCAAACCTTAAAGCGATTCAGGTAGCGCCAGCTTAATCGGTTCAGCAGAGCAGGCCAACCCCTTTTCGCCCtt                                                                                                                                >  2:554041/1‑149 (MQ=255)
         ccgccCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGCAGGTCCGCCCGTGTTCTCCgat                                                                                                                                >  2:554033/1‑147 (MQ=255)
          cgccCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTTTTTTTTTTTAACCTTCAAGCGATTCGGATAGCGCCAGCTTAATCGGTTTAACAGCGAAGGTTACCCCCTTTTCGCCGTTg                                                                                                                               >  2:782778/1‑149 (MQ=255)
              ccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTAATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcg                                                                                                                           >  2:859064/1‑149 (MQ=255)
                   gggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGaca                                                                                                                      >  1:861941/1‑149 (MQ=255)
                    ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaa                                                                                                                     <  1:152109/149‑1 (MQ=255)
                      cagcagGTCTGTGAAACAGGATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaaca                                                                                                                   <  1:212077/149‑1 (MQ=255)
                        gcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATa                                                                                                                 >  1:760992/1‑149 (MQ=255)
                        gcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCCTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATa                                                                                                                 <  2:6710/149‑1 (MQ=255)
                          agGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAAc                                                                                                               >  1:440222/1‑149 (MQ=255)
                           ggTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCACCCGCAGACCTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACg                                                                                                              <  1:689429/149‑1 (MQ=255)
                            gTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGc                                                                                                             >  1:37717/1‑149 (MQ=255)
                             tCTTTGGAAGAGTATGGGGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTGTTCTCTCCGTCGAGCGATTCAGATAGCGCTAGCTTGATCGGTTCAACAGCGAAGGTCAGCGGCTCTTCGCCGTTGTCCGCGACAACATAACGCa                                                                                                            <  2:7933/149‑1 (MQ=255)
                               tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCTGATGGCGGCCGCAAAAGCTGGTCAACGGAGAAAGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGt                                                                                                          <  2:662261/149‑1 (MQ=255)
                                 tgAAACAGTATAGAGATTGATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGATCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGc                                                                                                        >  2:810572/1‑149 (MQ=255)
                                        gTATAGAGAGTCATCGGCACAAAGGNTTTGCTTTTTGTCATTTATTCCAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTTGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                 <  1:619192/149‑1 (MQ=255)
                                               gaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCgg                                                                                          <  1:409187/149‑1 (MQ=255)
                                                aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                         <  1:580612/149‑1 (MQ=255)
                                                    ttaaaaacaAAGGTTTATATTTTTGTCTTTTATACTCACTTTCAATCGATTCTGATAGCGTCAGCTTAATCGGTTCAACAGCGAAGTTCAACCTCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaa                                                                                     <  2:646537/143‑1 (MQ=255)
                                                     tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGAAGTGCACCTTCTGTCTCGGTGtaat                                                                                    >  2:314411/1‑149 (MQ=255)
                                                          acaAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgt                                                                               >  1:642872/1‑149 (MQ=255)
                                                             aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGAAATAACgtttg                                                                            >  1:514493/1‑149 (MQ=255)
                                                                  tttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGAGCCAGCTTAATCGGTTCAACAGCGAAGGTAAACCCCTTTTCGCCTTTGTCCGCGCCAACATAACGCAGTGCAACTTCTGTCTCGGTATAATAACATTTGttttt                                                                       >  2:47429/1‑149 (MQ=255)
                                                                   ttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAAGCCCTTTTCGCCGTTGTCCGCGACAACATCACGCAGTGCACCTTCTGTTTTGGTGATATAACGCTTGTTTTTc                                                                      >  1:844058/1‑149 (MQ=255)
                                                                        ttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgc                                                                 <  1:524152/149‑1 (MQ=255)
                                                                             tCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACCACATCACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTta                                                            <  1:510074/149‑1 (MQ=255)
                                                                               aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCTCCGTTTTTCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTTCCCGCCGTtagt                                                          >  2:294250/1‑149 (MQ=255)
                                                                                tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACCTAACGCAGTGAACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCACGCCGTtcgta                                                         >  2:568507/1‑149 (MQ=255)
                                                                                tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtattc                                                         >  2:568503/1‑146 (MQ=255)
                                                                                tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTCTTTCCCCGCCGTtagta                                                         >  2:856330/1‑149 (MQ=255)
                                                                                    ttaaaaCCTTCAAGCGATTCATGTACCGCTAGTCTATTCGTTTCATCAGCCAATCTCACCTTCCTTTCGCTCTGGTCCGCGACAACATATCGCTGTGCCCCCTCTGTCTCCGTCTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGa                                                     <  2:332780/146‑1 (MQ=255)
                                                                                    ttaaaaCCTTAAATCGATTCATATAGCGCAAGTTTATTCGTTTTAACATCTAAGGTCAACCCCATTCCGCTGTTGTCGCCAAAAACATAACGCTTTGCCCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGa                                                     <  2:871009/146‑1 (MQ=255)
                                                                                    ttCAAACCTTCAAGCGATTCCAATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCACCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGa                                                     <  1:521979/149‑1 (MQ=255)
                                                                                    ttCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTACTAACGTTTGTTTTTCCCCGCCGTTAGTAGCGa                                                     >  2:819843/1‑149 (MQ=255)
                                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCTAAGGTCAACCCCTTTTCGCCGTTGCCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTCGTAGCGACTg                                                  >  1:429277/1‑149 (MQ=255)
                                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTGGCGCCTg                                                  >  1:528808/1‑149 (MQ=255)
                                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAAATGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                  >  2:57656/1‑149 (MQ=255)
                                                                                        aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                 >  2:547167/1‑149 (MQ=255)
                                                                                           cTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGt                                              >  1:888196/1‑149 (MQ=255)
                                                                                           cTTCAAGCGATTCAGATAGAGCAAGATTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTGCTCCCTGCCGTTACGAGCGACTGAAGt                                              >  2:806990/1‑149 (MQ=255)
                                                                                                     tttaGATAGCTCCAGCTTTATCGGTTCAACAGCGTAGGTTTTTCTCTTTTTGTCCTTGTCTGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTCTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                    <  2:223128/146‑1 (MQ=255)
                                                                                                     ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTCGCGCCTGCAGTTTTTTCTGGc                                    >  1:257444/1‑149 (MQ=255)
                                                                                                      tCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCACCGCAGTTAGTAGCTACTCCAGTTTTTTCTAGcc                                   >  2:375492/1‑148 (MQ=255)
                                                                                                          aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                               >  2:229284/1‑149 (MQ=255)
                                                                                                              gggcAGCTGAATCGGTTTGACTCGGCAGGTCAGCTCTTTTTTGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTGTGTCTCGGTGTAATGACGTTTGTTTTTGCCCGCCGTTAGTAGCGCCTGCAGTTTGTTCTGGCATCGTGttt                           <  2:690730/146‑1 (MQ=255)
                                                                                                              cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGGAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                           >  2:715663/1‑149 (MQ=255)
                                                                                                               gtgaGCTTAATCGGTTCAACATCGGATGTCAGCCTTTTCTCTCTGTTGTCCGCGACAACATACCTCCGTGTCCCTTCTGTCTCGGTGTAATAACGTTAGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                          <  2:42572/146‑1 (MQ=255)
                                                                                                                    ctgaaTCGGTTTAAAGGGGAAGGTGCACCATTTTTCGCCGTTGTCCGCGTCAACATAACGCTGTGCCCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                     <  2:670012/146‑1 (MQ=255)
                                                                                                                    cTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtctc                                                                                            <  2:455340/78‑1 (MQ=255)
                                                                                                                    cTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtctc                                                                                            >  1:455340/1‑78 (MQ=255)
                                                                                                                    cTTAATCGGTTCAACAGCGAACGCCAACCCCTGTTCTTAGTCGTCTGCTTTATCTTAACGCAGTGCACCTTCGGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                     <  1:677793/149‑1 (MQ=255)
                                                                                                                     ttAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTACTAACGTTTGTTTTTCCCCGCCGTTAATTGCGACTCCAGCTTTTTCGGGCGCAGTGTTTTGCTCag                    >  2:745954/1‑148 (MQ=255)
                                                                                                                     ttAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAACTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGACACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                    >  1:602694/1‑149 (MQ=255)
                                                                                                                     ttAATCGGTTCAACAGCGAAGGTCAACCCCTTGTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTATGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                    <  1:50760/149‑1 (MQ=255)
                                                                                                                      tAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                   >  1:898497/1‑149 (MQ=255)
                                                                                                                        aTCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                 <  1:755218/149‑1 (MQ=255)
                                                                                                                         tCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCAGTGACAACATAACGCCGGGACGCAGCTGGCTCGTTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGtcattcat                >  2:662705/1‑149 (MQ=255)
                                                                                                                          cgtcacAACTGCTAAGTTCAATCCCTTTTCACCGTTGTCTGTTACAACATAATGCAGTGCACCTTCTGTCTCGGTGTCATAAAGTTTATTTTTCCCCGCAGTTAGTGGCGACTCCTGTTCTTTCTTGCCTTGTGTTTTGGTCATTAATg               <  2:560106/144‑1 (MQ=255)
                                                                                                                           gtgtCAACAGCGCAGGTCAACCCCTTTCCGCCGTTGCCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg              <  1:500982/146‑1 (MQ=255)
                                                                                                                            gTTCAACAGCGAAGGTCAACCCCCTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg             <  1:782778/149‑1 (MQ=255)
                                                                                                                             ttCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACATTCTGTCTCAGTGTAATAACGTTTGTTTTTCCCTGCCGTTATTAGCTACTGCAGTTGGTTATGGCTCTGTGGTTTGGTCATTAATgggg            >  2:226599/1‑149 (MQ=255)
                                                                                                                                 aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAg        >  2:454192/1‑149 (MQ=255)
                                                                                                                                  catagAAGGTCAACCCTTTTTCTCCGTTGCCCGCGACAACATAACGCAGTGCACTTTCTGTCTTTGTGTAATAACGTTTTTTTTTCTCCTCCGTTAGTAGCGACTGCATTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt       <  1:180839/145‑1 (MQ=255)
                                                                                                                                  cAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaataa                                                                                  <  2:250466/74‑1 (MQ=255)
                                                                                                                                  cAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaataa                                                                                  >  1:250466/1‑74 (MQ=255)
                                                                                                                                  cAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt       <  1:454192/149‑1 (MQ=255)
                                                                                                                                      gAAGGACAGTCTCTTTGCGCCGTTGTCCGCGACAACATATCGCAGTGCACCTTCTATCCCGGTATAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc   <  2:671938/149‑1 (MQ=255)
                                                                                                                                       aaGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCTGTTATTTGCTACTGCAGTGTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg  <  1:465427/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCG  >  NC_000913_3_pae_tpiA/4111454‑4111737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑