Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,225,529 G→T 6.2% S69S (TCG→TCT cdd → cytidine/deoxycytidine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,225,5290GT 6.2% 1210.5 / 35.2 321S69S (TCG→TCTcddcytidine/deoxycytidine deaminase
Reads supporting (aligned to +/- strand):  ref base G (145/156);  new base T (4/16);  total (149/172)
Fisher's exact test for biased strand distribution p-value = 1.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 322 total aligned reads displayed.

GGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCT  >  CP009273/2225382‑2225678
                                                                                                                                                   |                                                                                                                                                     
  aaCCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAAtt                                                                                                                                                   <  1:4185023/150‑1 (MQ=255)
    ccTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAAtttt                                                                                                                                                 >  1:2841510/1‑150 (MQ=255)
     cTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGATGGCATTCGCACTACTTCCTCTGGCGGCGGCCTGTGCGCGTACTCCATTGTCGAATTTTa                                                                                                                                                >  2:1007873/1‑150 (MQ=255)
        ttCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATg                                                                                                                                             <  1:364167/150‑1 (MQ=255)
        ttCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGAAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATg                                                                                                                                             >  2:3318644/1‑150 (MQ=255)
              cAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGgcg                                                                                                                                       >  2:2960554/1‑150 (MQ=255)
              cAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGgcg                                                                                                                                       <  2:4009812/150‑1 (MQ=255)
              cAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGgcg                                                                                                                                       >  2:2839745/1‑150 (MQ=255)
                   aaGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAAtt                                                                                                                                  <  2:729825/150‑1 (MQ=255)
                    aGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTg                                                                                                                                 >  1:1775710/1‑150 (MQ=255)
                     gTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgc                                                                                                                                >  2:1454704/1‑150 (MQ=255)
                     gTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgc                                                                                                                                >  1:1411352/1‑150 (MQ=255)
                         ttCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcgcgc                                                                                                                            >  2:906330/1‑150 (MQ=255)
                         ttCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcgcgc                                                                                                                            >  1:3387873/1‑150 (MQ=255)
                         ttCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcgcgc                                                                                                                            >  1:2990433/1‑150 (MQ=255)
                              cGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCACTTGCGCGCGgtgt                                                                                                                       >  2:1798351/1‑150 (MQ=255)
                                 tttGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCACTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAg                                                                                                                    <  2:4473187/150‑1 (MQ=255)
                                  ttgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGc                                                                                                                   <  1:2276137/150‑1 (MQ=255)
                                   tgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCg                                                                                                                  <  1:2729640/150‑1 (MQ=255)
                                   tgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCg                                                                                                                  <  1:2833214/150‑1 (MQ=255)
                                   tgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCg                                                                                                                  <  1:1246936/150‑1 (MQ=255)
                                      tgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGaa                                                                                                               >  2:2237558/1‑150 (MQ=255)
                                       gACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAAc                                                                                                              >  1:1791539/1‑150 (MQ=255)
                                       gACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCAGTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAAc                                                                                                              >  2:554644/1‑150 (MQ=255)
                                         ccGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCt                                                                                                            >  2:1162303/1‑150 (MQ=255)
                                         ccGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCt                                                                                                            >  2:1162189/1‑150 (MQ=255)
                                         ccGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCt                                                                                                            >  2:2372185/1‑150 (MQ=255)
                                         ccGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCt                                                                                                            <  2:2505570/150‑1 (MQ=255)
                                          cGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTg                                                                                                           <  1:4299939/150‑1 (MQ=255)
                                             ggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTa                                                                                                        <  1:3144431/150‑1 (MQ=255)
                                              gAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTAt                                                                                                       <  1:3516391/150‑1 (MQ=255)
                                                gCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTAttt                                                                                                     >  2:2315392/1‑150 (MQ=255)
                                                   aGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCgg                                                                                                  >  2:1170355/1‑150 (MQ=255)
                                                      ctcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                                                                                               >  2:4383199/1‑150 (MQ=255)
                                                      ctcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                                                                                               <  1:2904410/150‑1 (MQ=255)
                                                         aTCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCaa                                                                                            <  2:2924413/150‑1 (MQ=255)
                                                          tCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAat                                                                                           <  1:1212040/150‑1 (MQ=255)
                                                            gCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCTAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAatat                                                                                         <  2:3762484/150‑1 (MQ=255)
                                                             cTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCTAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATg                                                                                        <  1:1427204/150‑1 (MQ=255)
                                                              tGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATgg                                                                                       <  1:594402/150‑1 (MQ=255)
                                                              tGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATgg                                                                                       <  1:267336/150‑1 (MQ=255)
                                                              tGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATgg                                                                                       <  1:44540/150‑1 (MQ=255)
                                                               gAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGa                                                                                      >  2:3911332/1‑150 (MQ=255)
                                                                  gagCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCTAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAAtt                                                                                   <  2:1163854/150‑1 (MQ=255)
                                                                  gagCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAAtt                                                                                   <  1:1487414/150‑1 (MQ=255)
                                                                   agCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAAttt                                                                                  >  2:166973/1‑150 (MQ=255)
                                                                    gcgcAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTa                                                                                 <  2:2006720/150‑1 (MQ=255)
                                                                    gcgcAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTa                                                                                 <  2:3409840/150‑1 (MQ=255)
                                                                      gcAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAtt                                                                               <  2:1882695/150‑1 (MQ=255)
                                                                       cAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTg                                                                              >  1:2092064/1‑150 (MQ=255)
                                                                       cAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTg                                                                              >  2:779326/1‑150 (MQ=255)
                                                                       cAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTg                                                                              <  1:3318609/150‑1 (MQ=255)
                                                                           ggggCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGc                                                                          >  2:3235090/1‑150 (MQ=255)
                                                                            gggCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCTGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCg                                                                         <  2:4292377/150‑1 (MQ=255)
                                                                             ggCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGa                                                                        <  2:4029569/150‑1 (MQ=255)
                                                                               cTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACa                                                                      <  1:1929324/150‑1 (MQ=255)
                                                                                tGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACaa                                                                     >  1:4310414/1‑150 (MQ=255)
                                                                                tGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACaa                                                                     >  1:1335966/1‑150 (MQ=255)
                                                                                 ggACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCTAATTTTAATGTTGGCGCAACTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAAt                                                                    <  1:935875/150‑1 (MQ=255)
                                                                                 ggACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAAt                                                                    <  1:2006517/150‑1 (MQ=255)
                                                                                     gAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATgcag                                                                >  1:3328410/1‑150 (MQ=255)
                                                                                       aGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCGGGTATTTCGGTGCCAATATGGAATTTAGTGGTGCGACACTgcagca                                                              >  2:2550470/1‑150 (MQ=255)
                                                                                           gcgcTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAAcc                                                          <  1:2977398/150‑1 (MQ=255)
                                                                                            cgcTGGGATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCTAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCg                                                         <  1:4407304/150‑1 (MQ=255)
                                                                                              cTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGtt                                                       <  1:1140110/150‑1 (MQ=255)
                                                                                               tGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTc                                                      <  2:396304/150‑1 (MQ=255)
                                                                                               tGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTc                                                      <  2:3893055/150‑1 (MQ=255)
                                                                                               tGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTc                                                      <  2:2691414/150‑1 (MQ=255)
                                                                                                ggCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCa                                                     >  2:4506738/1‑150 (MQ=255)
                                                                                                   aTTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGc                                                  >  2:2272686/1‑150 (MQ=255)
                                                                                                    ttCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGcc                                                 >  2:44356/1‑150 (MQ=255)
                                                                                                      cGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGa                                               <  1:1759541/150‑1 (MQ=255)
                                                                                                      cGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGa                                               <  1:4013249/150‑1 (MQ=255)
                                                                                                      cGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAAGTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGa                                               <  2:4320710/150‑1 (MQ=255)
                                                                                                         actactTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACa                                            <  2:1363455/150‑1 (MQ=255)
                                                                                                         actactTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACa                                            >  2:2142647/1‑150 (MQ=255)
                                                                                                           tactTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTCTTTCGGTGCCACTATGGAATTTATGGGTGCGACAATGCAGCAAACCGTTCATGCCGAACaaa                                          >  2:4051702/1‑150 (MQ=255)
                                                                                                            actTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAg                                         >  2:3561111/1‑150 (MQ=255)
                                                                                                            actTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAg                                         >  2:4246140/1‑150 (MQ=255)
                                                                                                                  gCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGAt                                   <  2:805847/150‑1 (MQ=255)
                                                                                                                   cTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATc                                  <  1:1335235/150‑1 (MQ=255)
                                                                                                                    tGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCa                                 >  1:2425613/1‑150 (MQ=255)
                                                                                                                    tGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGGGTGCGCGGAACCTGGTATTTCGGTGCCAATCTGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACACAGCGCGATCa                                 >  2:3827444/1‑150 (MQ=255)
                                                                                                                    tGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAAGTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCa                                 >  1:1797100/1‑150 (MQ=255)
                                                                                                                     ggcggcggcCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAg                                >  2:1645620/1‑150 (MQ=255)
                                                                                                                     ggcggcggcCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAg                                >  2:1646098/1‑150 (MQ=255)
                                                                                                                          cggcCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACg                           <  2:3346319/150‑1 (MQ=255)
                                                                                                                          cggcCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACg                           >  2:64636/1‑150 (MQ=255)
                                                                                                                          cggcCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACg                           >  2:3137518/1‑150 (MQ=255)
                                                                                                                                  gcgcGTACGCCCTTGTCGAATTTTAATGTTGGCGCAAATGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTg                   >  2:3665022/1‑150 (MQ=255)
                                                                                                                                   cgcgTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGa                  >  1:1027677/1‑150 (MQ=255)
                                                                                                                                      gTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAgtg               >  1:806033/1‑150 (MQ=255)
                                                                                                                                         cGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAgtggtg            >  2:546945/1‑150 (MQ=255)
                                                                                                                                          gCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGa           >  1:555800/1‑150 (MQ=255)
                                                                                                                                          gCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGa           >  1:1577630/1‑150 (MQ=255)
                                                                                                                                           ccATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGaa          >  2:4193498/1‑150 (MQ=255)
                                                                                                                                             aTTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGaaaa        >  1:3920833/1‑150 (MQ=255)
                                                                                                                                              ttGTCGAATTTTAATGTTGGCGCAATTGCGCGCTGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGATCATGCCGAACAAAGCGCGATCAGACACGCCTGGTTGAGTGGTGaaaaa       >  2:3884628/1‑150 (MQ=255)
                                                                                                                                                gTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAgc     >  2:1060201/1‑150 (MQ=255)
                                                                                                                                                 tCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAgcg    >  1:4508541/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCT  >  CP009273/2225382‑2225678

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: