Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 2,073,851 Δ1 bp intergenic (+109/+159) lysV → / ← xapR tRNA‑Lys/DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,073,8510T.100.0% 437.4 / NA 96intergenic (+109/+159)lysV/xapRtRNA‑Lys/DNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (21/75);  total (21/75)

GAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGTGCTGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTT  >  AP012306/2073798‑2073987
                                                     |                                                                                                                                        
gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGa                                       >  1:3631357/1‑150 (MQ=0)
 aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGTCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      <  2:3756806/150‑1 (MQ=11)
 aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      >  1:4132036/1‑150 (MQ=2)
 aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      <  2:126339/150‑1 (MQ=2)
 aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      <  1:679455/150‑1 (MQ=2)
  aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa                                     <  2:928083/150‑1 (MQ=255)
  aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa                                     >  2:2990189/1‑150 (MQ=255)
  aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa                                     >  1:3773502/1‑150 (MQ=255)
   ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    <  1:2737544/150‑1 (MQ=255)
   ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    <  1:839064/150‑1 (MQ=255)
    cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   <  2:1441591/150‑1 (MQ=255)
    cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   <  1:3567416/150‑1 (MQ=255)
    cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   <  1:3027456/150‑1 (MQ=255)
    cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   >  1:2843311/1‑150 (MQ=255)
    cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   <  2:1891681/150‑1 (MQ=255)
    cTGCAGCAGGTTTGAGTGGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   <  2:776387/150‑1 (MQ=11)
     tGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTg                                  >  2:2407400/1‑150 (MQ=255)
     tGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTg                                  <  2:2514048/150‑1 (MQ=255)
     tGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTg                                  <  1:1453377/150‑1 (MQ=255)
      gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  1:1245411/150‑1 (MQ=255)
      gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  2:990427/150‑1 (MQ=255)
      gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  1:3034885/150‑1 (MQ=255)
       cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGcc                                <  1:1639380/150‑1 (MQ=255)
       cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGcc                                <  1:2120409/150‑1 (MQ=255)
       cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGcc                                <  2:837332/150‑1 (MQ=255)
        agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               >  1:4047074/1‑150 (MQ=255)
         gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTg                              <  1:403242/150‑1 (MQ=255)
          cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                             <  2:4229843/150‑1 (MQ=255)
          cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                             >  1:493343/1‑150 (MQ=255)
          cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                             <  1:2576742/150‑1 (MQ=255)
           agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAGGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGAt                            <  1:4258212/150‑1 (MQ=11)
           agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGAt                            <  2:1191603/150‑1 (MQ=255)
            ggTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATg                           <  2:950545/150‑1 (MQ=255)
            ggTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATg                           <  2:1009045/150‑1 (MQ=255)
             gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          >  2:451887/1‑150 (MQ=255)
             gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          <  2:268081/150‑1 (MQ=255)
             gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          >  1:3071859/1‑150 (MQ=255)
              tttGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                         >  1:1546392/1‑150 (MQ=255)
               ttgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcgc                        <  2:2986939/150‑1 (MQ=255)
               ttgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcgc                        <  2:3862776/150‑1 (MQ=255)
               ttgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcgc                        <  2:655183/150‑1 (MQ=255)
               ttgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcgc                        <  1:2190540/150‑1 (MQ=255)
                tgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCTCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATTCGCt                       >  2:586617/1‑150 (MQ=14)
                 gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      <  2:1405517/150‑1 (MQ=255)
                 gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      <  2:1319646/150‑1 (MQ=255)
                 gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      >  2:4374376/1‑150 (MQ=255)
                  agttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTc                     >  2:4399243/1‑150 (MQ=255)
                  agttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTc                     <  1:2407400/150‑1 (MQ=255)
                  agttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTc                     >  2:2794068/1‑150 (MQ=255)
                  agttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTc                     <  2:328732/150‑1 (MQ=255)
                   gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    <  1:4539333/150‑1 (MQ=255)
                   gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    >  2:2050374/1‑150 (MQ=255)
                    ttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCAc                   <  2:1406926/150‑1 (MQ=255)
                    ttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCAc                   >  1:164516/1‑150 (MQ=255)
                    ttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCAc                   <  1:3367522/150‑1 (MQ=255)
                    ttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCAc                   <  2:2167847/150‑1 (MQ=255)
                     tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  2:3217658/150‑1 (MQ=255)
                     tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  1:3237051/150‑1 (MQ=255)
                     tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  2:3968136/150‑1 (MQ=255)
                     tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  2:1424797/150‑1 (MQ=255)
                       agCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                <  1:1740465/150‑1 (MQ=255)
                       agCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                <  2:4359389/150‑1 (MQ=255)
                       agCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                >  2:2449851/1‑150 (MQ=255)
                        gcgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCTGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  2:1335594/150‑1 (MQ=17)
                        gcgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  1:1181883/150‑1 (MQ=255)
                        gcgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  1:4258957/150‑1 (MQ=255)
                        gcgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  2:1924873/150‑1 (MQ=255)
                        gcgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  1:4335360/150‑1 (MQ=255)
                         cgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGc              <  1:2426035/150‑1 (MQ=255)
                          gcAGCGAAACAACGGAGCCGCTC‑‑GCGGGGGCGGCCCGAAGGGCGAGGGAAGCGAGGCATCCTGCACGACGCACCAATGGAAAAAAGCGCCCTAAAGCGGCTTTTTTGCTATCTGCGAATTGCGAAATTGCCTGATGCGCTTCACTTATCa             >  2:893692/1‑150 (MQ=14)
                          gcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             >  2:3601774/1‑150 (MQ=255)
                          gcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             <  2:3609297/150‑1 (MQ=255)
                           cAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            <  1:310626/150‑1 (MQ=255)
                           cAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            <  1:1395112/150‑1 (MQ=255)
                           cAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            <  2:3036279/150‑1 (MQ=255)
                            aGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGa           <  2:624938/150‑1 (MQ=255)
                             gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          <  1:988124/150‑1 (MQ=255)
                             gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          <  1:4400320/150‑1 (MQ=255)
                             gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          <  2:2142228/150‑1 (MQ=255)
                              cGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGact         <  2:2357713/150‑1 (MQ=255)
                               gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAATAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        >  2:2491154/1‑150 (MQ=17)
                               gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  1:515722/150‑1 (MQ=255)
                               gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  2:3334214/150‑1 (MQ=255)
                                aaaCAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       <  2:2923146/150‑1 (MQ=255)
                                 aacaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactact      <  1:3478526/150‑1 (MQ=255)
                                 aacaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactact      <  2:751266/150‑1 (MQ=255)
                                  acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  2:1710776/150‑1 (MQ=255)
                                   caacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    <  2:3676103/150‑1 (MQ=255)
                                   caacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    <  1:371796/150‑1 (MQ=255)
                                   caacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    >  2:748269/1‑150 (MQ=255)
                                   caacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    <  1:1505253/150‑1 (MQ=255)
                                    aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  1:152166/150‑1 (MQ=255)
                                    aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  2:1316954/150‑1 (MQ=255)
                                    aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  2:1971073/150‑1 (MQ=255)
                                    aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  2:4230789/150‑1 (MQ=255)
                                     acGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  1:2799052/150‑1 (MQ=255)
                                     acGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  2:2200828/150‑1 (MQ=255)
                                     acGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  2:3490536/150‑1 (MQ=255)
                                                     |                                                                                                                                        
GAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGTGCTGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTT  >  AP012306/2073798‑2073987

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