Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,882,020 G→T 16.4% intergenic (+151/‑70) tnaC → / → tnaA tryptophanase leader peptide/tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,882,0200GT16.4% 702.5 / 84.9 225intergenic (+151/‑70)tnaC/tnaAtryptophanase leader peptide/tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (106/82);  new base T (20/17);  total (126/99)
Fisher's exact test for biased strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 231 total aligned reads displayed.

TTGCCCTTCTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCG  >  CP009273/3881871‑3882169
                                                                                                                                                     |                                                                                                                                                     
ttGCCCTTCTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCg                                                                                                                                                       <  1:342742/150‑1 (MQ=255)
ttGCCCTTCTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCg                                                                                                                                                       <  1:2363657/150‑1 (MQ=255)
  gCCCTTCTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGaa                                                                                                                                                     >  1:1752321/1‑150 (MQ=255)
   cccTTCTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATATCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAAt                                                                                                                                                    <  1:1023139/150‑1 (MQ=255)
     cTTCTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTa                                                                                                                                                  >  2:2819639/1‑150 (MQ=255)
     cTTCTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTa                                                                                                                                                  >  1:3689973/1‑150 (MQ=255)
       tCTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAAt                                                                                                                                                <  1:391832/150‑1 (MQ=255)
        cTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATc                                                                                                                                               <  2:454706/150‑1 (MQ=255)
         tGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCg                                                                                                                                              <  1:1639375/150‑1 (MQ=255)
           tAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGt                                                                                                                                            <  2:1118589/150‑1 (MQ=255)
            aGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGta                                                                                                                                           <  2:682761/150‑1 (MQ=255)
                aTCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGc                                                                                                                                       >  2:2356243/1‑150 (MQ=255)
                aTCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGc                                                                                                                                       >  2:1278099/1‑150 (MQ=255)
                aTCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGc                                                                                                                                       >  1:4413018/1‑150 (MQ=255)
                       gagCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTa                                                                                                                                <  2:3009524/150‑1 (MQ=255)
                         gCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAat                                                                                                                              <  2:3890039/150‑1 (MQ=255)
                         gCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAat                                                                                                                              <  1:750608/150‑1 (MQ=255)
                            aaaCCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATAtt                                                                                                                           <  2:3719294/150‑1 (MQ=255)
                                cGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTcaca                                                                                                                       >  2:1677366/1‑150 (MQ=255)
                                 gATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAg                                                                                                                      >  1:282711/1‑150 (MQ=255)
                                  aTTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGTTGTAATATTCTCAgg                                                                                                                     >  2:3954018/1‑150 (MQ=255)
                                     aGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGAt                                                                                                                  <  2:1068790/150‑1 (MQ=255)
                                       aTTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCa                                                                                                                >  2:1168722/1‑150 (MQ=255)
                                        ttCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCAc                                                                                                               >  1:2672684/1‑150 (MQ=255)
                                         tCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACt                                                                                                              <  2:4619622/150‑1 (MQ=255)
                                           aaTGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGt                                                                                                            <  1:1317249/150‑1 (MQ=255)
                                              gtgATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAAt                                                                                                         >  2:826086/1‑150 (MQ=255)
                                                 aTCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTaa                                                                                                      <  1:3615786/150‑1 (MQ=255)
                                                   cTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTaaaa                                                                                                    >  2:3883190/1‑150 (MQ=255)
                                                         gtttgCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATg                                                                                              >  1:2491534/1‑150 (MQ=255)
                                                         gtttgCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATg                                                                                              >  1:3498228/1‑150 (MQ=255)
                                                           ttgCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGaa                                                                                            >  1:3305183/1‑150 (MQ=255)
                                                           ttgCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGaa                                                                                            >  1:4180217/1‑150 (MQ=255)
                                                              cTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGa                                                                                         >  1:3738119/1‑150 (MQ=255)
                                                              cTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGa                                                                                         <  1:484057/150‑1 (MQ=255)
                                                               tataTCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGAt                                                                                        >  2:3401319/1‑150 (MQ=255)
                                                               tataTCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGAt                                                                                        >  2:1082396/1‑150 (MQ=255)
                                                                      tAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTaa                                                                                 >  2:3841261/1‑150 (MQ=255)
                                                                       aaTTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAAt                                                                                <  1:1502515/150‑1 (MQ=255)
                                                                        aTTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATg                                                                               >  2:3357350/1‑150 (MQ=255)
                                                                         ttttGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATgg                                                                              >  1:50221/1‑150 (MQ=255)
                                                                          tttGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGa                                                                             <  2:1497394/150‑1 (MQ=255)
                                                                           ttGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGaa                                                                            <  1:708213/150‑1 (MQ=255)
                                                                            tGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGaaa                                                                           >  1:4600166/1‑150 (MQ=255)
                                                                             gCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATAACGTTAATTAAAATAAATGAAGGATTATGTAATGGaaaa                                                                          >  2:2565837/1‑150 (MQ=255)
                                                                              ccTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAAc                                                                         <  2:144479/150‑1 (MQ=255)
                                                                              ccTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAAc                                                                         >  1:1739391/1‑150 (MQ=255)
                                                                              ccTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAAc                                                                         >  1:1583342/1‑150 (MQ=255)
                                                                                  ttGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTa                                                                     <  1:3262718/150‑1 (MQ=255)
                                                                                  ttGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGCAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTa                                                                     <  2:1093395/150‑1 (MQ=255)
                                                                                    gCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTaaa                                                                   <  1:3516455/150‑1 (MQ=255)
                                                                                     cAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAAc                                                                  <  2:4026560/150‑1 (MQ=255)
                                                                                     cAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAAc                                                                  >  1:4429098/1‑150 (MQ=255)
                                                                                      aaaGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACa                                                                 >  1:278377/1‑150 (MQ=255)
                                                                                       aaGGTCATATCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACAt                                                                >  2:2642083/1‑150 (MQ=255)
                                                                                         ggTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACAtct                                                              <  1:2060551/150‑1 (MQ=255)
                                                                                         ggTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACAtct                                                              <  1:4085924/150‑1 (MQ=255)
                                                                                          gTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACAtctc                                                             >  1:4138663/1‑150 (MQ=255)
                                                                                          gTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACAtctc                                                             >  1:2427001/1‑150 (MQ=255)
                                                                                          gTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACAtctc                                                             <  2:4206841/150‑1 (MQ=255)
                                                                                           tCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTcc                                                            <  1:3209539/150‑1 (MQ=255)
                                                                                           tCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTcc                                                            <  1:3564842/150‑1 (MQ=255)
                                                                                            cATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTccc                                                           <  1:957697/150‑1 (MQ=255)
                                                                                            cATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTccc                                                           >  2:2888539/1‑150 (MQ=255)
                                                                                              tctctcGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTg                                                         <  2:1807933/150‑1 (MQ=255)
                                                                                                 ctcGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAAc                                                      >  1:295909/1‑150 (MQ=255)
                                                                                                 ctcGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAAc                                                      >  2:1940605/1‑150 (MQ=255)
                                                                                                 ctcGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAAc                                                      >  2:1940825/1‑150 (MQ=255)
                                                                                                 ctcGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAAc                                                      >  2:4183337/1‑150 (MQ=255)
                                                                                                  tcGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAAcc                                                     <  1:4187778/150‑1 (MQ=255)
                                                                                                   cgattatttaCAAGATTTAGAAAAAGATGGTCCTTGCATATATATCAAACGAATTAAAAGAAATAGCAGATGTAATAAACACAGGGATCACAGAAATTAAAATAAATGAAGGATTATGTAATGGAAAACATAAAACATCTCCCTGAAcct                                                    >  2:4370539/4‑149 (MQ=255)
                                                                                                   cGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCg                                                    <  1:2372080/150‑1 (MQ=255)
                                                                                                     tttatttaCTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGtt                                                  <  2:3440296/150‑1 (MQ=255)
                                                                                                     tttatttaCTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGtt                                                  <  2:3992289/150‑1 (MQ=255)
                                                                                                       tatttaCTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTcc                                                >  1:3695574/1‑150 (MQ=255)
                                                                                                       tatttaCTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTcc                                                >  1:3817920/1‑150 (MQ=255)
                                                                                                          ttaCTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCa                                             >  1:2274871/1‑150 (MQ=255)
                                                                                                           taCTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCAt                                            >  2:1525708/1‑150 (MQ=255)
                                                                                                                gTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCgtg                                       >  2:1372491/1‑150 (MQ=255)
                                                                                                                gTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCgtg                                       >  2:1372548/1‑150 (MQ=255)
                                                                                                                gTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCgtg                                       >  1:2446895/1‑150 (MQ=255)
                                                                                                                gTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCgtg                                       >  2:1026521/1‑150 (MQ=255)
                                                                                                                  tttAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGtt                                     <  1:2413614/150‑1 (MQ=255)
                                                                                                                     aGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTAtt                                  >  2:3647127/1‑150 (MQ=255)
                                                                                                                       tAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGa                                >  2:2653184/1‑150 (MQ=255)
                                                                                                                           tgatgGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCa                            <  1:977821/150‑1 (MQ=255)
                                                                                                                               ggTGCTTGCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAACATAAATGCAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTaa                        <  1:2122719/150‑1 (MQ=255)
                                                                                                                                      gCATATATATCTGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTAc                 >  1:4520064/1‑150 (MQ=255)
                                                                                                                                       cATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTAcc                >  1:4014748/1‑150 (MQ=255)
                                                                                                                                         tatatatCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCAc              <  2:871679/150‑1 (MQ=255)
                                                                                                                                          atatatCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACt             >  2:217007/1‑150 (MQ=255)
                                                                                                                                           tatatCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTc            >  2:1112027/1‑150 (MQ=255)
                                                                                                                                                cTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCt       <  1:3401319/150‑1 (MQ=255)
                                                                                                                                                cTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCt       >  1:1904272/1‑150 (MQ=255)
                                                                                                                                                 tGGCTAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCtt      >  1:1333182/1‑150 (MQ=255)
                                                                                                                                                 tGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCtt      >  2:1998740/1‑150 (MQ=255)
                                                                                                                                                   gCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTAt    <  1:4292356/150‑1 (MQ=255)
                                                                                                                                                   gCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTAt    <  1:1168722/150‑1 (MQ=255)
                                                                                                                                                     gAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCg  >  2:230185/1‑150 (MQ=255)
                                                                                                                                                     gAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCg  >  1:2746966/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGCCCTTCTGTAGCCATCACCAGAGCCAAACCGATTAGATTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCG  >  CP009273/3881871‑3882169

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: