Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,486,574 C→A A152E (GCA→GAA)  crp → cAMP‑activated global transcription factor, mediator of catabolite repression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,486,5740CA97.5% 407.1 / 2.1 118A152E (GCA→GAA) crpcAMP‑activated global transcription factor, mediator of catabolite repression
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (55/60);  total (56/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 118 total aligned reads displayed.

GCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACG  >  NC_000913_3_pae_pgi/3486430‑3486720
                                                                                                                                                |                                                                                                                                                  
gCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaac                                                                                                                                               >  1:821870/1‑150 (MQ=255)
gCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaac                                                                                                                                               <  1:377004/150‑1 (MQ=255)
gCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTATGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaac                                                                                                                                               >  3:522088/1‑150 (MQ=255)
 ccAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaaca                                                                                                                                              >  3:491055/1‑150 (MQ=255)
 ccAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaaca                                                                                                                                              <  3:861786/150‑1 (MQ=255)
   aaTTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaacaac                                                                                                                                            <  3:455408/150‑1 (MQ=255)
      tGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGGACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAg                                                                                                                                         <  4:912716/150‑1 (MQ=255)
         ttCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACg                                                                                                                                      >  1:834722/1‑150 (MQ=255)
          tCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGc                                                                                                                                     >  2:309485/1‑150 (MQ=255)
          tCAGGTAAACCCGGACATTCTGATGCGTTTATCTGCACAGATGGCGCGTCGTCTGAAAGTCACTTGAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCGGCATTGAACAGACTCTGCGTAACATGGAAAATCACCGAGACGc                                                                                                                                     >  2:829480/1‑150 (MQ=255)
            aGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTa                                                                                                                                   >  2:16029/1‑150 (MQ=255)
            aGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTa                                                                                                                                   <  4:165720/150‑1 (MQ=255)
             ggTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAt                                                                                                                                  >  2:367147/1‑150 (MQ=255)
               tAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGa                                                                                                                                <  1:1148740/150‑1 (MQ=255)
               tAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGa                                                                                                                                <  4:383512/150‑1 (MQ=255)
                   cccGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCa                                                                                                                            <  1:1533667/150‑1 (MQ=255)
                   cccGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCa                                                                                                                            <  2:658416/150‑1 (MQ=255)
                    ccGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCAc                                                                                                                           >  3:538771/1‑150 (MQ=255)
                           ttCTGCTGCGTTTGTCTGCACAGATGAGGTGTCGTATGCAAGTAACTTCAGAGAAAGCGGCCAACCTGGCGTGCCGAGACGTGAGGGGCCGCACTGCACAGACTCCCCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACg                                                                                                                    <  4:136719/150‑1 (MQ=255)
                            tCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACgg                                                                                                                   <  2:354409/150‑1 (MQ=255)
                              tGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCg                                                                                                                        <  1:936479/143‑1 (MQ=255)
                              tGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCg                                                                                                                        >  2:936479/1‑143 (MQ=255)
                              tGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTa                                                                                                                 >  2:420342/1‑150 (MQ=255)
                               gATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTAt                                                                                                                >  3:992582/1‑150 (MQ=255)
                                aTGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATg                                                                                                               >  3:449934/1‑150 (MQ=255)
                                aTGCGTTTGTCTGCACAGAGGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATg                                                                                                               <  2:1468744/150‑1 (MQ=255)
                                   cGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCaa                                                                                                            <  4:182034/150‑1 (MQ=255)
                                    gTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCaaa                                                                                                           <  2:648084/150‑1 (MQ=255)
                                    gTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACGTTAGATAAAGTTTTCAACCTGGCGTTACTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCaaa                                                                                                           >  2:875576/1‑150 (MQ=255)
                                            gCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAAtt                                                                                                   >  3:195522/1‑150 (MQ=255)
                                             cacaGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTa                                                                                                  <  1:1429949/150‑1 (MQ=255)
                                                aGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTAccc                                                                                               <  4:877609/150‑1 (MQ=255)
                                                  aTGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaacaa                                                                                                                                             <  4:605749/102‑1 (MQ=255)
                                                  aTGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaacaa                                                                                                                                             >  3:605749/1‑102 (MQ=255)
                                                       gcgTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCgg                                                                                                                       <  4:202985/119‑1 (MQ=255)
                                                       gcgTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCgg                                                                                                                       >  3:202985/1‑119 (MQ=255)
                                                       gcgTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAAGTACCCGTCAGGa                                                                                        >  2:682441/1‑150 (MQ=255)
                                                             tcTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTgg                                                                                  >  2:813011/1‑150 (MQ=255)
                                                              cTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGt                                                                                 >  4:511326/1‑150 (MQ=255)
                                                               tGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTc                                                                                <  1:126141/150‑1 (MQ=255)
                                                                gCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCa                                                                               >  4:308315/1‑150 (MQ=255)
                                                                 cAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAg                                                                              >  4:925724/1‑150 (MQ=255)
                                                                  aaGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGGCTCTGCTGAATCTGCAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAAACAAAATTACCCGTCAGGAAATTGGTCAGa                                                                             <  4:121033/150‑1 (MQ=255)
                                                                      cACTTCAGAGAGAGTGGGCGACCGGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAATCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGt                                                                         <  2:267550/150‑1 (MQ=255)
                                                                         ttCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                      >  4:921294/1‑150 (MQ=255)
                                                                         ttCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                      <  1:823528/150‑1 (MQ=255)
                                                                         ttCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                      >  1:1163766/1‑150 (MQ=255)
                                                                         ttCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                      <  1:1159305/150‑1 (MQ=255)
                                                                         ttCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                      >  1:1103528/1‑150 (MQ=255)
                                                                            agagAAAGTGGGCAACCTGGCGTTCCTCGAGGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAGAAAAGACGAGTGGCTATGACCCACCTGGAGGGGATGCGAATCAAAATTACTCCCCAGGAAATTGGTCAGGTTGTCGGCTg                                                                   >  1:152201/1‑150 (MQ=255)
                                                                              agaAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCaa                                                                                                            <  2:1418545/107‑1 (MQ=255)
                                                                              agaAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGtt                                                                 >  3:508651/1‑150 (MQ=255)
                                                                                aaaGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTtct                                                               >  2:131091/1‑150 (MQ=255)
                                                                                   gTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGACATTGGTCAGATTGTCGGCTGTTCTCGt                                                            >  4:859829/1‑150 (MQ=255)
                                                                                            cTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGGTCCGGCTGAATTGGGAAAAACAACCAGACGCTATGACTCAACCGGCCGGTATGCAAATCAAAATTACACGTCCGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTg                                                   >  4:800316/1‑150 (MQ=255)
                                                                                            cTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTg                                                   <  2:691883/150‑1 (MQ=255)
                                                                                               gCGTGCCCCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGGCGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGa                                                <  2:1245715/150‑1 (MQ=255)
                                                                                                 gTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGAc                                                                                                                               <  3:178453/69‑1 (MQ=255)
                                                                                                     cTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                      >  4:706634/1‑122 (MQ=255)
                                                                                                     cTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                      <  3:706634/122‑1 (MQ=255)
                                                                                                      tCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTc                                         >  3:601224/1‑150 (MQ=255)
                                                                                                       cGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCt                                        <  4:522088/150‑1 (MQ=255)
                                                                                                          cGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGaa                                     <  1:1056067/150‑1 (MQ=255)
                                                                                                            tGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGa                                   >  1:1132105/1‑150 (MQ=255)
                                                                                                             gACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGCCTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCTGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGAt                                  <  2:949784/150‑1 (MQ=255)
                                                                                                             gACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGAt                                  <  1:813011/150‑1 (MQ=255)
                                                                                                                ggtccGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCt                               >  2:1092347/4‑150 (MQ=255)
                                                                                                                gggCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAACCCGTGGGACGCCTTCTGAAGATGCt                               >  4:367468/1‑150 (MQ=255)
                                                                                                                 ggCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGCCTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTg                              <  1:1273871/150‑1 (MQ=255)
                                                                                                                 ggCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTg                              <  3:921294/150‑1 (MQ=255)
                                                                                                                  gCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTgg                             <  3:568258/150‑1 (MQ=255)
                                                                                                                   ccGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTgg                                                  >  4:132237/1‑128 (MQ=255)
                                                                                                                   ccGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTgg                                                  <  3:132237/128‑1 (MQ=255)
                                                                                                                    cGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGaa                           >  1:881470/1‑150 (MQ=255)
                                                                                                                    cGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGaa                           <  3:882318/150‑1 (MQ=255)
                                                                                                                    cGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGaa                           <  3:517956/150‑1 (MQ=255)
                                                                                                                     gCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAg                          <  4:538771/150‑1 (MQ=255)
                                                                                                                       aTTGCACAGACTCTGCTGAAGCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGAt                        <  2:185897/150‑1 (MQ=255)
                                                                                                                         tGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCa                      <  1:636432/150‑1 (MQ=255)
                                                                                                                           tacaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGa                    >  2:1190573/2‑150 (MQ=255)
                                                                                                                           cacaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGt                                                            <  1:1479576/110‑1 (MQ=255)
                                                                                                                            acaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAt                                                                                                                                  <  4:140516/39‑1 (MQ=255)
                                                                                                                            acaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAt                                                                                                                                  >  3:140516/1‑39 (MQ=255)
                                                                                                                            acaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGaa                   >  2:694320/1‑150 (MQ=255)
                                                                                                                            acaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGaa                   >  1:202047/1‑150 (MQ=255)
                                                                                                                             caGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAAc                  <  3:221945/150‑1 (MQ=255)
                                                                                                                                 ctctGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGa              <  2:288297/150‑1 (MQ=255)
                                                                                                                                  tctGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGAt             >  2:290138/1‑150 (MQ=255)
                                                                                                                                  tctGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGAt             <  1:420342/150‑1 (MQ=255)
                                                                                                                                   ctgctgAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGAtc            >  2:1576341/1‑150 (MQ=255)
                                                                                                                                      ctgAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTcc         <  3:76595/150‑1 (MQ=255)
                                                                                                                                      ctgAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTcc         <  1:1190573/150‑1 (MQ=255)
                                                                                                                                      ctgAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTcc         <  2:202047/150‑1 (MQ=255)
                                                                                                                                       tgAATCTGGAAAAACACCCAGACGCTCTGACTCGCCCGGACGGTATGCAAATCAGACTCGCCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCg        <  1:1357684/150‑1 (MQ=255)
                                                                                                                                       tgAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCg        <  2:821870/150‑1 (MQ=255)
                                                                                                                                           tCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAAc                  <  2:1246310/136‑1 (MQ=255)
                                                                                                                                           tCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAAc                  >  1:1246310/1‑136 (MQ=255)
                                                                                                                                           tCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGcaca    >  1:1206648/1‑150 (MQ=255)
                                                                                                                                           tCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGcaca    >  3:503695/1‑150 (MQ=255)
                                                                                                                                             tGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTAGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACg  >  1:86080/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACG  >  NC_000913_3_pae_pgi/3486430‑3486720

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: