Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,843,504 G→A P236L (CCG→CTG)  hycG ← hydrogenase 3 and formate hydrogenase complex, HycG subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,843,5040GA98.5% 262.1 / NA 67P236L (CCG→CTG) hycGhydrogenase 3 and formate hydrogenase complex, HycG subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/33);  minor base T (0/1);  total (33/34)
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

CACCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCGGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCA  >  NC_000913/2843363‑2843651
                                                                                                                                             |                                                                                                                                                   
caccTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGtttt                                                                                                                                            <  1:481014/151‑1 (MQ=255)
 accTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGtttt                                                                                                                                            >  2:133412/1‑150 (MQ=255)
    tgctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGc                                                                                                                                        >  1:644755/1‑151 (MQ=255)
     gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGc                                                                                                                                        <  1:381556/150‑1 (MQ=255)
     gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCt                                                                                                                                       <  1:834897/151‑1 (MQ=255)
     gctgcGCCTCGGCGGGCGTAGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCt                                                                                                                                       <  2:869044/151‑1 (MQ=255)
          gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTccag                                                                                                                                  <  1:672073/151‑1 (MQ=255)
              cggcggGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTccagcc                                                                                                                                <  2:535183/149‑1 (MQ=255)
                gcgggcgTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCg                                                                                                                            <  1:575848/151‑1 (MQ=255)
                 cgggcgTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACGCCACCTGTTCACATATCGGATACTCGCCTCTTCAACATCATTATTCAGATGGCTGACAATTTCGTTCAGACGCTGATCGTTTCCCGCTTCCAGATAGCGt                                                                                                                           <  2:688101/151‑1 (MQ=255)
                     cgTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCaa                                                                                                                       <  1:821765/151‑1 (MQ=255)
                        ggCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCt                                                                                                                    <  2:989306/151‑1 (MQ=255)
                          cATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGt                                                                                                                  >  1:124041/1‑151 (MQ=255)
                           aTCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCAGCGGATACGCGCCTCTTCAACAACATGATTCATATGGCTGACAATCTCGTTCAGACGGAGACAGTTTTCCGCTTCTAGCCTGCAGGCAACCTCtt                                                                                                                 >  2:890961/1‑151 (MQ=255)
                             cGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGttc                                                                                                                <  1:76476/150‑1 (MQ=255)
                              gTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGttctt                                                                                                              <  1:295565/151‑1 (MQ=255)
                              gTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGttctt                                                                                                              <  1:421488/151‑1 (MQ=255)
                               ttCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGttcttc                                                                                                             >  2:293957/1‑151 (MQ=255)
                                tctcATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGttcttc                                                                                                             <  2:401731/150‑1 (MQ=255)
                                   catcaATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCgccc                                                                                                         <  1:237192/151‑1 (MQ=255)
                                      caATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCgccctgc                                                                                                      >  1:211404/1‑151 (MQ=255)
                                      caATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCgccctg                                                                                                       <  1:727232/150‑1 (MQ=255)
                                       aaTAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCgccctgcc                                                                                                     >  1:41835/1‑151 (MQ=255)
                                          aaaTTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCgccctgccct                                                                                                   >  1:852516/1‑150 (MQ=255)
                                                aCGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACtgt                                                                                             <  2:533454/150‑1 (MQ=255)
                                                 cGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTgtgt                                                                                           <  1:355786/151‑1 (MQ=255)
                                                 cGGCCCAGTTGACTGAACACCACCTTTTCACTCATCGGATCCGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGGGCAACCTGTTCTTCGCCCTGCCCTAACtgt                                                                                             >  2:1054401/1‑149 (MQ=255)
                                                  ggCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTgtgt                                                                                           >  1:515229/1‑150 (MQ=255)
                                                   gCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTaa                                                                                         <  1:284355/151‑1 (MQ=255)
                                                    cTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTa                                                                                          >  1:363603/1‑149 (MQ=255)
                                                       aGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCACGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGCGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAgg                                                                                       >  2:1050680/1‑149 (MQ=255)
                                                        gTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTaa                                                                                    <  1:926495/151‑1 (MQ=255)
                                                             cTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGCAatcatc                                                                               >  1:662997/1‑151 (MQ=255)
                                                              tGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCg                                                                              >  2:602984/1‑151 (MQ=255)
                                                               gAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAatcatc                                                                               <  2:295503/149‑1 (MQ=255)
                                                               gAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAatcatc                                                                               <  1:1041050/149‑1 (MQ=255)
                                                               gAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCg                                                                              <  2:406229/150‑1 (MQ=255)
                                                                  caccacCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTACGGTAATCATCGGCa                                                                           >  2:567807/1‑150 (MQ=255)
                                                                  caccacCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCaa                                                                          <  2:342442/151‑1 (MQ=255)
                                                                    ccaccTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATc                                                                        >  1:328720/1‑151 (MQ=255)
                                                                       ccTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTg                                                                      <  1:814028/150‑1 (MQ=255)
                                                                          tttCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTgacg                                                                   >  1:160682/1‑150 (MQ=255)
                                                                           ttCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTgacga                                                                  <  1:1053766/150‑1 (MQ=255)
                                                                            tCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTgacgac                                                                 <  2:629323/150‑1 (MQ=255)
                                                                             cACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCg                                                               >  1:579083/1‑151 (MQ=255)
                                                                                 cATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAAc                                                           >  1:1033837/1‑151 (MQ=255)
                                                                                    cGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGa                                                         <  2:1044345/150‑1 (MQ=255)
                                                                                    cGGAGACGCGCCTCTTCAACAACATGATTCAGATGGCTGGCAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGa                                                         <  2:170745/150‑1 (MQ=255)
                                                                                     ggATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGAt                                                        >  1:89100/1‑150 (MQ=255)
                                                                                           gcgcCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGc                                                 <  1:1051891/151‑1 (MQ=255)
                                                                                            cgcCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGcc                                                <  2:592112/151‑1 (MQ=255)
                                                                                              ccTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGcc                                                >  2:446743/1‑149 (MQ=255)
                                                                                                 ctTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCg                                           >  2:619951/1‑151 (MQ=255)
                                                                                                 ctTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCg                                           >  2:619925/1‑151 (MQ=255)
                                                                                                    caacaaCATGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACg                                        >  1:164750/1‑151 (MQ=255)
                                                                                                       caacaTGATTCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGc                                     >  1:578769/1‑151 (MQ=255)
                                                                                                                tCAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATc                            >  2:266703/1‑151 (MQ=255)
                                                                                                                 cAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATcc                           >  1:929910/1‑151 (MQ=255)
                                                                                                                 cAGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATc                            <  2:260223/150‑1 (MQ=255)
                                                                                                                  aGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATcc                           >  2:710350/1‑150 (MQ=255)
                                                                                                                  aGATGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCa                          <  1:801461/151‑1 (MQ=255)
                                                                                                                    aTGGCTGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACt                        >  2:1009023/1‑151 (MQ=255)
                                                                                                                      ggCTGACAATCTCATTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACtt                       <  2:573619/150‑1 (MQ=255)
                                                                                                                         tGACAATCTCGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCAc                   >  1:943330/1‑151 (MQ=255)
                                                                                                                                 tcGTTCAGACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCg           >  2:648587/1‑151 (MQ=255)
                                                                                                                                        gACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCAc    >  2:258366/1‑151 (MQ=255)
                                                                                                                                        gACGCAGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCAc    >  2:767836/1‑151 (MQ=255)
                                                                                                                                          cgcaGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCa  <  2:1929/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CACCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCCTCTTCAACAACATGATTCAGATGGCTGACAATCTCGTTCAGACGCGGATCGTTTTCCGCTTCCAGCCAGCGTGCAACCTGTTCTTCGCCCTGCCCTAACTGTGTAAGGTAATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCA  >  NC_000913/2843363‑2843651

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