Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,550 C→A 8.1% R174L (CGT→CTT)  pepE ← (alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,219,5500CA 8.1% 1068.0 / 41.2 288R174L (CGT→CTT) pepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (134/130);  new base A (14/10);  total (148/140)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

CCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGG  >  CP009273/4219401‑4219699
                                                                                                                                                     |                                                                                                                                                     
 cAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACg                                                                                                                                                      >  2:3046825/1‑150 (MQ=255)
     ggAACCGCTTCTTCACCAGCCTTAAACACATAATTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACAGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCAGGCGCGACGACCAGCAGTTCGCGAATACgctgc                                                                                                                                                  >  1:4210914/1‑150 (MQ=255)
     ggAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgc                                                                                                                                                  >  1:2120618/1‑150 (MQ=255)
      gAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgctgct                                                                                                                                                 >  1:2465556/1‑150 (MQ=255)
       aaCCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTc                                                                                                                                                >  2:3307484/1‑150 (MQ=255)
       aaCCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTc                                                                                                                                                >  2:1032105/1‑150 (MQ=255)
       aaCCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTc                                                                                                                                                <  1:4466056/150‑1 (MQ=255)
          cGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACg                                                                                                                                             <  1:1236097/150‑1 (MQ=255)
          cGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACg                                                                                                                                             <  2:4509922/150‑1 (MQ=255)
           gCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACgg                                                                                                                                            >  2:3190062/1‑150 (MQ=255)
           gCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACgg                                                                                                                                            <  1:3194291/150‑1 (MQ=255)
               cttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGttt                                                                                                                                        >  1:4276869/1‑150 (MQ=255)
                     cAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCtt                                                                                                                                  <  2:681694/150‑1 (MQ=255)
                      aGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCttt                                                                                                                                 <  1:2897768/150‑1 (MQ=255)
                       gCCTTAAACACATCAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTa                                                                                                                                <  2:722624/150‑1 (MQ=255)
                       gCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTa                                                                                                                                <  2:3010595/150‑1 (MQ=255)
                         cTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATg                                                                                                                              >  2:1537713/1‑150 (MQ=255)
                          ttAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATgg                                                                                                                             >  1:3079091/1‑150 (MQ=255)
                                acaTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTcc                                                                                                                       >  2:4227160/1‑150 (MQ=255)
                                        ggtgTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACATTTATGGCCTTCCGGCAgcgc                                                                                                               >  2:596412/1‑150 (MQ=255)
                                        ggtgTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAgcgc                                                                                                               >  2:501569/1‑150 (MQ=255)
                                         gtgtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGcgcg                                                                                                              >  1:2309186/1‑150 (MQ=255)
                                         gtgtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGcgcg                                                                                                              >  1:2309067/1‑150 (MQ=255)
                                             tCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTggg                                                                                                          >  2:3520203/1‑150 (MQ=255)
                                             tCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTgg                                                                                                          >  2:1248337/1‑150 (MQ=255)
                                             tCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTgg                                                                                                          >  2:2374711/1‑150 (MQ=255)
                                               gggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTg                                                                                                        >  2:1403634/1‑150 (MQ=255)
                                               gggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTg                                                                                                        >  2:3277671/1‑150 (MQ=255)
                                                 gCCACCCAGCGTAGCGTGACCATTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCGTCCGGCAGCGCGTGGGTGaa                                                                                                      <  1:404149/150‑1 (MQ=255)
                                                  ccACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAg                                                                                                     >  2:4353341/1‑150 (MQ=255)
                                                   cACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGt                                                                                                    <  1:2277936/150‑1 (MQ=255)
                                                    aCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTg                                                                                                   >  1:4280682/1‑150 (MQ=255)
                                                    aCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTg                                                                                                   <  2:2958268/150‑1 (MQ=255)
                                                    aCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCGTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTg                                                                                                   <  2:1898973/150‑1 (MQ=255)
                                                       cAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCgg                                                                                                >  2:4473296/1‑150 (MQ=255)
                                                        agcgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCggg                                                                                               >  2:949526/1‑150 (MQ=255)
                                                        agcgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCggg                                                                                               >  2:1958446/1‑150 (MQ=255)
                                                          cgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGtt                                                                                             <  1:1830970/150‑1 (MQ=255)
                                                          cgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGtt                                                                                             <  1:205359/150‑1 (MQ=255)
                                                            tagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGa                                                                                           >  2:1807815/1‑150 (MQ=255)
                                                            tagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGACGTGCGGGTTGa                                                                                           >  2:2770315/1‑150 (MQ=255)
                                                            tagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTGGa                                                                                           >  2:2769908/1‑150 (MQ=255)
                                                             agcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAt                                                                                          >  1:1905292/1‑150 (MQ=255)
                                                             agcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAt                                                                                          >  1:735701/1‑150 (MQ=255)
                                                               cgtGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGAttt                                                                                        <  1:92724/150‑1 (MQ=255)
                                                                 tGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGc                                                                                      >  2:1781665/1‑150 (MQ=255)
                                                                   aCCTTTGCTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAg                                                                                    <  1:2116716/150‑1 (MQ=255)
                                                                           tCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACa                                                                            <  2:419540/150‑1 (MQ=255)
                                                                           tCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACa                                                                            <  2:2710155/150‑1 (MQ=255)
                                                                            cACTGTGATCCAGTTACCTTCCGGTAGACCAATACTCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAg                                                                           <  1:1023606/150‑1 (MQ=255)
                                                                            cACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAg                                                                           <  2:1620295/150‑1 (MQ=255)
                                                                             aCTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGa                                                                          >  1:745104/1‑150 (MQ=255)
                                                                             aCTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGa                                                                          >  1:3025624/1‑150 (MQ=255)
                                                                               tgtgATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGAtt                                                                        >  1:2979300/1‑150 (MQ=255)
                                                                                gtgATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGAttt                                                                       <  1:1486122/150‑1 (MQ=255)
                                                                                   aTCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgc                                                                    <  1:589321/150‑1 (MQ=255)
                                                                                   aTCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgc                                                                    >  2:4155126/1‑150 (MQ=255)
                                                                                   aTCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgc                                                                    >  1:251304/1‑150 (MQ=255)
                                                                                   aTCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgc                                                                    <  1:3603304/150‑1 (MQ=255)
                                                                                    tCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgcg                                                                   >  2:2210186/1‑150 (MQ=255)
                                                                                    tCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgcg                                                                   >  2:4447221/1‑150 (MQ=255)
                                                                                    tCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAgcg                                                                   <  2:2878785/150‑1 (MQ=255)
                                                                                      cAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCa                                                                 >  2:2099423/1‑150 (MQ=255)
                                                                                        gTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATc                                                               >  2:656412/1‑150 (MQ=255)
                                                                                        gTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATc                                                               >  2:2919991/1‑150 (MQ=255)
                                                                                         ttACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCg                                                              <  2:3069763/150‑1 (MQ=255)
                                                                                             cTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAAc                                                          <  1:3289687/150‑1 (MQ=255)
                                                                                                    tAGACCAATACTCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCgg                                                   <  1:1268790/150‑1 (MQ=255)
                                                                                                      gACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGAt                                                 <  1:1610951/150‑1 (MQ=255)
                                                                                                      gACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGAt                                                 <  1:3870958/150‑1 (MQ=255)
                                                                                                      gACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGAt                                                 >  1:1993616/1‑150 (MQ=255)
                                                                                                       aCCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATc                                                <  2:3785258/150‑1 (MQ=255)
                                                                                                       aCCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATc                                                <  1:3388591/150‑1 (MQ=255)
                                                                                                       aCCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATc                                                <  1:3833688/150‑1 (MQ=255)
                                                                                                         cAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGATCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGTTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGa                                              >  2:3969951/1‑150 (MQ=255)
                                                                                                          aataatCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGAc                                             <  2:2636884/150‑1 (MQ=255)
                                                                                                             aatCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAAt                                          >  1:1205772/1‑150 (MQ=255)
                                                                                                                cGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCgg                                       >  1:7930/1‑150 (MQ=255)
                                                                                                                cGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCgg                                       <  2:3678492/150‑1 (MQ=255)
                                                                                                                  tCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCa                                     >  1:3009699/1‑150 (MQ=255)
                                                                                                                  tCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCa                                     >  1:3216966/1‑150 (MQ=255)
                                                                                                                   cAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCat                                    <  2:870958/150‑1 (MQ=255)
                                                                                                                      ttCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATc                                 <  1:152807/150‑1 (MQ=255)
                                                                                                                         tGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGtt                              >  1:1683701/1‑150 (MQ=255)
                                                                                                                              cgacgaCCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTtggtgg                         <  1:215470/150‑1 (MQ=255)
                                                                                                                                  gacCAGCAGTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACg                     <  1:501569/150‑1 (MQ=255)
                                                                                                                                   acCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGa                    <  2:3421086/150‑1 (MQ=255)
                                                                                                                                   acCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGa                    >  2:1170857/1‑150 (MQ=255)
                                                                                                                                    ccAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGaa                   <  2:4280682/150‑1 (MQ=255)
                                                                                                                                     cagcagTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAAt                  >  2:2512605/1‑150 (MQ=255)
                                                                                                                                      agcagTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATa                 >  1:2102735/1‑150 (MQ=255)
                                                                                                                                        cagTTCGCGAATAAGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGt               <  2:2057342/150‑1 (MQ=255)
                                                                                                                                          gTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTg             <  1:4445609/150‑1 (MQ=255)
                                                                                                                                             cgcgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgc          >  1:3728775/1‑150 (MQ=255)
                                                                                                                                             cgcgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgc          >  1:2711118/1‑150 (MQ=255)
                                                                                                                                             cgcgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgc          >  2:4102735/1‑150 (MQ=255)
                                                                                                                                             cgcgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgc          >  2:1851541/1‑150 (MQ=255)
                                                                                                                                              gcgAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgca         <  1:3964776/150‑1 (MQ=255)
                                                                                                                                                gAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGgcaag       >  2:2325307/1‑150 (MQ=255)
                                                                                                                                                    aCGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGcaagcaag   >  1:2953625/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGG  >  CP009273/4219401‑4219699

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: