Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,111,596 G→A 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,5960GA97.0% 323.9 / ‑2.1 99L179L (CTG→TTG) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (48/48);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGAT  >  NC_000913_3_pae_tpiA/4111454‑4111740
                                                                                                                                              |                                                                                                                                                
ggAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCttt                                                                                                                                           >  1:1506921/1‑150 (MQ=255)
    gATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGc                                                                                                                                       <  4:369706/150‑1 (MQ=255)
    gATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGc                                                                                                                                       <  4:370416/150‑1 (MQ=255)
      tGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCg                                                                                                                                     >  1:321879/1‑150 (MQ=255)
         ccgccccgccGGGCAGCAGGTCTGTGAAACACGATAGAGTTTCCTCGGCGCAACGGCTTTGCTTTTTGTCCGTTATTCAACCCTTAAAGCGATTCAGATAGCGCCAGCTTAATCTGTTCAACAGCGAAGGTCAACCCCTTTTCGCCTTTg                                                                                                                                  >  4:141747/1‑150 (MQ=255)
         ccgccCCGCAGGGCAGCAGGTCTGTGAAACATTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTg                                                                                                                                  >  2:1349487/1‑150 (MQ=255)
         ccgccCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTg                                                                                                                                  >  3:399251/1‑150 (MQ=255)
          cgccCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGt                                                                                                                                 >  4:43589/1‑150 (MQ=255)
           gccCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTc                                                                                                                                >  2:1502040/1‑150 (MQ=255)
            ccccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTcc                                                                                                                               >  1:687894/1‑150 (MQ=255)
            ccccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTcc                                                                                                                               >  4:898461/1‑150 (MQ=255)
              ccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcgc                                                                                                                             <  2:1092878/150‑1 (MQ=255)
               cGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGtt                                                                                                                                   >  2:87075/1‑143 (MQ=255)
               cGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGtt                                                                                                                                   <  1:87075/143‑1 (MQ=255)
                gCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGa                                                                                                                           <  3:513073/150‑1 (MQ=255)
                gCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGa                                                                                                                           <  2:933613/150‑1 (MQ=255)
                   gggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaa                                                                                                                        <  3:327585/150‑1 (MQ=255)
                    ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaac                                                                                                                       <  4:399251/150‑1 (MQ=255)
                       agcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATa                                                                                                                    >  3:652613/1‑150 (MQ=255)
                           ggTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTGTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGc                                                                                                                <  3:632837/150‑1 (MQ=255)
                              cTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGt                                                                                                             <  4:343712/150‑1 (MQ=255)
                              cTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGt                                                                                                             >  4:711272/1‑150 (MQ=255)
                               tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTg                                                                                                            >  1:37066/1‑150 (MQ=255)
                               tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTg                                                                                                            >  4:53996/1‑150 (MQ=255)
                               tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTg                                                                                                            <  2:1094398/150‑1 (MQ=255)
                                gtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGc                                                                                                           <  3:898461/150‑1 (MQ=255)
                                   aaaCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAcc                                                                                                        >  3:721869/1‑150 (MQ=255)
                                    aaCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCt                                                                                                       <  3:43589/150‑1 (MQ=255)
                                     aCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCtt                                                                                                      <  4:261298/150‑1 (MQ=255)
                                       aGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                    <  2:687894/150‑1 (MQ=255)
                                         tataGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGt                                                                                                  <  1:1502040/150‑1 (MQ=255)
                                             gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCg                                                                                              <  2:1490114/150‑1 (MQ=255)
                                             gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCg                                                                                              >  1:119559/1‑150 (MQ=255)
                                              agaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCgg                                                                                             >  3:549708/1‑150 (MQ=255)
                                                aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTc                                                                                                                                         <  2:961190/104‑1 (MQ=255)
                                                aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTc                                                                                                                                         >  1:961190/1‑104 (MQ=255)
                                                    aTCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                            >  3:638661/1‑145 (MQ=255)
                                                    aTCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                            <  4:638661/145‑1 (MQ=255)
                                                       ggCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                    <  1:1501964/150‑1 (MQ=255)
                                                            aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaa                                                                                                                        >  3:586848/1‑109 (MQ=255)
                                                            aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaa                                                                                                                        <  4:586848/109‑1 (MQ=255)
                                                             aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTc                                                                                                                                <  3:565143/100‑1 (MQ=255)
                                                             aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTc                                                                                                                                >  4:565143/1‑100 (MQ=255)
                                                             aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                              >  1:504169/1‑150 (MQ=255)
                                                             aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                              <  1:1349487/150‑1 (MQ=255)
                                                                  tttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTCGTTTTTc                                                                         >  3:879394/1‑150 (MQ=255)
                                                                       tttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGcc                                                                                                                                      >  1:957353/1‑84 (MQ=255)
                                                                       tttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGcc                                                                                                                                      <  2:957353/84‑1 (MQ=255)
                                                                       tttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgc                                                                    >  3:247656/1‑150 (MQ=255)
                                                                               aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAg                                                                                                              >  2:136164/1‑100 (MQ=255)
                                                                               aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAg                                                                                                              <  1:136164/100‑1 (MQ=255)
                                                                                tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagtag                                                           >  2:1521518/1‑150 (MQ=255)
                                                                                     tCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTCGTAGCGACt                                                      >  2:1002597/1‑150 (MQ=255)
                                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                    >  1:1443904/1‑150 (MQ=255)
                                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                    >  2:1397330/1‑150 (MQ=255)
                                                                                        aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                       >  3:955753/1‑146 (MQ=255)
                                                                                        aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                       <  4:955753/146‑1 (MQ=255)
                                                                                        aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCa                                                   >  3:953965/1‑150 (MQ=255)
                                                                                        aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCa                                                   >  1:1514588/1‑150 (MQ=255)
                                                                                          ccTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGt                                                 >  4:54060/1‑150 (MQ=255)
                                                                                            ttCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttt                                               <  2:1443904/150‑1 (MQ=255)
                                                                                             tCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                              <  4:721869/150‑1 (MQ=255)
                                                                                               aaGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                            <  4:652613/150‑1 (MQ=255)
                                                                                                aGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                           <  4:953965/150‑1 (MQ=255)
                                                                                                 gCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                          >  4:410102/1‑150 (MQ=255)
                                                                                                     ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGGTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                      >  1:421342/1‑150 (MQ=255)
                                                                                                         gATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                  <  2:37066/150‑1 (MQ=255)
                                                                                                             gcgcCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                              <  3:53996/150‑1 (MQ=255)
                                                                                                               gcCAGCTTAATCGGTTCAAGATCGAAGGTCAACCCCGTTCCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                            <  3:35529/150‑1 (MQ=255)
                                                                                                                ccAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTGGCCGTTGTCCGCGACAACATAACGCAGTGCAGCTTCTGTCTCGGTGTAACAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTgg                           >  3:546577/1‑150 (MQ=255)
                                                                                                                ccAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTgg                           >  1:497370/1‑150 (MQ=255)
                                                                                                                ccAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTgg                           >  1:702685/1‑150 (MQ=255)
                                                                                                                 cAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                          <  4:174343/150‑1 (MQ=255)
                                                                                                                    cTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGa                                                        <  4:99925/117‑1 (MQ=255)
                                                                                                                    cTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGa                                                        >  3:99925/1‑117 (MQ=255)
                                                                                                                      tAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                              >  2:1494311/1‑125 (MQ=255)
                                                                                                                      tAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                              <  1:1494311/125‑1 (MQ=255)
                                                                                                                        aTCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                   <  1:1397330/150‑1 (MQ=255)
                                                                                                                         tCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAc                                                                                                         <  3:820008/63‑1 (MQ=255)
                                                                                                                         tCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAc                                                                                                         >  4:820008/1‑63 (MQ=255)
                                                                                                                          cGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg                 >  1:1495815/1‑150 (MQ=255)
                                                                                                                            gTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg               <  3:410102/150‑1 (MQ=255)
                                                                                                                             ttCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCg                                                                                              >  2:1457735/1‑70 (MQ=255)
                                                                                                                             ttCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCg                                                                                              <  1:1457735/70‑1 (MQ=255)
                                                                                                                             ttCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt              <  4:549708/150‑1 (MQ=255)
                                                                                                                                aaCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAg           >  2:158661/1‑150 (MQ=255)
                                                                                                                                 aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                           >  3:405390/1‑117 (MQ=255)
                                                                                                                                 aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                           <  4:405390/117‑1 (MQ=255)
                                                                                                                                   aGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                     >  1:1087765/1‑105 (MQ=255)
                                                                                                                                   aGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                     <  2:1087765/105‑1 (MQ=255)
                                                                                                                                    gCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGtgtg       >  4:892825/1‑150 (MQ=255)
                                                                                                                                      gAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                           <  1:808986/80‑1 (MQ=255)
                                                                                                                                      gAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg     <  2:497370/150‑1 (MQ=255)
                                                                                                                                      gAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg     <  3:412836/150‑1 (MQ=255)
                                                                                                                                      gAAGGTCAACCCCTTTACGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                           >  2:808986/1‑80 (MQ=255)
                                                                                                                                       aaGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                            >  1:604666/1‑62 (MQ=255)
                                                                                                                                       aaGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                            <  2:604666/62‑1 (MQ=255)
                                                                                                                                       aaGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCgg    <  3:669128/150‑1 (MQ=255)
                                                                                                                                         ggTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGAt  <  2:1514588/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGAT  >  NC_000913_3_pae_tpiA/4111454‑4111740

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