Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,601,700 T→A 6.7% L167* (TTA→TAA)  hyfA → hydrogenase 4, 4Fe‑4S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,601,7000TA 6.7% 203.0 / 5.1 60L167* (TTA→TAA) hyfAhydrogenase 4, 4Fe‑4S subunit
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base A (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

TCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGC  >  NC_000913/2601552‑2601848
                                                                                                                                                    |                                                                                                                                                    
tcactcaAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGtt                                                                                                                                                      >  1:271740/1‑149 (MQ=255)
tcactcaAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGtt                                                                                                                                                      >  1:177733/1‑149 (MQ=255)
  actcaAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTAc                                                                                                                                                    >  2:383107/1‑149 (MQ=255)
  actcaAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGCCGTAAc                                                                                                                                                    >  1:348577/1‑149 (MQ=255)
    tcaAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACgg                                                                                                                                                  >  2:242160/1‑149 (MQ=255)
    tcaAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACgg                                                                                                                                                  <  3:200389/149‑1 (MQ=255)
        aGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGa                                                                                                                                              <  1:100106/149‑1 (MQ=255)
         gACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGCCGTAACGCCTGAt                                                                                                                                             <  2:348577/149‑1 (MQ=255)
            ggCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCAc                                                                                                                                          <  3:173198/149‑1 (MQ=255)
                  gaaAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGa                                                                                                                                    >  2:434871/1‑149 (MQ=255)
                     aaCGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAg                                                                                                                                 >  3:57654/1‑149 (MQ=255)
                        gcgcCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCt                                                                                                                              <  2:193458/149‑1 (MQ=255)
                         cgcCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTg                                                                                                                             <  3:380704/149‑1 (MQ=255)
                             aacacaacaTGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAAc                                                                                                                         <  4:541858/149‑1 (MQ=255)
                                acaacaTGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTc                                                                                                                      >  1:471902/1‑149 (MQ=255)
                                  aacaTGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGCCGTAACGCCTGAtt                                                                                                                                            >  1:308759/1‑124 (MQ=255)
                                  aacaTGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGCCGTAACGCCTGAtt                                                                                                                                            <  2:308759/125‑2 (MQ=255)
                                   acaTGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGa                                                                                                                   >  2:574067/1‑149 (MQ=255)
                                   acaTGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGa                                                                                                                   <  2:106384/149‑1 (MQ=255)
                                     aTGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATg                                                                                                                 >  1:202926/1‑149 (MQ=255)
                                        cTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATgaaa                                                                                                              <  4:460033/149‑1 (MQ=255)
                                         tttGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATgaaag                                                                                                             >  2:473881/1‑149 (MQ=255)
                                                   agggAACCTGGGGTCCAGACCGTCACGGTGTACTGCGACCTGTGGGACTTCTTGCCAGAAGTTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAgcg                                                                                                   >  4:99171/2‑149 (MQ=255)
                                                     ggAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGcc                                                                                                 <  4:57654/149‑1 (MQ=255)
                                                        aCCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTg                                                                                              <  1:473881/149‑1 (MQ=255)
                                                          cTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGcc                                                                                            >  3:371686/1‑149 (MQ=255)
                                                                  cAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTg                                                                                    >  4:198845/1‑149 (MQ=255)
                                                                    gACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGt                                                                                  >  4:226645/1‑149 (MQ=255)
                                                                        gTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGGTCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGc                                                                              <  2:271740/149‑1 (MQ=255)
                                                                         tCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGcc                                                                             >  3:422722/1‑149 (MQ=255)
                                                                          cgcgGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCa                                                                            >  2:110901/1‑149 (MQ=255)
                                                                            cgGTGAAATGCGATCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGACTCCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAAt                                                                          <  3:262268/149‑1 (MQ=255)
                                                                                 aaaTGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCgg                                                                     <  4:543084/149‑1 (MQ=255)
                                                                                  aaTGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCggg                                                                    <  1:110901/149‑1 (MQ=255)
                                                                                  aaTGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCggg                                                                    <  2:39678/149‑1 (MQ=255)
                                                                                    tGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGgga                                                                  <  4:422722/149‑1 (MQ=255)
                                                                                     gCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGggag                                                                 >  1:446972/1‑149 (MQ=255)
                                                                                       gACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAGTCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGggagga                                                               >  3:511774/1‑149 (MQ=255)
                                                                                        aCCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGAt                                                              <  4:585535/149‑1 (MQ=255)
                                                                                        aCCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGAt                                                              <  2:202926/149‑1 (MQ=255)
                                                                                          cTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAGTCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTCGCCAATGGCGGGGAGGATcc                                                            >  2:125150/1‑149 (MQ=255)
                                                                                           tgtgtgATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATccc                                                           <  1:194437/149‑1 (MQ=255)
                                                                                             tgtgATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAgc                                                                                                    >  4:211437/1‑106 (MQ=255)
                                                                                             tgtgATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAgc                                                                                                    <  3:211437/106‑1 (MQ=255)
                                                                                               tgATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCttt                                                       >  1:124915/1‑149 (MQ=255)
                                                                                                 aTTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTcc                                                     <  2:177733/149‑1 (MQ=255)
                                                                                                  tttCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTccc                                                    <  1:424379/149‑1 (MQ=255)
                                                                                                   ttCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCct                                                   >  2:273795/1‑149 (MQ=255)
                                                                                                    tCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCg                                                                      >  2:154737/1‑129 (MQ=255)
                                                                                                    tCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCg                                                                      <  1:154737/129‑1 (MQ=255)
                                                                                                     cTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCctca                                                 >  4:171351/1‑149 (MQ=255)
                                                                                                            gAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAGTCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTCGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGa                                          <  1:125150/149‑1 (MQ=255)
                                                                                                            gAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGa                                          <  3:105856/149‑1 (MQ=255)
                                                                                                            gAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCATTCCCTCACTCAGGa                                          <  4:230310/149‑1 (MQ=255)
                                                                                                             aaGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAg                                         <  2:66250/149‑1 (MQ=255)
                                                                                                                             tCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGc                         >  2:579146/1‑149 (MQ=255)
                                                                                                                               gcgcTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGccc                       >  4:353050/1‑149 (MQ=255)
                                                                                                                                       ccGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAAGAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAatta               >  3:338874/1‑149 (MQ=255)
                                                                                                                                       ccGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAatta               >  1:210109/1‑149 (MQ=255)
                                                                                                                                        cGAATCAGGCGTTACGGCCGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAattat              <  3:424342/149‑1 (MQ=255)
                                                                                                                                         gAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAattatt             <  2:446972/149‑1 (MQ=255)
                                                                                                                                             cAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAAcc         >  3:135739/1‑149 (MQ=255)
                                                                                                                                                    tACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGc  <  4:297103/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGC  >  NC_000913/2601552‑2601848

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