Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,256,450 C→T 9.1% V68I (GTT→ATT)  cyuA ← putative L‑cysteine desulfidase CyuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,256,4500CT 9.1% 192.6 / 7.4 66V68I (GTT→ATT) cyuAputative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (4/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

GGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTTCCTCCT  >  NC_000913/3256301‑3256599
                                                                                                                                                     |                                                                                                                                                     
ggAGACTTTCCCCGCCGCCAGCAATGAATTATCATCGGCAATTTACTGAGCTGTTGCGTATATAAGAAATTACAGCCCGGCGATGAAAATTACAAATAACGCCCCCAGCGCCGCCGCAATTGGCAGCCCCACCATTCCCGTGCCGGGAAc                                                                                                                                                       >  2:243555/1‑150 (MQ=255)
 gagaCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACg                                                                                                                                                      >  2:79733/1‑150 (MQ=255)
  agaCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACgg                                                                                                                                                     <  2:243194/150‑1 (MQ=255)
      tttCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCGTTCCCGGCACGGGAATGGTggg                                                                                                                                                 <  1:1038009/150‑3 (MQ=255)
      tttCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCGTTCCCGGCACGGGAATGGTggg                                                                                                                                                 <  1:387801/150‑3 (MQ=255)
       ttCCCCGCCGCAAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACg                                                                                                                                                >  2:772763/1‑150 (MQ=255)
          cccGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGccc                                                                                                                                             <  2:1188637/150‑1 (MQ=255)
           ccgccgCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCa                                                                                                                                            <  1:1223805/150‑1 (MQ=255)
                gccAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCg                                                                                                                                       >  2:1178443/1‑150 (MQ=255)
                   agcagTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGttc                                                                                                                                    >  2:1000682/1‑150 (MQ=255)
                      agTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGttcttc                                                                                                                                 <  2:403479/150‑1 (MQ=255)
                        tGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcat                                                                                                                               <  1:7303/150‑1 (MQ=255)
                        tGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcat                                                                                                                               <  1:12492/150‑1 (MQ=255)
                         gCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                                                                                              <  2:1281599/150‑1 (MQ=255)
                         gCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                                                                                              <  1:328185/150‑1 (MQ=255)
                              ggCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGCCGCCCAGACCGTTCTTCATCAGAtt                                                                                                                         >  1:1140313/1‑150 (MQ=255)
                                    ggCCATTGCCTGCGCTGTTGCGTCTATCCGCACTTCCAGCCCGGCGTTGGCATTTCCCCCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                                   >  1:21898/1‑150 (MQ=255)
                                    ggCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                                   >  1:21900/1‑150 (MQ=255)
                                    ggCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                                   >  2:768937/1‑150 (MQ=255)
                                      cAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGaaa                                                                                                                 >  2:259271/1‑150 (MQ=255)
                                          tGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAACCCCa                                                                                                             >  1:1321953/1‑150 (MQ=255)
                                           gCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAg                                                                                                            <  1:249267/150‑1 (MQ=255)
                                              tGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCt                                                                                                         <  1:156479/150‑1 (MQ=255)
                                                 gcTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCt                                                                                                      <  1:667519/150‑1 (MQ=255)
                                                   tGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTac                                                                                                    <  1:49708/150‑1 (MQ=255)
                                                     ttGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTacac                                                                                                  <  1:772763/150‑1 (MQ=255)
                                                      tGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCTTTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACg                                                                                                 <  2:1325883/150‑1 (MQ=255)
                                                       gCGTCTTTCGGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGt                                                                                                <  1:852162/150‑1 (MQ=255)
                                                             ttCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAAcc                                                                                          <  1:92514/150‑1 (MQ=255)
                                                                aGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGa                                                                                       <  1:1025644/150‑1 (MQ=255)
                                                                 gCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCACCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGCAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGAc                                                                                      >  2:66291/1‑150 (MQ=255)
                                                                      tCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTcc                                                                                 <  1:1000682/150‑1 (MQ=255)
                                                                             cGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTg                                                                          >  2:223866/1‑150 (MQ=255)
                                                                             cGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTg                                                                          <  1:1291592/150‑1 (MQ=255)
                                                                                        tttCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCACACCGGTCGTCATCAGAGTTGGCCAAACCCCGGCTTCGACACGTTCCACCGGACCTTCCACTTCAGCCGCAGCAccc                                                               >  2:567975/1‑150 (MQ=255)
                                                                                          tCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGGCCTTCCAGTTCTGCCGCAGCAAccgc                                                             >  1:1200586/1‑150 (MQ=255)
                                                                                          tCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAccgc                                                             >  1:617967/1‑150 (MQ=255)
                                                                                            cACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAccgccg                                                           <  2:604304/150‑1 (MQ=255)
                                                                                                    gCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagc                                                   >  2:494988/1‑150 (MQ=255)
                                                                                                     cccccAGCGCCGCCGCAATTGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGTCTCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCATCCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcg                                                  <  1:243555/150‑1 (MQ=255)
                                                                                                     cccccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcg                                                  <  1:259271/150‑1 (MQ=255)
                                                                                                      ccccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGCCCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgc                                                 >  1:1349416/1‑150 (MQ=255)
                                                                                                      ccccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgc                                                 >  2:473548/1‑150 (MQ=255)
                                                                                                       cccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgcc                                                >  1:408698/1‑150 (MQ=255)
                                                                                                          aGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGt                                             <  2:393726/150‑1 (MQ=255)
                                                                                                            cgccgccgcAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGa                                           >  1:1120234/1‑150 (MQ=255)
                                                                                                                gccgcAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGCCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATc                                       <  1:518351/150‑1 (MQ=255)
                                                                                                                       tCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGGCCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAg                                <  2:1200586/150‑1 (MQ=255)
                                                                                                                        cGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGt                               >  1:789720/1‑150 (MQ=255)
                                                                                                                        cGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGt                               >  1:1158502/1‑150 (MQ=255)
                                                                                                                         ggcggcCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTa                              <  2:1282004/150‑1 (MQ=255)
                                                                                                                             gCCCCACCATTCCCGTGCCGGGAACGGTGCCGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATc                          >  1:339919/1‑150 (MQ=255)
                                                                                                                              ccccACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATcc                         >  2:1130705/1‑150 (MQ=255)
                                                                                                                              ccccACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATcc                         >  2:336820/1‑150 (MQ=255)
                                                                                                                               cccACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATccc                        >  1:1229696/1‑150 (MQ=255)
                                                                                                                                 caccgTTCCCGGCACGGGAATGGTGGCGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAg                      >  2:81906/1‑150 (MQ=255)
                                                                                                                                 caccgTTCCCGGCACGGGAATGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAg                      >  1:271641/1‑150 (MQ=255)
                                                                                                                                  accgTTCCCGGCACGGGAATGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAgc                     >  2:1416714/1‑150 (MQ=255)
                                                                                                                                  accgTTCCCGGCACGGGAATGGTGACGCCCAGACCGTTCTTCATCAGAATTGGAGAACACCAGGCTTCTACACGTTCAACCGGCACTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCCGCGCCAGTGAAAGCGGTTCAGTACATCCCAgc                     >  2:654341/1‑150 (MQ=255)
                                                                                                                                  accaTTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAgc                     <  2:21900/150‑1 (MQ=255)
                                                                                                                                  accaTTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTGCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAgc                     <  2:21898/150‑1 (MQ=255)
                                                                                                                                        cccGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACCCGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGt               >  2:267195/1‑150 (MQ=255)
                                                                                                                                           gggCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGtttt            <  1:768313/148‑1 (MQ=255)
                                                                                                                                           gTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGtttt            <  1:66291/150‑1 (MQ=255)
                                                                                                                                                gggAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTtc       >  2:1262071/1‑150 (MQ=255)
                                                                                                                                                 ttaaCGGTTACGCCCAGACCGTTCTTCATCAGATTTGTCGAAACCGAGGATTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTtcc      <  1:718745/148‑1 (MQ=255)
                                                                                                                                                 ggAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTtcc      <  1:1178443/150‑1 (MQ=255)
                                                                                                                                                     cGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCCCCCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTTccgcct  >  2:1018107/1‑150 (MQ=255)
                                                                                                                                                     cGGTGAAGCCCAGACCGTTATTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCCGCGACGGTTTTACTTcctcct  >  2:1023240/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTTCCTCCT  >  NC_000913/3256301‑3256599

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