Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,601,696 G→C 6.3% A166P (GCG→CCG)  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,6960GC 6.3% 219.4 / 5.7 63A166P (GCG→CCG) hyfAhydrogenase 4, 4Fe‑4S subunit
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base C (2/2);  total (32/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 = 7.71e-01

CTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACC  >  NC_000913/2601551‑2601841
                                                                                                                                                 |                                                                                                                                                 
ctcactcaAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGt                                                                                                                                                >  2:346413/1‑149 (MQ=255)
 tcactcaAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGtt                                                                                                                                               >  1:177733/1‑149 (MQ=255)
 tcactcaAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGtt                                                                                                                                               >  1:271740/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                                                             <  2:39678/149‑1 (MQ=255)
                                                                                   aaTGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCggg                                                             <  1:110901/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                                                                                             <  3:211437/106‑1 (MQ=255)
                                                                                              tgtgATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAgc                                                                                             >  4:211437/1‑106 (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                                                               <  1:154737/129‑1 (MQ=255)
                                                                                                     tCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCg                                                               >  2:154737/1‑129 (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)
                                                                                                                                                 |                                                                                                                                                 
CTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGTGTCCAGACCGTCGCGGTGAAATGCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACC  >  NC_000913/2601551‑2601841

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: ‑