Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 913,000 C→T G272D (GGC→GAC)  hcp ← hybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA913,0000CT98.6% 247.7 / NA 70G272D (GGC→GAC) hcphybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/41);  minor base A (0/1);  total (28/42)
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

TCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGTATTTACC  >  NC_000913_3_pae_tpiA/912856‑913145
                                                                                                                                                |                                                                                                                                                 
tCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTgg                                                                                                                                              <  3:351546/150‑1 (MQ=255)
tCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTgg                                                                                                                                              <  2:297303/150‑1 (MQ=255)
tCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTgg                                                                                                                                              >  4:144717/1‑150 (MQ=255)
gcGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTgg                                                                                                                                              <  2:176429/149‑1 (MQ=255)
        gCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGAc                                                                                                                                      <  1:327134/150‑1 (MQ=255)
          aGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACAt                                                                                                                                    <  1:578523/150‑1 (MQ=255)
                   tcatcaCGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCg                                                                                                                           >  2:187483/1‑150 (MQ=255)
                    catcaCGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGt                                                                                                                          <  1:304466/150‑1 (MQ=255)
                    catcaCGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCAAAGTGGGTGTCGACATTAACGCCCGt                                                                                                                          <  4:238542/150‑1 (MQ=255)
                         cGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCtt                                                                                                                     <  3:63863/150‑1 (MQ=255)
                             ggggCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGt                                                                                                                 <  3:307745/150‑1 (MQ=255)
                             ggggCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGt                                                                                                                 <  2:703684/150‑1 (MQ=255)
                                   aGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCg                                                                                                                   <  3:210918/142‑1 (MQ=255)
                                   aGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCg                                                                                                                   >  4:210918/1‑142 (MQ=255)
                                   aGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTc                                                                                                           >  2:265021/1‑150 (MQ=255)
                                      gAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCcag                                                                                                        <  1:475264/150‑1 (MQ=255)
                                           gAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGc                                                                                                                        >  4:439478/1‑129 (MQ=255)
                                           gAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGc                                                                                                                        <  3:439478/129‑1 (MQ=255)
                                                   tCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAgaga                                                                                           >  1:585941/1‑150 (MQ=255)
                                                        ttGCTGATTCTGCCAGCCGCTGCCGTAGTTACCAACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCttt                                                                                      <  1:724500/150‑1 (MQ=255)
                                                         tGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTg                                                                                     <  2:270784/150‑1 (MQ=255)
                                                                    caaGCCGCTGTCGTAGATACCGACCAGATGCATGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGAcc                                                                          >  2:44096/3‑150 (MQ=255)
                                                                     cAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCg                                                                         >  4:167154/1‑150 (MQ=255)
                                                                      aGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCgg                                                                        >  4:278750/1‑150 (MQ=255)
                                                                         cGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCTTGGGTGTAGACATTAACGCGCGTTCCTTCGGTCTGTTCCAGCAGGTTGTAGAGGTCTTTGAGATCGTGACCGGaaa                                                                     >  1:705294/1‑150 (MQ=255)
                                                                            tGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTGGAGATCTTTGAGATCGTGACCGGAAATCa                                                                  <  4:310574/150‑1 (MQ=255)
                                                                              ccGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTag                                                                                              <  1:46825/120‑1 (MQ=255)
                                                                              ccGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTag                                                                                              >  2:46825/1‑120 (MQ=255)
                                                                              ccGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGa                                                                <  4:357303/150‑1 (MQ=255)
                                                                              ccGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGa                                                                >  2:528222/1‑150 (MQ=255)
                                                                               cGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGaa                                                               >  3:63563/1‑150 (MQ=255)
                                                                                gTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTgg                                                                                                                                              >  1:419855/1‑70 (MQ=255)
                                                                                gTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTgg                                                                                                                                              <  2:419855/70‑1 (MQ=255)
                                                                                 tAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCt                                                                                                               >  1:554974/1‑100 (MQ=255)
                                                                                 tAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCt                                                                                                               <  2:554974/100‑1 (MQ=255)
                                                                                 tAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATg                                                             <  3:514/150‑1 (MQ=255)
                                                                                  aGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGc                                                            <  1:37031/150‑1 (MQ=255)
                                                                                   gTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCa                                                           <  1:265021/150‑1 (MQ=255)
                                                                                      aCCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCAttt                                                        >  4:105336/1‑150 (MQ=255)
                                                                                       ccGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTgg                                                                                                                                              <  4:432244/63‑1 (MQ=255)
                                                                                       ccGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTgg                                                                                                                                              >  3:432244/1‑63 (MQ=255)
                                                                                       ccGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTa                                                       >  1:291209/1‑150 (MQ=255)
                                                                                         gACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTAcc                                                     >  2:197808/1‑150 (MQ=255)
                                                                                                  cTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTc                                              >  3:377669/1‑148 (MQ=255)
                                                                                                  cTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTc                                              <  4:377669/148‑1 (MQ=255)
                                                                                                     gAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCATCGCttt                                         <  3:278750/150‑1 (MQ=255)
                                                                                                       aTTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGa                                       <  2:291209/150‑1 (MQ=255)
                                                                                                        tttACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGAt                                                                                 <  1:393108/107‑1 (MQ=255)
                                                                                                        tttACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGAt                                                                                 >  2:393108/1‑107 (MQ=255)
                                                                                                        tttACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGAc                                      <  2:585941/150‑1 (MQ=255)
                                                                                                             gCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCttt                                                                                      <  3:381551/97‑1 (MQ=255)
                                                                                                             gCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCttt                                                                                      >  4:381551/1‑97 (MQ=255)
                                                                                                             gCAGCTCCGGGTAGCCATGCGCAGGCAGCATTGCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGa                                 <  4:20860/150‑1 (MQ=255)
                                                                                                              cAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGAc                                <  1:252212/150‑1 (MQ=255)
                                                                                                                    cGGGTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGGTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGGCCTGGGt                          >  3:263964/1‑150 (MQ=255)
                                                                                                                      ggTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAgg                                                                                                    <  2:312069/74‑1 (MQ=255)
                                                                                                                      ggTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAgg                                                                                                    >  1:312069/1‑74 (MQ=255)
                                                                                                                      ggTAGCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTg                        <  1:528222/150‑1 (MQ=255)
                                                                                                                          gCCATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGt                    <  1:110251/150‑1 (MQ=255)
                                                                                                                           ccATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAAt                                                                    >  1:80634/1‑101 (MQ=255)
                                                                                                                           ccATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAAt                                                                    <  2:80634/101‑1 (MQ=255)
                                                                                                                            cATGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCg                  >  1:834485/1‑150 (MQ=255)
                                                                                                                             aTGCGCAGGCAGCATTTCGTCGTGGGTGTAGACATAAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCgg                 <  2:220940/150‑1 (MQ=255)
                                                                                                                                cgcAGGCAGCATTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTg              <  3:144717/150‑1 (MQ=255)
                                                                                                                                     gcagcaTTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGGCCTGGGTTGGCGTCGGGTGACCGt         >  1:730985/1‑150 (MQ=255)
                                                                                                                                     gcagcaTTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGt         <  4:510484/150‑1 (MQ=255)
                                                                                                                                         caTTTCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGTAttt     >  1:187344/1‑150 (MQ=255)
                                                                                                                                          aTTTCGTCGTGGGTGTAGACATTACCGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGACATCTTCGAGATCGTGACCGGAAATCAGAATGCTTTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGTATTTa    <  2:586188/150‑1 (MQ=255)
                                                                                                                                            ttCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGTATTTAcc  <  4:328365/150‑1 (MQ=255)
                                                                                                                                            ttCGTCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGTATTTAcc  <  4:328351/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TCGATGATGCAGTTCGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGTATTTACC  >  NC_000913_3_pae_tpiA/912856‑913145

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: