Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,602,017 (C)7→8 coding (200/2019 nt) hyfB → hydrogenase 4 component B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,602,0101.C90.0% 182.5 / 13.9 60coding (193/2019 nt)hyfBhydrogenase 4 component B
Reads supporting (aligned to +/- strand):  ref base . (1/5);  new base C (28/26);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

TGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAG  >  NC_000913/2601873‑2602138
                                                                                                                                          |                                                                                                                                  
tGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑cc                                                                                                                                <  1:934469/140‑1 (MQ=255)
  gCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑cccc                                                                                                                              >  2:455707/1‑140 (MQ=255)
        cTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCctt                                                                                                                        <  2:942542/140‑1 (MQ=255)
            ttttACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑ccccc                                                                                                                             <  2:955206/131‑1 (MQ=255)
            ttttACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑ccccc                                                                                                                             >  1:955206/1‑131 (MQ=255)
                aCTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAAttt                                                                                                                 <  2:995523/140‑1 (MQ=255)
                aCTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCTTTTGAATTTg                                                                                                                <  1:197576/140‑1 (MQ=255)
                     gTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGc                                                                                                               >  2:524012/1‑137 (MQ=255)
                     gTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGc                                                                                                               <  1:524012/137‑1 (MQ=255)
                           aCAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAc                                                                                                      <  2:59439/140‑1 (MQ=255)
                            cAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAcc                                                                                                     <  1:455707/140‑1 (MQ=255)
                               aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTg                                                                                                  >  1:32309/1‑140 (MQ=255)
                               aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTg                                                                                                  >  2:524632/1‑140 (MQ=255)
                               aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTg                                                                                                  >  1:821969/1‑140 (MQ=255)
                               aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTg                                                                                                  >  1:64125/1‑140 (MQ=255)
                                cTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGc                                                                                                 >  1:189157/1‑140 (MQ=255)
                                 tGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCg                                                                                                <  1:465779/140‑1 (MQ=255)
                                    cTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAAt                                                                                             >  2:434994/1‑140 (MQ=255)
                                     tATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGTCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATg                                                                                            <  1:390614/140‑1 (MQ=255)
                                         aaGCTTTCCGGCATCACATCGGTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATa                                                                                        >  2:768732/1‑140 (MQ=255)
                                         aaGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                                    <  2:22525/128‑1 (MQ=255)
                                         aaGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                                    >  1:22525/1‑128 (MQ=255)
                                             tttCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCct                                                                                    >  1:402474/1‑140 (MQ=255)
                                               tCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                                    >  2:647882/1‑122 (MQ=255)
                                               tCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                                    <  1:647882/122‑1 (MQ=255)
                                                ccGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCtctc                                                                                 <  2:189157/140‑1 (MQ=255)
                                                 cGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCg                                                                                >  1:262194/1‑140 (MQ=255)
                                                  ggCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTAACCCTGCGAATGGATAGCCTCTCGGc                                                                              <  1:536787/140‑1 (MQ=255)
                                                   gCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGc                                                                              <  1:775749/140‑1 (MQ=255)
                                                   gCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGc                                                                              <  2:262194/140‑1 (MQ=255)
                                                        acaTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTa                                                                         <  1:667269/140‑1 (MQ=255)
                                                        acaTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTa                                                                         <  2:821969/140‑1 (MQ=255)
                                                           tCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATgg                                                                      <  2:1030273/140‑1 (MQ=255)
                                                            cGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGt                                                                     >  2:686449/1‑140 (MQ=255)
                                                             gCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTTCGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTg                                                                    <  1:898231/140‑1 (MQ=255)
                                                                   gggCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGtt                                                              <  2:32309/140‑1 (MQ=255)
                                                                         cgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTcc                                                        >  2:90658/1‑140 (MQ=255)
                                                                         cgTGATTGGCATCATCAGCGGAATTACGCAATTACATACTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCTAATGGATGGCCTCTCGGCATTTATGGTGCTGGTTGTCTcc                                                        >  2:960190/1‑140 (MQ=255)
                                                                            gATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGc                                                    <  2:765171/140‑1 (MQ=255)
                                                                                  catcatCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtg                                               <  2:250047/140‑1 (MQ=255)
                                                                                       tcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt                                          >  2:813212/1‑140 (MQ=255)
                                                                                       tcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt                                          <  1:813212/140‑1 (MQ=255)
                                                                                           cGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCCGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGt                                      <  2:617501/140‑1 (MQ=255)
                                                                                            ggAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGtt                                     <  2:402474/140‑1 (MQ=255)
                                                                                            ggAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGtt                                     <  1:768732/140‑1 (MQ=255)
                                                                                                  aCGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGTTAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTctctct                               >  2:1040571/1‑140 (MQ=255)
                                                                                                  aCGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTctctct                               <  2:64125/140‑1 (MQ=255)
                                                                                                  aCGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CGCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTctctct                               <  1:638792/140‑1 (MQ=255)
                                                                                                       aTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCa                          <  1:90658/140‑1 (MQ=255)
                                                                                                        ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCAt                         <  1:523312/139‑1 (MQ=255)
                                                                                                        ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCAt                         >  2:523312/1‑139 (MQ=255)
                                                                                                             tGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACt                    >  2:854434/1‑140 (MQ=255)
                                                                                                                  gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCa                          >  1:1058002/1‑129 (MQ=255)
                                                                                                                  gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCa                          <  2:1058002/129‑1 (MQ=255)
                                                                                                                  gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCa                          >  1:673750/1‑129 (MQ=255)
                                                                                                                  gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTCCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCa                          <  2:673750/129‑1 (MQ=255)
                                                                                                                   tgtAACTTTAGTCGCCCGTTTTGCCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCAGGTGGTGGTTTGTTCTCTCTATTCATTGACTTatat               >  2:468277/1‑140 (MQ=255)
                                                                                                                       aCTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGcgcg          >  1:837067/1‑140 (MQ=255)
                                                                                                                            aGTCGCCCGTTTTGCCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATa      >  2:371332/1‑140 (MQ=255)
                                                                                                                            aGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATAc     >  1:311383/1‑140 (MQ=255)
                                                                                                                            aGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATAc     >  2:527952/1‑140 (MQ=255)
                                                                                                                            aGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATAc     >  2:317568/1‑140 (MQ=255)
                                                                                                                            aGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATAc     >  2:1087235/1‑140 (MQ=255)
                                                                                                                             gTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACg    >  2:124582/1‑140 (MQ=255)
                                                                                                                               cGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAg  >  2:130656/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                  
TGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAG  >  NC_000913/2601873‑2602138

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