Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,790 A→T 69.9% L94Q (CTG→CAG)  pepE ← (alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,219,7900AT69.9% 469.4 / 279.2 316L94Q (CTG→CAG) pepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (43/52);  new base T (103/118);  total (146/170)
Fisher's exact test for biased strand distribution p-value = 9.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 316 total aligned reads displayed.

CTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGG  >  CP009273/4219644‑4219937
                                                                                                                                                  |                                                                                                                                                   
   gCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACt                                                                                                                                               <  1:2562209/150‑1 (MQ=255)
     ggATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTgg                                                                                                                                             >  2:1947006/1‑150 (MQ=255)
      gATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGa                                                                                                                                            <  1:2408636/150‑1 (MQ=255)
       aTCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGaa                                                                                                                                           <  1:3228767/150‑1 (MQ=255)
        tCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGaaa                                                                                                                                          >  2:1140765/1‑150 (MQ=255)
        tCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCCGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGaaa                                                                                                                                          <  1:352421/150‑1 (MQ=255)
         cGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCACGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCACCTGGAAAg                                                                                                                                         >  2:3678245/1‑150 (MQ=255)
          gACAATCGGCATCTCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGt                                                                                                                                        <  1:990792/150‑1 (MQ=255)
          gACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGt                                                                                                                                        >  2:4402734/1‑150 (MQ=255)
           aCAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTa                                                                                                                                       <  1:3719463/150‑1 (MQ=255)
           aCAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTa                                                                                                                                       >  1:3037213/1‑150 (MQ=255)
           aCAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTa                                                                                                                                       <  1:1504286/150‑1 (MQ=255)
           aCAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACATCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTa                                                                                                                                       <  1:3686095/150‑1 (MQ=255)
           aCAATCGGCATATCGTTGGTGGTACGAATAGGTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAAGTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTa                                                                                                                                       >  2:3969080/1‑150 (MQ=255)
             aaTCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTAtt                                                                                                                                     <  2:1488446/150‑1 (MQ=255)
             aaTCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTAtt                                                                                                                                     >  2:380377/1‑150 (MQ=255)
                cGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTccc                                                                                                                                  >  1:4637895/1‑150 (MQ=255)
                cGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTccc                                                                                                                                  >  2:117042/1‑150 (MQ=255)
                 ggCATATCGTGGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCccg                                                                                                                                 <  2:4198255/150‑1 (MQ=255)
                  gCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCccgc                                                                                                                                >  1:1193980/1‑150 (MQ=255)
                   cATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCccgcc                                                                                                                               <  1:263009/150‑1 (MQ=255)
                    atatCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCccgccg                                                                                                                              >  1:704014/1‑150 (MQ=255)
                         gTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGAc                                                                                                                         <  2:1242307/150‑1 (MQ=255)
                           tggtggTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACaa                                                                                                                       >  1:3582299/1‑150 (MQ=255)
                             gtggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATc                                                                                                                     <  1:1526693/150‑1 (MQ=255)
                               ggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCAc                                                                                                                   <  2:3905115/150‑1 (MQ=255)
                                    gAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATCGACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTACCGCCGCCGACAATCACGAtct                                                                                                              <  2:2688040/150‑1 (MQ=255)
                                    gAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTACTTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGAtct                                                                                                              >  1:2884246/1‑150 (MQ=255)
                                    gAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGAtct                                                                                                              >  1:1278626/1‑150 (MQ=255)
                                    gAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGAtct                                                                                                              <  2:4600377/150‑1 (MQ=255)
                                     aaTAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGAtctc                                                                                                             >  2:1796976/1‑150 (MQ=255)
                                       tAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAg                                                                                                           <  1:4273747/150‑1 (MQ=255)
                                        aGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGc                                                                                                          <  1:4078688/150‑1 (MQ=255)
                                         gTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCa                                                                                                         <  2:2758353/150‑1 (MQ=255)
                                          ttGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAt                                                                                                        <  2:3240686/150‑1 (MQ=255)
                                           tGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAtt                                                                                                       <  1:3837628/150‑1 (MQ=255)
                                               caagcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCa                                                                                                   >  1:3400400/1‑150 (MQ=255)
                                               caagcaagGTTAGCGCCTGCGCTCCAGCCAATAGACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCa                                                                                                   <  2:655056/150‑1 (MQ=255)
                                                 agcaagGTTAGCGCCTGCGCTCCAGCCAATCTACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAAt                                                                                                 <  1:1895816/150‑1 (MQ=255)
                                                 agcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAAt                                                                                                 >  1:4118327/1‑150 (MQ=255)
                                                  gcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATc                                                                                                >  1:1826476/1‑150 (MQ=255)
                                                  gcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATc                                                                                                <  2:4644746/150‑1 (MQ=255)
                                                  gcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATc                                                                                                <  2:2997443/150‑1 (MQ=255)
                                                    aagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgc                                                                                              >  1:2183050/1‑150 (MQ=255)
                                                    aagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgc                                                                                              >  1:1032282/1‑150 (MQ=255)
                                                     agGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgcg                                                                                             <  1:599843/150‑1 (MQ=255)
                                                       gTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGc                                                                                           >  2:1447088/1‑150 (MQ=255)
                                                       gTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGc                                                                                           >  2:311996/1‑150 (MQ=255)
                                                           gcgcCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACg                                                                                       <  1:5056782/150‑1 (MQ=255)
                                                           gcgcCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACg                                                                                       <  2:506384/150‑1 (MQ=255)
                                                               cTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGAt                                                                                   <  2:2571213/150‑1 (MQ=255)
                                                                tGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATc                                                                                  <  2:2892692/150‑1 (MQ=255)
                                                                 gcgcTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATct                                                                                 <  1:4839130/150‑2 (MQ=255)
                                                                  cgcTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATcta                                                                                <  1:1022123/150‑3 (MQ=255)
                                                                    cTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACa                                                                              <  1:2005899/150‑1 (MQ=255)
                                                                      ccAGCCACTATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAAc                                                                            <  2:3564024/150‑1 (MQ=255)
                                                                      ccAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCACTCGCGGCAACGGGATCGACAAc                                                                            >  2:1688333/1‑150 (MQ=255)
                                                                      ccAGCCAATAGACAGAGCGCCACGTTGAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAAc                                                                            <  2:577978/150‑1 (MQ=255)
                                                                        aGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCACCGGGATCTACAACGc                                                                          >  2:3520759/1‑150 (MQ=255)
                                                                        aGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGc                                                                          >  1:4543905/1‑150 (MQ=255)
                                                                             atatACAGAGCGCCACGTTTAACCACGGCAGTAATTGGTGCCAGCGGCCCGCGCTCGCGGCTCTGTTTCTGCAACTCGAAAGTCTTCCCGCCGCCGACAGTCACGGTCTCAGCCATTTCTCTCGCGGCACCGGGGTCTACAGCGCTATGa                                                                     >  2:4983470/1‑150 (MQ=255)
                                                                               ataCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAAt                                                                   >  2:4476843/1‑150 (MQ=255)
                                                                                taCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATa                                                                  <  2:3060401/150‑1 (MQ=255)
                                                                                 aCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATAc                                                                 <  2:4620104/150‑1 (MQ=255)
                                                                                    gagCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCgg                                                              >  2:228020/1‑150 (MQ=255)
                                                                                     agCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGt                                                             <  1:3350278/150‑1 (MQ=255)
                                                                                       cgcCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGa                                                           <  2:473959/150‑1 (MQ=255)
                                                                                         ccACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACa                                                         <  2:2121600/150‑1 (MQ=255)
                                                                                           aCGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGa                                                       <  1:2855423/150‑1 (MQ=255)
                                                                                                tAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAacac                                                  >  1:3996038/1‑150 (MQ=255)
                                                                                                  aCCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACAccc                                                >  2:1987968/1‑150 (MQ=255)
                                                                                                   ccACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCa                                               <  1:3960937/150‑1 (MQ=255)
                                                                                                   ccACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCa                                               <  1:2989696/150‑1 (MQ=255)
                                                                                                    cACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAg                                              <  2:3419704/150‑1 (MQ=255)
                                                                                                     aCGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGc                                             <  2:794746/150‑1 (MQ=255)
                                                                                                            tAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCga                                      >  1:1780826/1‑150 (MQ=255)
                                                                                                                tGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGGAACACCCAGCGGAGCGAGAAc                                  <  2:1856877/150‑1 (MQ=255)
                                                                                                                 ggTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAAcc                                 <  2:2444071/150‑1 (MQ=255)
                                                                                                                       agcagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                           >  1:4994536/1‑150 (MQ=255)
                                                                                                                          agcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                        >  2:333516/1‑150 (MQ=255)
                                                                                                                          agcCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                        >  1:415916/1‑150 (MQ=255)
                                                                                                                           gcCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGttt                       <  2:4302417/150‑1 (MQ=255)
                                                                                                                           gcCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTACAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGttt                       <  2:146591/150‑1 (MQ=255)
                                                                                                                            cccGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtttt                      <  1:5216375/150‑1 (MQ=255)
                                                                                                                             ccGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTc                     >  2:550608/1‑150 (MQ=255)
                                                                                                                             ccGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTc                     <  1:2239196/150‑1 (MQ=255)
                                                                                                                             ccGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTc                     <  1:3451627/150‑1 (MQ=255)
                                                                                                                                   tCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCtgtg               <  2:2206214/150‑1 (MQ=255)
                                                                                                                                     gcgGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTa             >  1:2799905/1‑150 (MQ=255)
                                                                                                                                         cACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatca         >  2:1283135/1‑150 (MQ=255)
                                                                                                                                         cACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatca         <  1:3147246/150‑1 (MQ=255)
                                                                                                                                          aCTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTTatcat        <  2:1032363/150‑1 (MQ=255)
                                                                                                                                          aCTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatcat        <  1:1715158/150‑1 (MQ=255)
                                                                                                                                            tGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc      <  2:4824578/150‑1 (MQ=255)
                                                                                                                                            tGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc      <  2:4269202/150‑1 (MQ=255)
                                                                                                                                            tGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc      <  2:704014/150‑1 (MQ=255)
                                                                                                                                             gTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATccc     >  2:2068642/1‑150 (MQ=255)
                                                                                                                                             gTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATccc     >  2:4541990/1‑150 (MQ=255)
                                                                                                                                              tttCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCa    <  2:1571850/150‑1 (MQ=255)
                                                                                                                                                tCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAgg  <  1:4402734/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGG  >  CP009273/4219644‑4219937

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: