Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,846,307 G→A G55G (GGC→GGT hycE ← hydrogenase 3, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,846,3070GA98.9% 357.5 / ‑4.1 94G55G (GGC→GGThycEhydrogenase 3, large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (43/50);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

CACGGACGGATATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCGCGGCGAGATAATGTT  >  NC_000913/2846160‑2846449
                                                                                                                                                   |                                                                                                                                              
cACGGACGGATATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACagccaacc                                                                                                                                              <  1:474900/150‑1 (MQ=255)
cACGGACGGATATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACagccaacc                                                                                                                                              >  1:1285862/1‑150 (MQ=255)
cACGGACGGATATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACagccaacc                                                                                                                                              >  2:395149/1‑150 (MQ=255)
    gacggaTATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCt                                                                                                                                         <  2:563783/151‑1 (MQ=255)
     acggaTATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTg                                                                                                                                        <  1:619870/151‑1 (MQ=255)
      cggaTATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGt                                                                                                                                       <  2:986495/151‑1 (MQ=255)
      cggaTATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGt                                                                                                                                       <  2:598978/151‑1 (MQ=255)
      cggaTATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGt                                                                                                                                       >  1:680938/1‑151 (MQ=255)
           atTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTgt                                                                                                                                   >  1:417164/1‑150 (MQ=255)
            ttCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTgt                                                                                                                                   >  2:1468088/1‑149 (MQ=255)
              ccGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTgtagta                                                                                                                               <  2:407239/151‑1 (MQ=255)
              ccGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTgtagta                                                                                                                               <  1:1494273/151‑1 (MQ=255)
                ggTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTaaa                                                                                                                             <  1:1125314/151‑1 (MQ=255)
                     gttgGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGaaa                                                                                                                         <  2:680996/150‑1 (MQ=255)
                      ttgGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACt                                                                                                                       >  1:785406/1‑151 (MQ=255)
                        ggCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTcc                                                                                                                     >  2:16907/1‑151 (MQ=255)
                         gCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAACGAAACTcca                                                                                                                    <  2:582728/151‑1 (MQ=255)
                         gCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTcc                                                                                                                     >  2:590225/1‑150 (MQ=255)
                         gCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTc                                                                                                                      <  1:1049991/149‑1 (MQ=255)
                           gTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTccacc                                                                                                                  >  2:835747/1‑151 (MQ=255)
                             cAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCCCCCTGTTTGTAGTAAAGAAACTccacca                                                                                                                 <  1:526049/150‑1 (MQ=255)
                             cAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTccacca                                                                                                                 <  2:1590262/150‑1 (MQ=255)
                                  tCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTcc                                                                                                            <  1:437387/150‑1 (MQ=255)
                                     aCGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCa                                                                                                        <  2:1534003/151‑1 (MQ=255)
                                     aCGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGc                                                                                                         >  1:1015501/1‑150 (MQ=255)
                                       gcgAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAgg                                                                                                      >  2:1176440/1‑151 (MQ=255)
                                        cgAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGt                                                                                                     >  1:704998/1‑151 (MQ=255)
                                         gAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTa                                                                                                    <  1:628223/151‑1 (MQ=255)
                                            ccGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGtt                                                                                                 >  1:349478/1‑151 (MQ=255)
                                             cGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTc                                                                                                >  2:585123/1‑151 (MQ=255)
                                              gTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCa                                                                                               <  1:1044578/151‑1 (MQ=255)
                                                   ccAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCtt                                                                                           <  1:1600371/150‑1 (MQ=255)
                                                    cAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCttt                                                                                          >  2:1658897/1‑150 (MQ=255)
                                                    cAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCttt                                                                                          >  2:1234313/1‑150 (MQ=255)
                                                      acacTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACg                                                                                       >  2:1079577/1‑151 (MQ=255)
                                                      acacTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACg                                                                                       <  1:704696/151‑1 (MQ=255)
                                                        acTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTTGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGt                                                                                     <  1:1511838/151‑1 (MQ=255)
                                                           tAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGa                                                                                   <  1:346489/150‑1 (MQ=255)
                                                            aGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACa                                                                                 <  2:719474/151‑1 (MQ=255)
                                                             gTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGAc                                                                                  >  1:1408380/1‑149 (MQ=255)
                                                              tGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGt                                                                               >  2:759985/1‑151 (MQ=255)
                                                              tGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGt                                                                               >  1:1455582/1‑151 (MQ=255)
                                                               gCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGt                                                                               <  2:161426/150‑1 (MQ=255)
                                                               gCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGGAGTTCACCTTTACGGTGACAGTc                                                                              <  1:263672/151‑1 (MQ=255)
                                                                 ccTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTc                                                                              >  1:1360705/1‑149 (MQ=255)
                                                                 ccTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAg                                                                            <  2:925959/151‑1 (MQ=255)
                                                                   ttCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCt                                                                          >  1:36975/1‑151 (MQ=255)
                                                                   ttCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCt                                                                          >  1:1084929/1‑151 (MQ=255)
                                                                       ccATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGAt                                                                       <  1:1457302/150‑1 (MQ=255)
                                                                       ccATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGAt                                                                       <  2:1000348/150‑1 (MQ=255)
                                                                           agACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTg                                                                  >  1:142999/2‑151 (MQ=255)
                                                                             aCAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGt                                                                >  2:1211323/1‑151 (MQ=255)
                                                                               aGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCt                                                              >  1:67363/1‑151 (MQ=255)
                                                                               aGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCt                                                              >  2:774388/1‑151 (MQ=255)
                                                                                  aCGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGcc                                                           >  1:806947/1‑151 (MQ=255)
                                                                                    gtagtaAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGcc                                                           >  2:188776/1‑149 (MQ=255)
                                                                                      agtaAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGc                                                       >  1:1619204/1‑151 (MQ=255)
                                                                                      agtaAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGc                                                       <  2:649892/151‑1 (MQ=255)
                                                                                         aaaCGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTg                                                    >  2:408816/1‑151 (MQ=255)
                                                                                         aaaCGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTg                                                    <  2:618951/151‑1 (MQ=255)
                                                                                           aCGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGt                                                  >  1:373919/1‑151 (MQ=255)
                                                                                            cGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGt                                                  <  2:1179637/150‑1 (MQ=255)
                                                                                              gCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTcct                                               <  1:1551448/151‑2 (MQ=255)
                                                                                               cGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTcc                                                >  1:167896/1‑149 (MQ=255)
                                                                                                   aTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCacg                                          >  2:940450/1‑151 (MQ=255)
                                                                                                   aTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCacg                                          <  1:1473999/151‑1 (MQ=255)
                                                                                                     gACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCacga                                         <  2:484156/150‑1 (MQ=255)
                                                                                                      aCCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCacgacg                                       <  2:1193004/151‑1 (MQ=255)
                                                                                                      aCCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCacgacg                                       <  2:835774/151‑1 (MQ=255)
                                                                                                        cATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGcc                                     >  1:1684516/1‑151 (MQ=255)
                                                                                                        cATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGcc                                     >  2:369594/1‑151 (MQ=255)
                                                                                                           tCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCg                                   <  1:70788/150‑1 (MQ=255)
                                                                                                              gattGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGa                                 >  1:809886/3‑149 (MQ=255)
                                                                                                              gTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAACCGCCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGaaa                               <  1:1528900/151‑1 (MQ=255)
                                                                                                                   cGTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCt                          >  1:1419448/1‑151 (MQ=255)
                                                                                                                    gTTCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTc                         >  2:1116487/1‑151 (MQ=255)
                                                                                                                     ttCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTc                         >  2:903371/1‑150 (MQ=255)
                                                                                                                     ttCGTCGTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTc                         >  2:114832/1‑150 (MQ=255)
                                                                                                                      tcgtcgTTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCTCCTCCGGCAGGTAGTTCTCCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCa                        <  1:1622829/150‑1 (MQ=255)
                                                                                                                       cgtcgtTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCAt                       <  1:1346508/150‑1 (MQ=255)
                                                                                                                          cgtTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCGGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCa                    <  1:395121/150‑1 (MQ=255)
                                                                                                                          cgtTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCa                    <  1:1433632/150‑1 (MQ=255)
                                                                                                                          cgtTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCa                    <  1:481284/150‑1 (MQ=255)
                                                                                                                          cgtTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCa                    <  1:1180491/150‑1 (MQ=255)
                                                                                                                           gtTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAgc                  <  2:543909/151‑1 (MQ=255)
                                                                                                                           gtTACCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAgc                  <  1:95285/151‑1 (MQ=255)
                                                                                                                              aCCAAACAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGcgcg               <  2:519145/151‑1 (MQ=255)
                                                                                                                                 aaaCAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCgcggcg            >  1:768689/1‑151 (MQ=255)
                                                                                                                                  aaCAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCgcggcg            >  1:182340/1‑150 (MQ=255)
                                                                                                                                  aaCAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCgcggcg            >  1:283904/1‑150 (MQ=255)
                                                                                                                                   aCAGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCGCGGCgag          <  2:1285955/151‑1 (MQ=255)
                                                                                                                                     aGCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCGCGGCGAGAt        <  2:1283755/151‑1 (MQ=255)
                                                                                                                                      gCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCGCGGCGAGAt        >  2:575941/1‑150 (MQ=255)
                                                                                                                                      gCACCGACAGCCAACCACCCTGTTTGTAGTAAAGAAACGCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCGCGGCGAGATa       <  2:758010/151‑1 (MQ=255)
                                                                                                                                         ccGACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCGCGGCGAGATAAt     <  1:642470/150‑1 (MQ=255)
                                                                                                                                           gACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCGCGGCGAGATAATGtt  <  1:646515/151‑1 (MQ=255)
                                                                                                                                           gACAGCCAACCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCGCGGCGAGATAATGtt  <  1:452524/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CACGGACGGATATTCCGGTTTGTTGGCGTCAACTTCGACGCGAACCGTAATCCAACACTTAGTGCCCTTCTCCATCGACAGCACGTAGTAAACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGGTAGTTCACCTTTACGGTGACAGTCAGCTGATCTTTGGTCTGCCAGGCGTGGTCCAGCACGACGCCCGGAAATGCCTCATTCAGCGCGGCGAGATAATGTT  >  NC_000913/2846160‑2846449

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