Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,355,635 G→T 8.1% D87Y (GAC→TAC)  puuD → gamma‑glutamyl‑gamma‑aminobutyrate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,355,6350GT 8.1% 1228.9 / 49.2 331D87Y (GAC→TAC) puuDgamma‑glutamyl‑gamma‑aminobutyrate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (169/135);  new base T (8/19);  total (177/154)
Fisher's exact test for biased strand distribution p-value = 1.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 336 total aligned reads displayed.

GCAGGCGGCTTGCCTATTGCGCTACCACATGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCC  >  CP009273/1355488‑1355783
                                                                                                                                                   |                                                                                                                                                    
   ggcggcTTGCCTATTGCGCTACCACATGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGc                                                                                                                                                  >  1:705111/1‑149 (MQ=255)
      ggcTTGCCTATTGCGCTACCACTTGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCGGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGGGAAAACGGCGATGAGCCTGACGCCGa                                                                                                                                               <  2:3428597/149‑1 (MQ=255)
        cTTGCCTTTTGCCCTCCCACACGCGCTGGCGGAACCGTCATTAATTGAACAACTTGTGCCTAAACTCGAGGTCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTTAAAACGGCGATGAGCCTGACGCCGATc                                                                                                                                             <  1:3017750/149‑1 (MQ=255)
         ttGCCTATTGCGCTACCACATGCGCTGGCGGAACCGTCATTACTTGAACAACCTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTTACGCCGATcc                                                                                                                                            <  2:604942/149‑1 (MQ=255)
           gCCTATTGCGCTACCACATGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCg                                                                                                                                          >  1:4432207/1‑149 (MQ=255)
           gCCTATTGCGCTACCACATGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCg                                                                                                                                          >  1:1876079/1‑149 (MQ=255)
                  gcgcTACCACATGCGCTGGCGGACCCGTAATTACTTGAACAACTTTTGCCGAAACTCGTTGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGa                                                                                                                                   <  1:2255512/149‑1 (MQ=255)
                  gcgcTACCACATGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGa                                                                                                                                   <  2:3851013/149‑1 (MQ=255)
                     cTACCACATGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGAtct                                                                                                                                <  2:440654/149‑1 (MQ=255)
                         cacaTGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGAGGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATTTGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTg                                                                                                                            <  1:4780796/149‑1 (MQ=255)
                          acaTGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAACACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGa                                                                                                                           <  2:3328043/149‑1 (MQ=255)
                          acaTGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTTACGCCGATCCCGGGCGTGATCTTCTGa                                                                                                                           <  1:369708/149‑1 (MQ=255)
                                  tGGCGTAACCTTAATTACTTTAACAACTTTGCCCTAAACTCGTTGGAATTTATCGTCCTGGTAGTCCCAGAAATGTGCAGCCGCACCTATTTGGTTAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCg                                                                                                                   <  1:1613693/149‑1 (MQ=255)
                                  gggCGTAACCTTCATTACTTGAACAAGTTTTGCCGAAACTCGTTGGCATTTATCTCCGTGGTAGTCCCAGCATGGTGCCGCCGCACATATTTGGTCAACACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCg                                                                                                                   <  1:4000509/148‑1 (MQ=255)
                                   ggCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCACTGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGa                                                                                                                  >  2:91038/1‑149 (MQ=255)
                                   ggCGGAACCGTCATTACGTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGa                                                                                                                  >  2:1648267/1‑149 (MQ=255)
                                    gCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGAt                                                                                                                 <  2:3188671/149‑1 (MQ=255)
                                        aaCCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGAtaatt                                                                                                             <  2:2699983/149‑1 (MQ=255)
                                         aCCGTCATTACTTGAACAACTTTTGCCTAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGAtaatta                                                                                                            <  2:609526/149‑1 (MQ=255)
                                         aCCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGAtaatta                                                                                                            <  2:1565620/149‑1 (MQ=255)
                                          ccGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGAtaattaa                                                                                                           >  2:4687285/1‑149 (MQ=255)
                                           cGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGAtaattaat                                                                                                          >  1:2738651/1‑149 (MQ=255)
                                            gTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATg                                                                                                         <  2:225055/149‑1 (MQ=255)
                                             tCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGc                                                                                                        >  2:4111171/1‑149 (MQ=255)
                                             tCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGc                                                                                                        >  1:774653/1‑149 (MQ=255)
                                             tCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGc                                                                                                        >  1:2021142/1‑149 (MQ=255)
                                                ttACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGc                                                                                                     >  2:4136063/1‑149 (MQ=255)
                                                  aCTTTAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACATATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCAc                                                                                                   <  1:596094/149‑1 (MQ=255)
                                                   cTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACt                                                                                                  >  2:2433376/1‑149 (MQ=255)
                                                        acaacTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTTACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGaaa                                                                                             <  1:2318633/149‑1 (MQ=255)
                                                        acaacTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATAGGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGaaa                                                                                             <  2:1073603/149‑1 (MQ=255)
                                                          aacTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTTACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAgg                                                                                           >  1:2810042/1‑149 (MQ=255)
                                                           acTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGgc                                                                                          <  1:4187109/149‑1 (MQ=255)
                                                                  cctaaaCTCTATGGCATTTATCTTCCGGTTAGTCCCAGCAATGTGCAGCCGCACATTTATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATccc                                                                                   <  1:950781/146‑1 (MQ=255)
                                                                  cctaaaCTCGTTGGCATTTATCTTCGTGGTTGTCCAAGAAATGTGCAGCCGAACATATATGGTTAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATccc                                                                                   <  1:2399195/146‑1 (MQ=255)
                                                                  ccGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATccc                                                                                   <  2:263164/149‑1 (MQ=255)
                                                                   cGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATcccc                                                                                  >  2:1582861/1‑149 (MQ=255)
                                                                     aaaCTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCAt                                                                                >  2:4703550/1‑149 (MQ=255)
                                                                     aaaCTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCAt                                                                                <  1:1142379/149‑1 (MQ=255)
                                                                        cTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCAtttt                                                                             >  1:2290311/1‑149 (MQ=255)
                                                                        cTCGATGGAATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACATATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCAtttt                                                                             <  2:4643837/149‑1 (MQ=255)
                                                                         tCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTc                                                                            >  2:931806/1‑149 (MQ=255)
                                                                           gttGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATTTGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGc                                                                          <  1:1477844/147‑1 (MQ=255)
                                                                             tGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCa                                                                        <  2:739814/149‑1 (MQ=255)
                                                                              ggCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTTACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCAt                                                                       <  2:2289424/149‑1 (MQ=255)
                                                                              ggCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCAt                                                                       >  1:1038173/1‑149 (MQ=255)
                                                                              ggCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCAt                                                                       >  2:3537710/1‑149 (MQ=255)
                                                                              ggCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCAt                                                                       >  2:1697850/1‑149 (MQ=255)
                                                                                cATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCt                                                                     <  2:3231102/149‑1 (MQ=255)
                                                                                  tttATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGc                                                                   >  2:1453281/1‑149 (MQ=255)
                                                                                  tttATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGc                                                                   >  1:182505/1‑149 (MQ=255)
                                                                                      tCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCggg                                                               >  1:4412862/1‑149 (MQ=255)
                                                                                      tCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCggg                                                               >  2:1047780/1‑149 (MQ=255)
                                                                                        ttCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGt                                                             >  1:2585986/1‑149 (MQ=255)
                                                                                         tCCTGGTAGTCCCAGCAATGTGCTGCCGCACCTTTTTGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGtt                                                            <  1:4703181/149‑1 (MQ=255)
                                                                                           cTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTa                                                          <  1:4370232/149‑1 (MQ=255)
                                                                                           cTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTa                                                          >  2:1891090/1‑149 (MQ=255)
                                                                                            tGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTAc                                                         >  1:1309137/1‑149 (MQ=255)
                                                                                             ggTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACa                                                        >  1:2907377/1‑149 (MQ=255)
                                                                                               tAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAg                                                      >  2:1869564/1‑149 (MQ=255)
                                                                                                aGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGa                                                     >  2:2910769/1‑149 (MQ=255)
                                                                                                 gTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGaa                                                    >  2:1992448/1‑149 (MQ=255)
                                                                                                 gTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGaa                                                    <  2:3903953/149‑1 (MQ=255)
                                                                                                 gTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGaa                                                    >  1:3510944/1‑149 (MQ=255)
                                                                                                 gTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGaa                                                    >  2:2870134/1‑149 (MQ=255)
                                                                                                   cccAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACt                                                  >  1:1532484/1‑149 (MQ=255)
                                                                                                   cccAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACt                                                  >  2:842829/1‑149 (MQ=255)
                                                                                                   cccAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACt                                                  <  1:2772060/149‑1 (MQ=255)
                                                                                                         aaTGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTTACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACtggtggt                                            <  2:3628483/149‑1 (MQ=255)
                                                                                                         aaTGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACtggtggt                                            >  1:3137295/1‑149 (MQ=255)
                                                                                                          aTGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACtggtggtg                                           >  2:1369993/1‑149 (MQ=255)
                                                                                                            gtgCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGc                                         <  2:311055/149‑1 (MQ=255)
                                                                                                               cAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAAc                                      <  2:1762461/149‑1 (MQ=255)
                                                                                                                 gCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCg                                    >  1:1886754/1‑149 (MQ=255)
                                                                                                                 gCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCg                                    >  1:632112/1‑149 (MQ=255)
                                                                                                                 gCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCg                                    >  1:4488270/1‑149 (MQ=255)
                                                                                                                        cTATATGGTGAAAACGGCGATGAGCCTTACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTc                             <  2:4055822/149‑1 (MQ=255)
                                                                                                                        cTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTc                             >  2:4505630/1‑149 (MQ=255)
                                                                                                                        cTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTc                             <  1:3323242/149‑1 (MQ=255)
                                                                                                                         tataTGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCg                            >  2:1832013/1‑149 (MQ=255)
                                                                                                                          atatGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGt                           >  2:4527543/1‑149 (MQ=255)
                                                                                                                          atatGGGGAACACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGt                           <  1:1375906/149‑1 (MQ=255)
                                                                                                                             tGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGc                        <  1:4114833/149‑1 (MQ=255)
                                                                                                                             tGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGc                        >  2:4360233/1‑149 (MQ=255)
                                                                                                                                 gAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCg                    >  1:1711417/1‑149 (MQ=255)
                                                                                                                                  aaaaCGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGc                   >  1:3100303/1‑149 (MQ=255)
                                                                                                                                  aaaaCGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGc                   >  1:895765/1‑149 (MQ=255)
                                                                                                                                  aaaaCGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGc                   >  1:3798991/1‑149 (MQ=255)
                                                                                                                                   aaaCGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCa                  >  1:3975136/1‑149 (MQ=255)
                                                                                                                                      cGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGc               >  2:3026253/1‑149 (MQ=255)
                                                                                                                                        gCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCtg             >  2:877480/1‑149 (MQ=255)
                                                                                                                                         cGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCtgt            <  2:356792/149‑1 (MQ=255)
                                                                                                                                          gATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCtgtg           <  2:301900/149‑1 (MQ=255)
                                                                                                                                          gATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCtgtg           >  2:814643/1‑149 (MQ=255)
                                                                                                                                           aTGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGc          >  2:409796/1‑149 (MQ=255)
                                                                                                                                           aTGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGc          >  2:1270766/1‑149 (MQ=255)
                                                                                                                                            tGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCg         >  2:2407150/1‑149 (MQ=255)
                                                                                                                                            tGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCg         >  1:865412/1‑149 (MQ=255)
                                                                                                                                               gCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAAc      >  1:3591757/1‑149 (MQ=255)
                                                                                                                                                   gACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGcc  <  2:4788681/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAGGCGGCTTGCCTATTGCGCTACCACATGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCC  >  CP009273/1355488‑1355783

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