Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,360,489 G→T 12.0% C228F (TGC→TTC)  puuE → 4‑aminobutyrate aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,360,4890GT12.0% 755.3 / 54.2 222C228F (TGC→TTC) puuE4‑aminobutyrate aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (119/76);  new base T (15/12);  total (134/88)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

AACACAGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGC  >  CP009273/1360340‑1360637
                                                                                                                                                     |                                                                                                                                                    
aaCACAGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCtgtg                                                                                                                                                      <  2:244428/150‑1 (MQ=255)
 acacATGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGACGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGGTGCCGCTATTCGCCGCCTGTGc                                                                                                                                                     >  2:1655768/1‑150 (MQ=255)
  cacaGGACTCCCTCGACGCCATCGAACGCTTGTTTATATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGcg                                                                                                                                                    <  1:2537567/150‑1 (MQ=255)
   acaGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGcga                                                                                                                                                   >  1:1670804/1‑150 (MQ=255)
    caGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGcgac                                                                                                                                                  <  2:291342/150‑1 (MQ=255)
     aGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTTcgacg                                                                                                                                                 >  2:1467028/1‑150 (MQ=255)
     aGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGAGTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCACAAGAGCTGGTTGCCGCTATTCGCCGCCTGTTcgacg                                                                                                                                                 >  2:3648150/1‑150 (MQ=255)
     aGGACTCCCTCGACGCCATCGAACGCATGTTTAAATCAGACAACGAAGCGAAGCAGGTGGCGGCGATTATATTCGAACCGGAGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCAGCCTGTTcgacg                                                                                                                                                 >  2:1467541/1‑150 (MQ=255)
       gACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTTCGACGAg                                                                                                                                               <  1:2689226/150‑1 (MQ=255)
           cccTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACg                                                                                                                                           >  2:2027530/1‑150 (MQ=255)
            ccTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGAGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCAGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACgg                                                                                                                                          >  1:3281998/1‑150 (MQ=255)
             cTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGt                                                                                                                                         >  1:2395361/1‑150 (MQ=255)
             cTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGt                                                                                                                                         >  1:2395770/1‑150 (MQ=255)
               cGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACTGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTAt                                                                                                                                       >  1:2409958/1‑150 (MQ=255)
                gACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCACCGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTAtt                                                                                                                                      <  1:806196/150‑1 (MQ=255)
                gACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTAtt                                                                                                                                      <  1:806275/150‑1 (MQ=255)
                    ccATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGtga                                                                                                                                  >  1:284896/1‑150 (MQ=255)
                     cATCGAACGCTTGTTTACATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGtgat                                                                                                                                 <  2:2952530/150‑1 (MQ=255)
                      aTCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGtgatg                                                                                                                                >  1:3183303/1‑150 (MQ=255)
                         gAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGAtt                                                                                                                             <  1:3471860/150‑1 (MQ=255)
                         gAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGAtt                                                                                                                             >  2:437020/1‑150 (MQ=255)
                         gAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGAtt                                                                                                                             <  1:442017/150‑1 (MQ=255)
                         gAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGAtt                                                                                                                             <  1:442006/150‑1 (MQ=255)
                            cGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCt                                                                                                                          >  1:3874056/1‑150 (MQ=255)
                            cGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCt                                                                                                                          <  2:3959045/150‑1 (MQ=255)
                            cGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCt                                                                                                                          >  2:4426276/1‑150 (MQ=255)
                               ttGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCtgat                                                                                                                       >  2:2352352/1‑150 (MQ=255)
                                tGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCtgatg                                                                                                                      >  2:3903033/1‑150 (MQ=255)
                                  tttAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGaa                                                                                                                    >  1:1090831/1‑150 (MQ=255)
                                    tAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGt                                                                                                                  >  2:389672/1‑150 (MQ=255)
                                     aaaTCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTg                                                                                                                 >  1:2418549/1‑150 (MQ=255)
                                       aTCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCa                                                                                                               <  2:2635278/150‑1 (MQ=255)
                                         cAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCaaa                                                                                                             >  1:4204804/1‑150 (MQ=255)
                                            aCATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTTCGACGAGCACGGTATTGTGATGATTGCTGATGTAGTGCAAAGCg                                                                                                          >  2:343878/1‑150 (MQ=255)
                                            aCATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCg                                                                                                          >  1:2695963/1‑150 (MQ=255)
                                              aTCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTTCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGc                                                                                                        <  1:2821947/150‑1 (MQ=255)
                                                 gaagcgaagcAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCttt                                                                                                     >  2:4486876/1‑150 (MQ=255)
                                                 gaagcgaagcAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGAGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCttt                                                                                                     >  2:3082618/1‑150 (MQ=255)
                                                 gaagcgaagcAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGAGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCttt                                                                                                     >  2:3083195/1‑150 (MQ=255)
                                                           aGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCg                                                                                           <  1:3900668/150‑1 (MQ=255)
                                                              tGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAATAGCTGGTTGCCGCTATTCGCCGCCTGTTCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTa                                                                                        >  2:105315/1‑150 (MQ=255)
                                                               ggcggcGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTaa                                                                                       >  2:3828766/1‑150 (MQ=255)
                                                                gcggcgATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAg                                                                                      <  2:61231/150‑1 (MQ=255)
                                                                 cggcgATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGc                                                                                     >  1:3037639/1‑150 (MQ=255)
                                                                  ggcgTTTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAATAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCt                                                                                    <  2:839549/150‑1 (MQ=255)
                                                                  ggcgATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCt                                                                                    <  1:1817327/150‑1 (MQ=255)
                                                                     gATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGtt                                                                                 >  2:1437670/1‑150 (MQ=255)
                                                                        tattTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTTCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGc                                                                              >  1:4260373/1‑150 (MQ=255)
                                                                         attTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGcc                                                                             >  1:180005/1‑150 (MQ=255)
                                                                          ttttCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCa                                                                            <  1:1044748/150‑1 (MQ=255)
                                                                            ttCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATg                                                                          <  2:3940163/150‑1 (MQ=255)
                                                                            ttCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATg                                                                          <  2:3940150/150‑1 (MQ=255)
                                                                             tCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATgg                                                                         >  2:811953/1‑150 (MQ=255)
                                                                              cGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTTCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGa                                                                        <  2:1207153/150‑1 (MQ=255)
                                                                              cGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGa                                                                        >  2:750096/1‑150 (MQ=255)
                                                                               gAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGAt                                                                       >  1:40147/1‑150 (MQ=255)
                                                                               gAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGAt                                                                       >  1:1247855/1‑150 (MQ=255)
                                                                               gAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGAt                                                                       >  2:2272757/1‑150 (MQ=255)
                                                                                aaCCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCAGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATc                                                                      <  1:4317054/150‑1 (MQ=255)
                                                                                aaCCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCAGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATc                                                                      <  1:1533356/150‑1 (MQ=255)
                                                                                  ccGGTGCTGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCAt                                                                    <  2:1677144/150‑1 (MQ=255)
                                                                                   cGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCAtt                                                                   >  2:603670/1‑150 (MQ=255)
                                                                                    ggTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTa                                                                  >  2:1737554/1‑150 (MQ=255)
                                                                                        cAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTTCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGcc                                                              <  1:4492254/150‑1 (MQ=255)
                                                                                        cAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGcc                                                              <  1:2027530/150‑1 (MQ=255)
                                                                                            gcgagggcgGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATa                                                          >  1:1531743/1‑150 (MQ=255)
                                                                                              gagggcgGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAg                                                        >  2:2750871/1‑150 (MQ=255)
                                                                                                gggcgGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGcc                                                      <  1:3664389/150‑1 (MQ=255)
                                                                                                 ggcgGTTTCAACCTCGCGCCAAACGAGCTGGTTGCCGCTATTCGCCGCCTGTTCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCg                                                     <  2:4205878/150‑1 (MQ=255)
                                                                                                   cgGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTTCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGa                                                   <  1:3270619/150‑1 (MQ=255)
                                                                                                    ggTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGAt                                                  >  1:2920640/1‑150 (MQ=255)
                                                                                                    ggTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGAt                                                  >  2:579926/1‑150 (MQ=255)
                                                                                                     gTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGAtt                                                 >  2:2831025/1‑150 (MQ=255)
                                                                                                     gTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGAtt                                                 <  2:2254389/150‑1 (MQ=255)
                                                                                                     gTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGAtt                                                 >  2:3256120/1‑150 (MQ=255)
                                                                                                       ttCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTa                                               >  1:482455/1‑150 (MQ=255)
                                                                                                            cGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGAc                                          <  2:3661003/150‑1 (MQ=255)
                                                                                                            cGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGAc                                          >  1:1482555/1‑150 (MQ=255)
                                                                                                            cGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGAc                                          >  1:2910291/1‑150 (MQ=255)
                                                                                                            cGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGAc                                          >  2:1817962/1‑150 (MQ=255)
                                                                                                            cGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGAATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCCTGGATCATTACGCCGATAAGCCGGATTTAATGAc                                          >  2:1817943/1‑150 (MQ=255)
                                                                                                               tGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCAAGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGAt                                       <  2:230383/150‑1 (MQ=255)
                                                                                                                    cAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGa                                  >  2:192926/1‑150 (MQ=255)
                                                                                                                         gagCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGc                             >  1:3328676/1‑150 (MQ=255)
                                                                                                                           gCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCt                           <  1:452317/150‑1 (MQ=255)
                                                                                                                               gTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTcgcg                       >  1:4353306/1‑150 (MQ=255)
                                                                                                                                      cTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGa                >  1:1109852/1‑150 (MQ=255)
                                                                                                                                       tATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGAt               <  2:1597393/150‑1 (MQ=255)
                                                                                                                                       tATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGAt               <  2:3859942/150‑1 (MQ=255)
                                                                                                                                        aTTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATg              <  1:4212570/150‑1 (MQ=255)
                                                                                                                                        aTTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATg              >  2:3316161/1‑150 (MQ=255)
                                                                                                                                          tCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGcc            <  1:3707085/150‑1 (MQ=255)
                                                                                                                                          tCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGcc            <  2:3520605/150‑1 (MQ=255)
                                                                                                                                           cgccgcCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCg           >  1:897601/1‑150 (MQ=255)
                                                                                                                                               gccTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAATAGCCTCGCGGGCGGGATGCCGCttt       >  1:2220003/1‑150 (MQ=255)
                                                                                                                                               gccTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCttt       <  2:254809/150‑1 (MQ=255)
                                                                                                                                                 cTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCg     <  1:3511634/150‑1 (MQ=255)
                                                                                                                                                 cTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCg     <  2:235487/150‑1 (MQ=255)
                                                                                                                                                    tgCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGc  >  1:2320242/1‑150 (MQ=255)
                                                                                                                                                    tgCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGc  >  1:2320361/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACACAGGACTCCCTCGACGCCATCGAACGCTTGTTTAAATCAGACATCGAAGCGAAGCAGGTGGCGGCGATTATTTTCGAACCGGTGCAGGGCGAGGGCGGTTTCAACGTTGCGCCAAAAGAGCTGGTTGCCGCTATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGC  >  CP009273/1360340‑1360637

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: