Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,892,363 A→G 15.4% Y291Y (TAT→TAC yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,892,3630AG15.4% 310.4 / 45.3 135Y291Y (TAT→TACyoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base A (56/58);  new base G (11/10);  total (67/68)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 138 total aligned reads displayed.

GTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGG  >  NC_000913/1892266‑1892462
                                                                                                 |                                                                                                   
gTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATa                                                                                                   >  1:3293152/1‑100 (MQ=255)
gTTTCGCCACGTCATAAAAAAGTTCCAGGGTGTCATCGAGCAATAAAAATGCCCGCTGAATCTGCGGACTAGCCAACAGCTCACCCAGGATACCGAGATa                                                                                                   >  1:2121481/1‑100 (MQ=255)
 tttCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATaa                                                                                                  >  1:3001575/1‑100 (MQ=255)
   tCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAAcc                                                                                                <  2:2692281/100‑1 (MQ=255)
    cGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCt                                                                                               >  1:2395174/1‑100 (MQ=255)
     gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTg                                                                                              <  1:1229636/100‑1 (MQ=255)
     gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTg                                                                                              <  2:2076966/100‑1 (MQ=255)
     gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTCGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                                                              >  2:2550320/1‑100 (MQ=255)
       cACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcg                                                                                            <  2:2780542/100‑3 (MQ=255)
       cACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGc                                                                                            <  1:1057972/100‑1 (MQ=255)
        aCGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgt                                                                                           <  1:3491001/100‑4 (MQ=255)
        aCGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgt                                                                                           <  1:2133930/100‑4 (MQ=255)
          gTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgtga                                                                                         <  1:860662/100‑6 (MQ=255)
           tCATAACAAATTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGc                                                                                            <  1:670820/96‑1 (MQ=255)
           tCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgtgag                                                                                        <  1:3186340/100‑7 (MQ=39)
           tCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAc                                                                                        <  1:2008899/100‑1 (MQ=255)
            cATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAcc                                                                                       >  2:3619895/1‑100 (MQ=255)
            cATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAc                                                                                        <  2:2309456/99‑1 (MQ=255)
             aTAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCg                                                                                      <  2:2016399/100‑1 (MQ=255)
              tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCg                                                                                      >  1:3270809/1‑99 (MQ=255)
              tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGa                                                                                     >  1:1351073/1‑100 (MQ=255)
               aaCAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAg                                                                                    >  2:778584/1‑100 (MQ=255)
                aCAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGt                                                                                   >  2:3607088/1‑100 (MQ=255)
                 cAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGtt                                                                                  >  2:2584613/1‑100 (MQ=255)
                   aaGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGc                                                                                <  2:866254/100‑1 (MQ=255)
                    aGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCa                                                                               <  2:1136874/100‑1 (MQ=255)
                    aGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCa                                                                               >  2:233582/1‑100 (MQ=255)
                       tCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGAc                                                                            <  1:3607115/100‑1 (MQ=255)
                       tCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGAc                                                                            <  1:2912527/100‑1 (MQ=255)
                        ctaGGGTGTCATGCTACAGTATAAATGGCCGCTGTGTTTGCGGATTGGTTGCCTGCTCACTCAGGTTACTGCGATAACCTCGCTCACCGAGTTGCAg                                                                              <  2:1752886/95‑1 (MQ=255)
                          aGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGaa                                                                         <  2:815181/100‑1 (MQ=255)
                           gggTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGaaa                                                                        <  2:3198930/100‑1 (MQ=255)
                             gtgtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAAt                                                                      >  1:1821166/1‑100 (MQ=255)
                              tgtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATc                                                                     <  1:2904848/100‑1 (MQ=255)
                               gtCATCGAGCAGTAAGAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATcc                                                                    <  2:3088915/100‑1 (MQ=255)
                                 cATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                                                                  >  1:2073938/1‑100 (MQ=255)
                                 cATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                                                                  >  2:541840/1‑100 (MQ=255)
                                 cATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                                                                  >  2:3000854/1‑100 (MQ=255)
                                  aTCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgc                                                                 >  1:3075488/1‑100 (MQ=255)
                                  aTCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgc                                                                 <  1:1520418/100‑1 (MQ=255)
                                  aTCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgc                                                                 <  1:1814223/100‑1 (MQ=255)
                                   tCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgcg                                                                <  1:419941/100‑1 (MQ=255)
                                   tCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgcg                                                                <  1:3534781/100‑1 (MQ=255)
                                    cGAGCAGTAAAAATGGCCGCTGAATTTGCGTATTGGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                                                                  <  2:3474498/97‑1 (MQ=255)
                                          gTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTg                                                         >  1:860677/1‑100 (MQ=255)
                                          gTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTg                                                         >  1:2888663/1‑100 (MQ=255)
                                          gTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTAACGCAGGTTACCGCGGTAACCTGGATCACCGAGTTGCAGCCGAACATCCTGCGCACTCTg                                                         >  2:1634475/1‑100 (MQ=255)
                                             aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGc                                                      <  2:3122719/100‑1 (MQ=255)
                                             aaaaTGCCCGCTGAACTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGc                                                      <  1:1017202/100‑1 (MQ=255)
                                              aaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCa                                                     <  2:2888640/100‑1 (MQ=255)
                                              aaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCa                                                     <  1:1155256/100‑1 (MQ=255)
                                               aaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCaa                                                    <  2:928018/100‑1 (MQ=255)
                                                aTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAg                                                   <  1:2541725/100‑1 (MQ=255)
                                                  gCCCGCTGAATTTGCGCATTCGCTAACAGCGCACGCAGGATACCGCGATAACCTGGTTCACCGAGTTGCGGACGAAAATCCTGCGCACTTTTGGCAAGAc                                                 <  2:1301013/100‑1 (MQ=255)
                                                    cgggtGAAGTTGCGGATTAGCTGACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGGTGCAGACGAGAATCCTGCGCACTCTGGGCAAGa                                                  <  2:285921/93‑1 (MQ=255)
                                                    ccGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGa                                               >  1:751424/1‑100 (MQ=255)
                                                     cGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGAt                                              >  1:2783223/1‑100 (MQ=255)
                                                       cTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa                                            >  2:170596/1‑100 (MQ=255)
                                                          aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAgcg                                         <  2:3102004/100‑1 (MQ=255)
                                                           aTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAgcgc                                        <  2:637734/100‑1 (MQ=255)
                                                              tGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACt                                     >  2:932992/1‑100 (MQ=255)
                                                              tGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACt                                     <  1:3209918/100‑1 (MQ=255)
                                                              tGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACt                                     >  1:2331014/1‑100 (MQ=255)
                                                                  gATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTg                                 <  2:2182088/100‑1 (MQ=255)
                                                                   aTTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGt                                >  1:980100/1‑100 (MQ=255)
                                                                      aGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAAc                             <  1:960320/100‑1 (MQ=255)
                                                                      aGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAAc                             >  1:466050/1‑100 (MQ=255)
                                                                      aGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAACGAGTTGCAGACGAAAATCCTGCGAAAAATGGGCAAGACGATCAGCGCACTTTTGTAAc                             >  2:2147193/1‑100 (MQ=255)
                                                                       gCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAAct                            >  1:2210977/1‑100 (MQ=255)
                                                                         tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAActgc                          >  1:1591721/1‑100 (MQ=255)
                                                                         tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAActgc                          >  1:1437594/1‑100 (MQ=255)
                                                                         tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAAACCTGCGCCCTCTGGGCCACACGATCAGCGCACTTTTGCACctgc                          >  2:591684/1‑100 (MQ=255)
                                                                              ggtgagccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgt                     >  2:2836222/8‑100 (MQ=38)
                                                                              gCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCCGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgt                     >  2:3158778/1‑100 (MQ=255)
                                                                                tgagccAGGTTATCGCGGTAACCTGGCTCACCGCGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgtgt                   >  2:2405193/6‑100 (MQ=255)
                                                                                tgagccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgtgt                   >  1:1787104/6‑100 (MQ=255)
                                                                                tCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgtgt                   <  1:2287308/100‑1 (MQ=255)
                                                                                tCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgtg                    <  2:2277678/99‑1 (MQ=255)
                                                                                 gagccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTc                  >  1:1542915/5‑100 (MQ=255)
                                                                                 cACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTc                  <  2:3531625/100‑1 (MQ=255)
                                                                                  agccgggTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACGCTGTGCAAGACGATCAGGGCAGTTTTGTAACTGCTGGGGGTcg                 >  1:165294/6‑99 (MQ=255)
                                                                                   cGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCtt                <  1:1668504/100‑1 (MQ=255)
                                                                                    ccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCttt               >  1:44809/2‑100 (MQ=255)
                                                                                    ccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCttt               >  1:51674/2‑100 (MQ=255)
                                                                                    ccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCttt               >  1:1672916/2‑100 (MQ=255)
                                                                                     cAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCtttt              <  2:1303040/100‑1 (MQ=255)
                                                                                      aGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTa             <  2:343137/100‑1 (MQ=255)
                                                                                      aGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTa             <  1:1170661/100‑1 (MQ=255)
                                                                                       ggTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTaa            <  1:3632531/100‑1 (MQ=255)
                                                                                         ttACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAAtt          >  1:2098692/1‑100 (MQ=255)
                                                                                          tACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTc         <  1:2087853/100‑1 (MQ=255)
                                                                                           aCCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCg        >  1:1897972/1‑100 (MQ=255)
                                                                                           aCCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCg        <  1:1120848/100‑1 (MQ=255)
                                                                                            ccGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCgg       >  1:386184/1‑100 (MQ=255)
                                                                                              gcgATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCgg       <  2:1633344/98‑1 (MQ=255)
                                                                                              gcgATAACCTGGCTCACCGAGGTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTg     <  1:298899/100‑1 (MQ=255)
                                                                                               cgATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGc    <  2:1437586/100‑1 (MQ=255)
                                                                                                 aTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCgg  <  1:2394155/100‑1 (MQ=255)
                                                                                                 aTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCgg  <  1:2712754/100‑1 (MQ=255)
                                                                                                 aTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCgg  >  1:664126/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGG  >  NC_000913/1892266‑1892462

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