Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 244,809 A→T 6.6% S161S (TCA→TCT)  yafJ → type 2 glutamine amidotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913244,8090AT 6.6% 229.0 / 3.9 83S161S (TCA→TCT) yafJtype 2 glutamine amidotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (42/33);  new base T (3/3);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTG  >  NC_000913/244662‑244953
                                                                                                                                                   |                                                                                                                                                
gggCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCAc                                                                                                                                                 <  1:216274/149‑1 (MQ=255)
gggCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCAc                                                                                                                                                 <  1:325833/149‑1 (MQ=255)
  gCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGTAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCAGAAATTAACGCAGCGTTACCCGCGCACACCGGGGAACTTGGCGGCGGTATTTTAATATTTCGCCTCACTg                                                                                                                                               >  2:103852/1‑149 (MQ=255)
  gCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAACGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAGCATGGCGGCGGTATTTAAATATATCGCCTCACTg                                                                                                                                               >  1:570347/1‑149 (MQ=255)
   cTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCCCttg                                                                                                                                              >  1:451885/1‑147 (MQ=255)
   cTACAAATCACTGGAAACCGGAAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAATCCTTATGATGGCTACTGCTTAAATTAACGCAGCGTAACCCGCGCACACCGGGCAACATGGAGGCGGTATTTAAATATATCGCCTCACTgg                                                                                                                                              >  1:1059872/1‑149 (MQ=255)
     aCAAATCACTGGAAACCGTCAACTTCCGCCCGGTAGTCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCTGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCg                                                                                                                                            <  1:411409/149‑1 (MQ=255)
     aCAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCg                                                                                                                                            <  1:405620/149‑1 (MQ=255)
     aCAAATCACTGGAAACCGGCAACTACCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCCCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATCGCGGCGGTATTTAAATATATCGCCTCACTGGCg                                                                                                                                            >  1:1133670/1‑149 (MQ=255)
       aaaTCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGa                                                                                                                                          <  1:881154/149‑1 (MQ=255)
        aaTCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGAt                                                                                                                                         <  1:1590938/149‑1 (MQ=255)
          tCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGa                                                                                                                                       >  1:611010/1‑149 (MQ=255)
           cACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGaa                                                                                                                                      <  2:1203206/149‑1 (MQ=255)
           cACTGGAAACCGGCAACTTCCCCCCGGAAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAAAGCACCGATACCCGCGCACACCGGGCAACATCGCGGCGGTATTTAAATATATCGCAACACTGGCGGATGaa                                                                                                                                      >  2:334819/1‑149 (MQ=255)
                    ccGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAg                                                                                                                             >  1:461955/1‑149 (MQ=255)
                     cGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGa                                                                                                                            >  1:530513/1‑149 (MQ=255)
                     cGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGa                                                                                                                            >  1:811820/1‑149 (MQ=255)
                       gCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAg                                                                                                                          <  2:1205507/149‑1 (MQ=255)
                       gCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAg                                                                                                                          <  2:1211944/149‑1 (MQ=255)
                           cTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCg                                                                                                                      <  2:1326531/149‑1 (MQ=255)
                            ttCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGt                                                                                                                     >  1:1247637/1‑149 (MQ=255)
                              ccGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGttt                                                                                                                   <  2:535837/149‑1 (MQ=255)
                               cGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGtttt                                                                                                                  >  2:858272/1‑149 (MQ=255)
                                gCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCTCTGGCGGATGAACTGCGGCAGAAGGCCGttttg                                                                                                                 >  1:120399/1‑148 (MQ=255)
                                 cccGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGTGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCa                                                                                                                >  2:988044/1‑149 (MQ=255)
                                       aGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTTCCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACAtgc                                                                                                          <  2:44546/149‑1 (MQ=255)
                                         gCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACAtgctg                                                                                                        <  1:181478/149‑1 (MQ=255)
                                           gAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATgctgct                                                                                                      >  1:1547308/1‑149 (MQ=255)
                                            taaCCGACAGCGAAAAATCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCtt                                                                                                     <  1:631679/148‑1 (MQ=255)
                                             aaCCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCttt                                                                                                    <  2:19694/149‑1 (MQ=255)
                                                cGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCgg                                                                                                 <  1:775725/149‑1 (MQ=255)
                                                  aCAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGAc                                                                                               <  2:1247637/149‑1 (MQ=255)
                                                    aGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACgg                                                                                             >  2:1447639/1‑149 (MQ=255)
                                                     gCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACggg                                                                                            >  1:784992/1‑149 (MQ=255)
                                                      cGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAGGGGCGTTTTCAACATGCTGCTTTCGGACGGgc                                                                                           >  2:72931/1‑149 (MQ=255)
                                                       gAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGAAGAACTGCGGCAGAATGGCGTTTTCAACATGCTGCTTGCGGACGGgcg                                                                                          >  1:590505/1‑149 (MQ=255)
                                                        aaaaaGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGgcgc                                                                                         >  1:119159/1‑149 (MQ=255)
                                                            aGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATg                                                                                     <  1:530033/149‑1 (MQ=255)
                                                            aGCCTGCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGTGCAACATGGCGGCGGTATTTAAATAAATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGGTTTCAACATGCTGCTTTCGGACGGGCGCTATg                                                                                     <  2:241430/149‑1 (MQ=255)
                                                             gCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCTCCTCAATGGCGGATGAACTGCGGCAGAAGGGCGTTTTCATCATGCTGCTTTCGGACGGGCGCTATGt                                                                                    >  1:62721/1‑149 (MQ=255)
                                                                atCTGCGGGCTCCAGCATAAAATAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAAt                                                                                 <  2:758954/148‑1 (MQ=255)
                                                                 tCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATATAAATATATCGCCTCCCTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCTACATGCTGCTTACGGACGGGCGCAATGTATTg                                                                                >  1:290566/1‑149 (MQ=255)
                                                                  ctgctgGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATgg                                                                               >  1:343136/1‑149 (MQ=255)
                                                                        gCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTAtt                                                                         <  2:811820/149‑1 (MQ=255)
                                                                          tCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGc                                                                       >  2:1072136/1‑149 (MQ=255)
                                                                              gCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGa                                                                   <  2:611010/149‑1 (MQ=255)
                                                                                aTAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCTCTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACt                                                                 >  2:500782/1‑149 (MQ=255)
                                                                                aTAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCTCTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACt                                                                 >  2:506897/1‑149 (MQ=255)
                                                                                aTAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGACCTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACt                                                                 >  2:1328641/1‑149 (MQ=255)
                                                                                aTAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACt                                                                 >  2:539310/1‑149 (MQ=255)
                                                                                 tAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTa                                                                >  1:118036/1‑149 (MQ=255)
                                                                                 tAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTa                                                                >  1:1532662/1‑149 (MQ=255)
                                                                                   aaTTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAAt                                                              >  2:1590451/1‑149 (MQ=255)
                                                                                    aTTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATTGCGTATTGCTCGACTAAtt                                                             >  2:1021430/1‑149 (MQ=255)
                                                                                    aTTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAAtt                                                             >  1:875505/1‑149 (MQ=255)
                                                                                    aTTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAAtt                                                             >  1:1171143/1‑149 (MQ=255)
                                                                                    aTTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGAAGAAGGGCGTTTTCAACATGCTGATTTCGGAAGGGCGCTATGTAATGGCGTATTTCTCGACTAAtt                                                             >  2:172188/1‑149 (MQ=255)
                                                                                     ttAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAAttt                                                            >  2:1029028/1‑149 (MQ=255)
                                                                                     ttAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAAttt                                                            >  1:783424/1‑149 (MQ=255)
                                                                                      tAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCGCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTa                                                           <  2:1406665/149‑1 (MQ=255)
                                                                                          gCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACt                                                       >  2:301342/1‑149 (MQ=255)
                                                                                                          acacCGGGCAACATGGCGGTGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCcgcgc                                       <  1:988044/149‑1 (MQ=255)
                                                                                                            acCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCgcgcgc                                     >  1:1056675/1‑149 (MQ=255)
                                                                                                               gggCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGt                                  <  2:570347/149‑1 (MQ=255)
                                                                                                                gccAACATGGCGGCGGTATTTAAATATATCGCCGCTCTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGtt                                 <  1:103852/147‑1 (MQ=255)
                                                                                                                 gCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGttt                                <  2:192918/149‑1 (MQ=255)
                                                                                                                   aaCATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTgg                              <  2:414660/149‑1 (MQ=255)
                                                                                                                   aaCATGGCGGCGGTATTTAAATATATCGCCTCACTAGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTgg                              >  1:788712/1‑149 (MQ=255)
                                                                                                                      aTGGCGGCGGTATTTAAATATATCGCCTCTCTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGt                           <  1:500782/149‑1 (MQ=255)
                                                                                                                      aTGGCGGCGGTATTTAAATATATCGCCTCTCTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGt                           <  1:506897/149‑1 (MQ=255)
                                                                                                                       ttgcggcgGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGGAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTg                          <  2:120399/147‑1 (MQ=255)
                                                                                                                        ggaggcggTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTTCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTgg                         <  1:172188/146‑1 (MQ=255)
                                                                                                                              ggTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGt                   >  2:628813/1‑149 (MQ=255)
                                                                                                                               gTATTTAAATATATCGCCTCACTGGAGGATGACCTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGtt                  >  1:1506989/1‑149 (MQ=255)
                                                                                                                                tATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTg                 >  1:681276/1‑149 (MQ=255)
                                                                                                                                tATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTg                 <  1:1249994/149‑1 (MQ=255)
                                                                                                                                  tttAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCt               <  2:530513/149‑1 (MQ=255)
                                                                                                                                  tttAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCt               <  1:301342/149‑1 (MQ=255)
                                                                                                                                  tttAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCt               >  1:140330/1‑149 (MQ=255)
                                                                                                                                    tAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTgg             >  2:1170747/1‑149 (MQ=255)
                                                                                                                                            tCGCCTCACTGGCGGATGAACTGCGGCAGTAGTGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGAt     <  2:590505/149‑1 (MQ=255)
                                                                                                                                               ccTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGAtgtg  >  2:449514/1‑149 (MQ=255)
                                                                                                                                               ccTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGAtgtg  >  2:1259707/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTG  >  NC_000913/244662‑244953

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