Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,609,458 T→C R324R (CGT→CGC)  hyfG → hydrogenase 4, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,609,4580TC99.0% 371.6 / ‑4.3 100R324R (CGT→CGC) hyfGhydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (52/47);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01

GATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACA  >  NC_000913/2609311‑2609604
                                                                                                                                                   |                                                                                                                                                  
gATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCaa                                                                                                                                                  >  2:1015327/1‑150 (MQ=255)
gATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCaa                                                                                                                                                  >  1:1196226/1‑150 (MQ=255)
gATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCaa                                                                                                                                                  >  2:761639/1‑150 (MQ=255)
 aTTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCaaaa                                                                                                                                                <  2:289644/151‑1 (MQ=255)
  ttCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCaaaa                                                                                                                                                <  1:316124/150‑1 (MQ=255)
   tCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAAcc                                                                                                                                              >  1:85336/1‑151 (MQ=255)
     tgctgGAAGTCGCACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCt                                                                                                                                             >  1:527355/1‑150 (MQ=255)
      gctgGAAGTCGAACGGCTGCACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTAc                                                                                                                                           >  1:1224862/1‑151 (MQ=255)
      gctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCt                                                                                                                                             <  1:119900/149‑1 (MQ=255)
      gctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTAc                                                                                                                                           >  1:469937/1‑151 (MQ=255)
         ggAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGt                                                                                                                                        >  1:744053/1‑151 (MQ=255)
           aaGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCt                                                                                                                                      >  1:797343/1‑151 (MQ=255)
             gTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGa                                                                                                                                    <  1:969411/151‑1 (MQ=255)
               cGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAAt                                                                                                                                  >  1:1295819/1‑151 (MQ=255)
                gAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATGTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGCCTGAGTc                                                                                                                                 >  1:1111685/1‑151 (MQ=255)
                 aaCGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCt                                                                                                                                <  1:1339277/151‑1 (MQ=255)
                   cGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTg                                                                                                                               >  2:453637/1‑150 (MQ=255)
                     gCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGAtt                                                                                                                            <  1:1592367/151‑1 (MQ=255)
                          acaGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATtggtgg                                                                                                                       <  2:1351078/151‑1 (MQ=255)
                          acaGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATtggtgg                                                                                                                       >  1:25089/1‑151 (MQ=255)
                           caGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTggtggt                                                                                                                      <  1:1400336/151‑1 (MQ=255)
                               cATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTc                                                                                                                  >  1:32811/1‑151 (MQ=255)
                                aTTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcg                                                                                                                 <  2:209989/151‑1 (MQ=255)
                                aTTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTca                                                                                                                 <  1:1300120/151‑2 (MQ=255)
                                 tttGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcg                                                                                                                 >  1:1466744/1‑150 (MQ=255)
                                  ttGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgcc                                                                                                               >  2:758445/1‑151 (MQ=255)
                                    gCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccgc                                                                                                             >  1:657510/1‑151 (MQ=255)
                                     cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCcgcg                                                                                                            >  1:378855/1‑151 (MQ=255)
                                      ttAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCcgcg                                                                                                            >  2:115006/1‑150 (MQ=255)
                                       tAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCcgcg                                                                                                            <  1:762055/149‑1 (MQ=255)
                                       tAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCcgcg                                                                                                            <  2:1688713/149‑1 (MQ=255)
                                        aaCCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGat                                                                                                          <  1:1536061/150‑1 (MQ=255)
                                        aaCCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGat                                                                                                          <  2:1154547/150‑1 (MQ=255)
                                         aCCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGatat                                                                                                        <  1:1162786/151‑1 (MQ=255)
                                          ccTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATAtt                                                                                                       <  2:7358/151‑1 (MQ=255)
                                            ttGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATtct                                                                                                     >  2:1219222/1‑151 (MQ=255)
                                             tGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATtc                                                                                                      <  1:933685/149‑1 (MQ=255)
                                                ccTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCaaa                                                                                                 <  2:237753/151‑1 (MQ=255)
                                                   ctcCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAag                                                                                                <  1:231089/149‑1 (MQ=255)
                                                   ctcCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAag                                                                                                <  2:378882/149‑1 (MQ=255)
                                                        gCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCaa                                                                                           <  1:1589139/149‑1 (MQ=255)
                                                        gCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAAc                                                                                          <  2:1224957/150‑1 (MQ=255)
                                                         ccATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGt                                                                                        >  1:107583/1‑151 (MQ=255)
                                                         ccATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGt                                                                                        >  1:1337491/1‑151 (MQ=255)
                                                         ccATCTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACAGTCTGAATCTGAGTGGTGGTGTTCGGCGCGATATTCGCaaaaa                                                                                               <  1:912276/144‑3 (MQ=255)
                                                          cATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTc                                                                                       >  1:678002/1‑151 (MQ=255)
                                                           aTTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCt                                                                                      <  1:1492771/151‑1 (MQ=255)
                                                            tttCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTg                                                                                     >  2:387730/1‑151 (MQ=255)
                                                               cGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCaa                                                                                  <  2:788586/151‑1 (MQ=255)
                                                               cGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCaa                                                                                  <  1:467537/151‑1 (MQ=255)
                                                               cGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCa                                                                                   >  1:707463/1‑150 (MQ=255)
                                                                 ttGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAAc                                                                                >  1:1297158/1‑151 (MQ=255)
                                                                 ttGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAAc                                                                                >  2:781697/1‑151 (MQ=255)
                                                                   ggTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGc                                                                              >  2:257680/1‑151 (MQ=255)
                                                                    gTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGc                                                                              >  1:1297736/1‑150 (MQ=255)
                                                                      tttGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTg                                                                            >  1:1296188/1‑150 (MQ=255)
                                                                      tttGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTg                                                                            <  2:1187439/150‑1 (MQ=255)
                                                                       ttGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGa                                                                           <  2:511920/150‑1 (MQ=255)
                                                                            aCCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTgg                                                                     <  1:1044608/151‑1 (MQ=255)
                                                                               ggCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTgc                                                                  >  1:789988/1‑151 (MQ=255)
                                                                                  tttATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGcgcg                                                               >  2:708973/1‑151 (MQ=255)
                                                                                     aTGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGcgcg                                                               >  2:1188860/1‑148 (MQ=255)
                                                                                       gCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATg                                                          <  1:1432395/151‑1 (MQ=255)
                                                                                         aaTTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcg                                                        >  1:336981/1‑151 (MQ=255)
                                                                                          aTTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcgc                                                       >  2:984631/1‑151 (MQ=255)
                                                                                            ttttCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcgcgc                                                     >  2:553761/1‑151 (MQ=255)
                                                                                               tCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCg                                                   >  2:1569594/1‑150 (MQ=255)
                                                                                                ccGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGAc                                                 <  1:853396/151‑1 (MQ=255)
                                                                                                 cgcgTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACg                                                >  2:454454/1‑151 (MQ=255)
                                                                                                   cgTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCTGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACgtg                                              <  1:1162656/151‑1 (MQ=255)
                                                                                                     tGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTc                                            >  1:968046/1‑151 (MQ=255)
                                                                                                      gCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTc                                            <  2:100825/150‑1 (MQ=255)
                                                                                                      gCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCg                                           <  2:647387/151‑1 (MQ=255)
                                                                                                      gCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCg                                           <  1:1015252/151‑1 (MQ=255)
                                                                                                       cGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCgg                                          <  1:101677/151‑1 (MQ=255)
                                                                                                       cGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCgg                                          >  1:165253/1‑151 (MQ=255)
                                                                                                       cGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCgg                                          >  1:1007836/1‑151 (MQ=255)
                                                                                                        gggAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGa                                         <  1:48037/151‑1 (MQ=255)
                                                                                                         ggAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCgg                                          >  1:1695985/1‑149 (MQ=255)
                                                                                                          gAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAg                                        >  1:1034534/1‑150 (MQ=255)
                                                                                                          gAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGc                                       <  1:684516/151‑1 (MQ=255)
                                                                                                          gAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGc                                       <  1:1272721/151‑1 (MQ=255)
                                                                                                           aaaaGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCt                                      <  1:1388198/151‑1 (MQ=255)
                                                                                                            aaaGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTg                                     <  2:469980/151‑1 (MQ=255)
                                                                                                            aaaGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTg                                     >  1:554075/1‑151 (MQ=255)
                                                                                                              aGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGt                                   >  1:613982/1‑151 (MQ=255)
                                                                                                                 cGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGt                                   >  1:1056292/1‑148 (MQ=255)
                                                                                                                 cGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGt                                   >  1:1124596/1‑148 (MQ=255)
                                                                                                                 cGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTag                                 <  2:1002035/150‑1 (MQ=255)
                                                                                                                   aTGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAgaga                              <  2:1092190/151‑1 (MQ=255)
                                                                                                                      aCGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGAtgc                           >  2:908299/1‑151 (MQ=255)
                                                                                                                           ggCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATgctgct                       <  1:1255093/150‑1 (MQ=255)
                                                                                                                              ggAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCtgcttgct                   <  1:235950/151‑1 (MQ=255)
                                                                                                                              ggAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCtgcttgct                   >  2:1006936/1‑151 (MQ=255)
                                                                                                                                 aTTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACg                >  2:1002669/1‑151 (MQ=255)
                                                                                                                                   tGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGcc              >  1:1494035/1‑151 (MQ=255)
                                                                                                                                      tGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGa            >  2:550517/1‑150 (MQ=255)
                                                                                                                                               cgcgCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaaca  <  2:1441407/151‑1 (MQ=255)
                                                                                                                                               cgcgCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaaca  <  2:797409/151‑1 (MQ=255)
                                                                                                                                               cgcgCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaac   <  2:571635/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACA  >  NC_000913/2609311‑2609604

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