Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,263,512:1 +G 5.7% coding (184/1005 nt) asd → aspartate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,263,5121.G 5.7% 1476.8 / 75.6 336coding (184/1005 nt)asdaspartate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (156/161);  new base G (10/9);  total (166/170)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 350 total aligned reads displayed.

CACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAAC  >  NC_002947/2263364‑2263661
                                                                                                                                                     |                                                                                                                                                     
cACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCtt                                                                                                                                                      >  1:4507331/1‑150 (MQ=255)
 aCCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTc                                                                                                                                                     <  1:2535206/150‑1 (MQ=255)
 aCCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTc                                                                                                                                                     <  1:5982521/150‑1 (MQ=255)
 aCCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTc                                                                                                                                                     <  2:5668139/150‑1 (MQ=255)
  ccGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCg                                                                                                                                                    >  2:4399366/1‑150 (MQ=255)
        gCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCg                                                                                                                                              <  1:7086860/150‑1 (MQ=255)
        gCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAAGCTGTCGACTTCg                                                                                                                                              <  1:6505715/150‑1 (MQ=255)
         cGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAAGCTGTCGACTTCGc                                                                                                                                             <  1:7388034/150‑1 (MQ=255)
         cGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAAGCTGTCGACTTCGc                                                                                                                                             <  1:718269/150‑1 (MQ=255)
         cGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAAGCTGTCGACTTCGc                                                                                                                                             <  1:5317845/150‑1 (MQ=255)
          gTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGcc                                                                                                                                            >  2:8296402/1‑150 (MQ=255)
           tCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGctc                                                                                                                                           <  1:1392689/150‑3 (MQ=255)
             ggTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAg                                                                                                                                         <  1:7852724/150‑1 (MQ=255)
                   gCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAg                                                                                                                                   >  1:8740902/1‑150 (MQ=255)
                   gCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAg                                                                                                                                   >  2:734997/1‑150 (MQ=255)
                   gCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAg                                                                                                                                   >  1:4071040/1‑150 (MQ=255)
                     cTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCt                                                                                                                                 <  2:5755065/150‑1 (MQ=255)
                      ttGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTg                                                                                                                                >  2:4745973/1‑150 (MQ=255)
                      ttGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTg                                                                                                                                <  1:8154833/150‑1 (MQ=255)
                       tGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTgg                                                                                                                               <  1:4525626/150‑1 (MQ=255)
                        ggTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGc                                                                                                                              <  1:182967/150‑1 (MQ=255)
                        ggTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGc                                                                                                                              <  2:3887708/150‑1 (MQ=255)
                        ggTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGc                                                                                                                              <  1:4881477/150‑1 (MQ=255)
                           gCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCctt                                                                                                                           >  2:6194917/1‑150 (MQ=255)
                           gCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCctt                                                                                                                           >  1:1601019/1‑150 (MQ=255)
                            cAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCcttc                                                                                                                          <  1:789830/150‑1 (MQ=255)
                              gATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCcttctt                                                                                                                        <  2:6481392/150‑1 (MQ=255)
                               aTCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCttcttc                                                                                                                       <  2:2928420/150‑1 (MQ=255)
                               aTCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCttcttc                                                                                                                       <  2:3766640/150‑1 (MQ=255)
                                   tCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCa                                                                                                                   >  1:3837813/1‑150 (MQ=255)
                                              gCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCgccgcc                                                                                                        >  2:8460923/1‑150 (MQ=255)
                                                    ccGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGc                                                                                                  <  2:5411033/150‑1 (MQ=255)
                                                      ggTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCg                                                                                                >  1:4259313/1‑150 (MQ=255)
                                                             aCTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCttt                                                                                         >  1:3358272/1‑150 (MQ=255)
                                                              ctctGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTg                                                                                        <  2:789569/150‑1 (MQ=255)
                                                                 tGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCg                                                                                     <  1:2472859/150‑1 (MQ=255)
                                                                  gCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCgg                                                                                    >  1:4670551/1‑150 (MQ=255)
                                                                  gCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCgg                                                                                    >  2:153991/1‑150 (MQ=255)
                                                                  gCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCAGTCAGCCGCAGCTTTGCCgg                                                                                    >  1:5863407/1‑150 (MQ=255)
                                                                   cACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGc                                                                                   <  2:1742275/150‑1 (MQ=255)
                                                                     ccTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCaa                                                                                 >  1:5471623/1‑150 (MQ=255)
                                                                      ctgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAg                                                                                <  1:4617057/150‑1 (MQ=255)
                                                                      ctgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAg                                                                                <  2:2619204/150‑1 (MQ=255)
                                                                              cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCa                                                                        >  2:315994/1‑150 (MQ=255)
                                                                              cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCa                                                                        >  1:3579759/1‑150 (MQ=255)
                                                                                   tGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCTGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCgg                                                                   >  2:8097280/1‑150 (MQ=255)
                                                                                   tGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCgg                                                                   <  1:6753531/150‑1 (MQ=255)
                                                                                   tGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCgg                                                                   <  1:8607693/150‑1 (MQ=255)
                                                                                    ggAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCggg                                                                  <  2:4091403/150‑1 (MQ=255)
                                                                                        tCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGTGGGCTGc                                                              >  1:4845846/1‑150 (MQ=255)
                                                                                        tCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGc                                                              >  2:699585/1‑150 (MQ=255)
                                                                                        tCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGc                                                              >  2:4974344/1‑150 (MQ=255)
                                                                                          ggcggGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCAc                                                            >  1:2865160/1‑150 (MQ=255)
                                                                                           gcggGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACg                                                           <  1:615563/150‑1 (MQ=255)
                                                                                            cggGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACgg                                                          <  2:5392101/150‑1 (MQ=255)
                                                                                                cagcagCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCAt                                                      <  2:5538289/150‑1 (MQ=255)
                                                                                                 agcagcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATc                                                     <  1:3439887/150‑1 (MQ=255)
                                                                                                    agcGTGATGTTCGCTGGCAAGTAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                                  <  2:4230283/150‑1 (MQ=255)
                                                                                                    agcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                                  <  2:2538289/150‑1 (MQ=255)
                                                                                                    agcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                                  >  2:8720100/1‑150 (MQ=255)
                                                                                                    agcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                                  <  1:1998447/150‑1 (MQ=255)
                                                                                                    agcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                                  <  2:3561764/150‑1 (MQ=255)
                                                                                                    agcGTGATGGTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                                  <  2:6585043/150‑1 (MQ=255)
                                                                                                             gtCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCcggc                                         <  2:8779235/149‑1 (MQ=255)
                                                                                                               cGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCcggcgg                                       >  2:4012442/1‑150 (MQ=255)
                                                                                                                gCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCggcggc                                      >  1:5681744/1‑150 (MQ=255)
                                                                                                                  tGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCt                                    >  1:769464/1‑150 (MQ=255)
                                                                                                                  tGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCt                                    >  1:5160627/1‑150 (MQ=255)
                                                                                                                    gCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTgg                                  >  2:7664492/1‑150 (MQ=255)
                                                                                                                      aaGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGgac                                >  1:8425572/1‑150 (MQ=255)
                                                                                                                            cTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgcc                          <  2:3245006/150‑1 (MQ=255)
                                                                                                                                  gTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggcc                    <  1:2684114/150‑1 (MQ=255)
                                                                                                                                  gTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggcc                    >  1:1020422/1‑150 (MQ=255)
                                                                                                                                  ctgggcgAAGTCGAAGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggcc                    >  1:5402575/5‑150 (MQ=255)
                                                                                                                                   tGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccc                   <  1:3147135/150‑1 (MQ=255)
                                                                                                                                     ggcgAAGTCGAAGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccctg                 >  2:5107259/2‑150 (MQ=255)
                                                                                                                                      gcgAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGccctggccctgg                <  2:1346299/150‑1 (MQ=255)
                                                                                                                                      gcgAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGccctggccctgg                >  2:1048646/1‑150 (MQ=255)
                                                                                                                                      gcgAAGTCGACAGCTTCGACTTCGCCAAGGTCAAGCTGGCCTTCTTCGCCACCTGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGAGGCCTGGACGACGccctggccctgg                >  1:1185053/1‑150 (MQ=255)
                                                                                                                                        gAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTg              >  1:745989/1‑150 (MQ=255)
                                                                                                                                          aGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGcc            <  1:4928758/150‑1 (MQ=255)
                                                                                                                                          aGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGcc            <  2:5201179/150‑1 (MQ=255)
                                                                                                                                            tCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCgg          >  2:4498769/1‑150 (MQ=255)
                                                                                                                                            tCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCgg          >  1:2154211/1‑150 (MQ=255)
                                                                                                                                            tCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCgg          >  2:1828725/1‑150 (MQ=255)
                                                                                                                                              gACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa        <  2:7052076/150‑1 (MQ=255)
                                                                                                                                              gACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa        >  1:7017765/1‑150 (MQ=255)
                                                                                                                                              gACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa        <  2:7599826/150‑1 (MQ=255)
                                                                                                                                              gACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa        <  1:8378797/150‑1 (MQ=255)
                                                                                                                                               aCAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAg       >  2:1451754/1‑150 (MQ=255)
                                                                                                                                                cAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGc      <  1:6189728/150‑1 (MQ=255)
                                                                                                                                                   cTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCaa   >  2:2372427/1‑150 (MQ=255)
                                                                                                                                                   cTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCaa   >  1:3800657/1‑150 (MQ=255)
                                                                                                                                                   cTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCaa   >  2:3472416/1‑150 (MQ=255)
                                                                                                                                                    ttCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCa                                                                                  <  1:7393522/70‑1 (MQ=255)
                                                                                                                                                    ttCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCa                                                                                  >  2:7393522/1‑70 (MQ=255)
                                                                                                                                                    ttCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAAc  >  2:6038206/1‑150 (MQ=255)
                                                                                                                                                    ttCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAAc  <  1:6451743/150‑1 (MQ=255)
                                                                                                                                                    ttCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAAc  >  1:2064554/1‑150 (MQ=255)
                                                                                                                                                    ttCGACCTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAAc  >  1:7501047/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAAC  >  NC_002947/2263364‑2263661

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