Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 372,055 C→A 6.8% H273N (CAT→AAT)  mhpC → 2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA372,0550CA 6.8% 152.6 / 3.7 60H273N (CAT→AAT) mhpC2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base A (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-02

CTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGA  >  NC_000913_3_pae_tpiA/371908‑372197
                                                                                                                                                    |                                                                                                                                              
cTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCNACTGTGGTCACTGGGCGCAGTGGGAACa                                                                                                                                               >  1:348579/1‑149 (MQ=255)
cTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACa                                                                                                                                               <  1:435014/149‑1 (MQ=255)
   gCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCTGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGc                                                                                                                                            <  1:244799/149‑1 (MQ=255)
     ggAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCg                                                                                                                                          >  2:619233/1‑149 (MQ=255)
     ggAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATAACTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCg                                                                                                                                          >  1:28557/1‑149 (MQ=255)
           cAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCAGGCGCtt                                                                                                                                    >  1:245379/1‑149 (MQ=255)
            aaaGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCttt                                                                                                                                   >  1:147547/1‑149 (MQ=255)
                   aaaCCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGACGACGCTTtcaatcaa                                                                                                                            >  2:914332/1‑149 (MQ=255)
                      cccTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTg                                                                                                                         >  2:96330/1‑149 (MQ=255)
                       ccTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTgg                                                                                                                        >  2:763589/1‑149 (MQ=255)
                          catcgTCTGGGGGCGCGATGCACGCCTTGTGCAGAGGGATGCTGGTCTGAGTCTTTCATCCGGCATTGCCGGTTCTGAACTGCATAACTGCCAAGACTATGGTCACTGAGCGAAGTGGGAACATGCCGACGCTTTCAATCAACTGGTga                                                                                                                     <  2:124872/145‑2 (MQ=255)
                                    gggCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCt                                                                                                           <  1:97221/149‑1 (MQ=255)
                                    gggCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCt                                                                                                           <  1:344955/149‑1 (MQ=255)
                                      gcgcAACGTCCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCg                                                                                                         <  1:334497/149‑1 (MQ=255)
                                      gcgcAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCg                                                                                                         <  1:559587/149‑1 (MQ=255)
                                                 gcatcgtgCCGAGGGATGCGGGTCTGCGGCTGATGTCCGGCATTGCCGGTTCGGAACTGCATATCTGCCGCGACTGTGGTCACTGGGTGCAGTGTGAACAGGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTaa                                                                                              <  2:933473/144‑1 (MQ=255)
                                                 gCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGCTCATGCCTACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTaa                                                                                              <  1:647332/149‑1 (MQ=255)
                                                 gCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTaa                                                                                              <  2:190932/149‑1 (MQ=255)
                                                 gCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTCAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTaa                                                                                              <  1:598189/149‑1 (MQ=255)
                                                       tgCCGATGGGTGGGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATg                                                                                        <  2:755504/149‑1 (MQ=255)
                                                       tgCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATg                                                                                        >  1:767086/1‑149 (MQ=255)
                                                        gCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATgg                                                                                       <  2:39030/149‑1 (MQ=255)
                                                          cGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCGTAAGGAATGGTc                                                                                     >  2:504139/1‑149 (MQ=255)
                                                           gatggatgCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGTATGGTCa                                                                                    >  1:232435/1‑149 (MQ=255)
                                                            atggatgCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCAt                                                                                   <  2:513439/149‑1 (MQ=255)
                                                                  gCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGaa                                                                             <  2:527709/149‑1 (MQ=255)
                                                                  gCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGaa                                                                             <  1:640161/149‑1 (MQ=255)
                                                                    gggTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGCATGGTCATGCCGAAGc                                                                           >  1:116776/1‑149 (MQ=255)
                                                                       tctgcgtctgcTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATa                                                                        >  2:741858/1‑149 (MQ=255)
                                                                               tgctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTATTGAg                                                                >  2:641211/1‑149 (MQ=255)
                                                                                       ggCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCGTGAAGACACATAATCTTGAGCAACTggc                                                        >  1:820489/1‑149 (MQ=255)
                                                                                        gCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGTGATCATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCCCGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTggcg                                                       <  2:461409/149‑1 (MQ=255)
                                                                                           ttGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGgcggcg                                                    >  2:717483/1‑149 (MQ=255)
                                                                                               cGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAAAATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGt                                                                                      >  1:787764/1‑111 (MQ=255)
                                                                                               cGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAAAATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGt                                                                                      >  1:329419/1‑111 (MQ=255)
                                                                                               cGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAAAATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGt                                                                                      <  2:329419/111‑1 (MQ=255)
                                                                                               cGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAAAATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGt                                                                                      <  2:787764/111‑1 (MQ=255)
                                                                                                  ttCTGAACTGCATATCTTGCGAGAGCAGTTGCACTGGGGGGTGTAGGACGCCTCCGACGCTTTCAATCAACTGGCGTTGAAGTGCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTAc                                             <  1:914332/149‑1 (MQ=255)
                                                                                                  ttCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCCAGACCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTAc                                             >  1:368302/1‑149 (MQ=255)
                                                                                                      gaacctcATTTCTTCCGCGACTTTGTTCACTGGGCGCAGTGGGAACATGTCGAGGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATCCTCTTGAGCAACTGGCGGCGGATTTAcgccg                                         <  1:246557/143‑1 (MQ=255)
                                                                                                      gAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAGGCATACTCTTGAGCAACTGGCGGCGGATTTAagccg                                         >  1:399479/1‑149 (MQ=255)
                                                                                                       aaCTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTAcgccgc                                        >  2:170300/1‑149 (MQ=255)
                                                                                                       aaCTGCATATCTGCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGTACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGGGCAACTGGCGGCGGATTTAcgccgc                                        <  1:74229/149‑1 (MQ=255)
                                                                                                         cTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTAcgccgcgc                                      <  2:861998/149‑1 (MQ=255)
                                                                                                           gCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGAGGCTTTGAATACACTGGTGCTGAATTTCCTCGCAAGCCCTTAAGGAAGGGTCATGACGAAGCATACTCTTTAGCAACTGGCGGCGGATTTAcgccgcgccg                                    >  1:135156/1‑149 (MQ=255)
                                                                                                                 gtgccGCGATGGTGGGCAGTGGGCGCAGTGGCAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGc                              <  2:916642/146‑1 (MQ=255)
                                                                                                                 cTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGc                              >  1:380746/1‑149 (MQ=255)
                                                                                                                  ttCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTAGCGGCGGATTTACGCCGCGCCGCAGAGCa                             >  2:51169/1‑149 (MQ=255)
                                                                                                                  ttCCGCGACTGTGCTCACTGGGCGCAGTGGGAGCAGGCCGACGCGTTCAATCTACGGGTGCCGAATTGCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCa                             <  1:894705/149‑1 (MQ=255)
                                                                                                                     cgcgACTGTGGTCACTGGGCGAAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACGCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAggg                          <  1:722414/149‑1 (MQ=255)
                                                                                                                      gcgACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGc                         >  1:878732/1‑149 (MQ=255)
                                                                                                                        gACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGCTTTACGCCGCGCCGCAGAGCAGGGCGa                       >  1:847937/1‑149 (MQ=255)
                                                                                                                         aCTGTCGTCACTGGGCGCAGTGGTAACAGGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGaa                      <  1:343557/149‑1 (MQ=255)
                                                                                                                           cTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTTATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAg                     <  2:231644/149‑1 (MQ=255)
                                                                                                                              tgGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGa                  >  1:101301/1‑149 (MQ=255)
                                                                                                                                gTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTACTCGCACGCACTTAAGGAATGGTCATGACGAAGCATACTCTTGGGCAACTGACGGCGGATTTACGACGCGCCGCAGACACGCGCCAAGCGAtt                >  2:233672/1‑149 (MQ=255)
                                                                                                                                 tCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTg               >  1:518844/1‑149 (MQ=255)
                                                                                                                                   aCTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCa             >  1:376500/1‑149 (MQ=255)
                                                                                                                                            aGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTgcgc    <  2:28557/149‑1 (MQ=255)
                                                                                                                                              tGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGa  >  2:374452/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGA  >  NC_000913_3_pae_tpiA/371908‑372197

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