Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,882,377 C→A 93.0% T96T (ACC→ACA tnaA → tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,882,3770CA93.0% 601.7 / 35.4 217T96T (ACC→ACAtnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (10/5);  new base A (82/118);  total (94/123)
Fisher's exact test for biased strand distribution p-value = 6.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 217 total aligned reads displayed.

TCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAG  >  CP009273/3882232‑3882521
                                                                                                                                                 |                                                                                                                                                
tCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAAt                                                                                                                                                <  2:1473637/148‑1 (MQ=255)
tCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAAt                                                                                                                                                >  2:1813959/1‑148 (MQ=255)
 cGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAAtt                                                                                                                                               <  2:1160484/148‑1 (MQ=255)
 cGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAAtt                                                                                                                                               <  1:184436/148‑1 (MQ=255)
  gATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTc                                                                                                                                              >  1:2582338/1‑148 (MQ=255)
   aTTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTcc                                                                                                                                             >  1:1503836/1‑148 (MQ=255)
    tttACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCg                                                                                                                                            <  1:3315953/148‑1 (MQ=255)
        cTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTc                                                                                                                                        <  1:2321470/148‑1 (MQ=255)
        cTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTc                                                                                                                                        <  2:190703/148‑1 (MQ=255)
          gACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCAc                                                                                                                                      <  1:2369203/148‑1 (MQ=255)
           aCCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCAcc                                                                                                                                     >  2:3037739/1‑148 (MQ=255)
              gACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAgg                                                                                                                                  >  1:126518/1‑148 (MQ=255)
              gACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAgg                                                                                                                                  <  2:931773/148‑1 (MQ=255)
               aCAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAggg                                                                                                                                 >  2:1656755/1‑148 (MQ=255)
               aCAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAggg                                                                                                                                 >  2:2061722/1‑148 (MQ=255)
                cAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGc                                                                                                                                >  1:2400203/1‑148 (MQ=255)
                cAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGc                                                                                                                                >  1:2397547/1‑148 (MQ=255)
                  gCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCg                                                                                                                              <  2:3241667/148‑1 (MQ=255)
                   cGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGt                                                                                                                             >  1:670309/1‑148 (MQ=255)
                   cGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGt                                                                                                                             >  1:2813902/1‑148 (MQ=255)
                         cGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTATCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACACTTCCGACTCACCAGGGCCGTGGCGCa                                                                                                                       <  1:1903050/148‑1 (MQ=255)
                          ggggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCag                                                                                                                      <  2:1262876/148‑1 (MQ=255)
                            ggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                                    <  2:1899272/148‑1 (MQ=255)
                               ggTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCaa                                                                                                                 <  1:3002927/148‑1 (MQ=255)
                               ggTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCaa                                                                                                                 <  2:2583791/148‑1 (MQ=255)
                                gTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCaaa                                                                                                                <  1:308482/148‑1 (MQ=255)
                                 tGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAAt                                                                                                               <  2:1841095/148‑1 (MQ=255)
                                   aCGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCt                                                                                                             >  1:30943/1‑148 (MQ=255)
                                   aCGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCt                                                                                                             <  2:1376600/148‑1 (MQ=255)
                                        gagCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATAtt                                                                                                        <  1:1050663/148‑1 (MQ=255)
                                         agCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTc                                                                                                       <  1:809696/148‑1 (MQ=255)
                                          gCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTcc                                                                                                      >  2:2391632/1‑148 (MQ=255)
                                           cATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCg                                                                                                     >  1:1044304/1‑148 (MQ=255)
                                           cATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCg                                                                                                     >  2:464379/1‑148 (MQ=255)
                                           cATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCg                                                                                                     >  1:626065/1‑148 (MQ=255)
                                             tGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGt                                                                                                   <  2:438129/148‑1 (MQ=255)
                                             tGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGt                                                                                                   >  2:877880/1‑148 (MQ=255)
                                              gCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTa                                                                                                  <  1:2154429/148‑1 (MQ=255)
                                                 ggCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTg                                                                                               >  2:296193/1‑148 (MQ=255)
                                                   cTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGAt                                                                                             <  2:876443/148‑1 (MQ=255)
                                                      cGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTaa                                                                                          >  1:1924892/1‑148 (MQ=255)
                                                       gatgatGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTaaa                                                                                         <  2:3058333/148‑1 (MQ=255)
                                                         tgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTaaaaa                                                                                       >  1:1960507/1‑148 (MQ=255)
                                                         tgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTaaaaa                                                                                       <  1:71010/148‑1 (MQ=255)
                                                          gatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTaaaaaa                                                                                      <  1:1448159/148‑1 (MQ=255)
                                                               gcggcgACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGa                                                                                 <  2:270014/148‑1 (MQ=255)
                                                                             tACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGcc                                                                   >  2:3055271/1‑148 (MQ=255)
                                                                             tACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGcc                                                                   <  1:2773797/148‑1 (MQ=255)
                                                                             tACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGcc                                                                   <  2:2968924/148‑1 (MQ=255)
                                                                                  cGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGAt                                                              >  2:574612/1‑148 (MQ=255)
                                                                                     cAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCgc                                                           >  2:3468991/1‑148 (MQ=255)
                                                                                             gCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAt                                                   <  2:1112868/148‑1 (MQ=255)
                                                                                             gCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAt                                                   <  2:3117131/148‑1 (MQ=255)
                                                                                             gCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAt                                                   <  2:2865238/148‑1 (MQ=255)
                                                                                             gCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAt                                                   >  2:1123641/1‑148 (MQ=255)
                                                                                                 ctaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggtg                                               <  2:3322607/148‑1 (MQ=255)
                                                                                                  taTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggtgg                                              >  1:1321757/1‑148 (MQ=255)
                                                                                                     gCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGt                                           >  1:1660354/1‑148 (MQ=255)
                                                                                                      cGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGtt                                          >  2:3513720/1‑148 (MQ=255)
                                                                                                      cGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGtt                                          >  2:427253/1‑148 (MQ=255)
                                                                                                          aGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTctct                                      <  2:989836/148‑1 (MQ=255)
                                                                                                           gCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTa                                     >  1:2455931/1‑148 (MQ=255)
                                                                                                           gCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTa                                     <  1:3470796/148‑1 (MQ=255)
                                                                                                            ccGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTaa                                    <  2:2684631/148‑1 (MQ=255)
                                                                                                             cGAGTCAGTGAAAAATATCTTTGTTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATATATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCTCAGCAAAATGGTGGCGTTCTCTAAc                                   <  2:36519/148‑1 (MQ=255)
                                                                                                             cGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc                                   >  1:570375/1‑148 (MQ=255)
                                                                                                              gAGTCAGTGAAAAATATCTTTGGTAATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACt                                  <  2:707233/148‑1 (MQ=255)
                                                                                                                  cAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAttt                              >  1:1438197/1‑148 (MQ=255)
                                                                                                                  cAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAttt                              <  2:2242428/148‑1 (MQ=255)
                                                                                                                  cAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAttt                              <  2:1609825/148‑1 (MQ=255)
                                                                                                                   aGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAACATGGTGGCGTTCTCTAACTATTTc                             <  1:3037739/148‑1 (MQ=255)
                                                                                                                    gTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCt                            <  2:2770191/148‑1 (MQ=255)
                                                                                                                     tGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGTCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCtt                           >  1:115148/1‑148 (MQ=255)
                                                                                                                      gAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCttt                          >  1:1516320/1‑148 (MQ=255)
                                                                                                                       aaaaaTATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTg                         >  1:2481337/1‑148 (MQ=255)
                                                                                                                        aaaaTATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                        >  2:3363317/1‑148 (MQ=255)
                                                                                                                         aaaTATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGAt                       >  1:961582/1‑148 (MQ=255)
                                                                                                                          aaTATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATa                      >  1:2065670/1‑148 (MQ=255)
                                                                                                                           atatCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAc                     <  2:2815230/148‑1 (MQ=255)
                                                                                                                           atatCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAc                     >  1:52687/1‑148 (MQ=255)
                                                                                                                           atatCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAc                     >  1:398342/1‑148 (MQ=255)
                                                                                                                            tatCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAcc                    >  1:174474/1‑148 (MQ=255)
                                                                                                                            tatCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAcc                    <  2:2379651/148‑1 (MQ=255)
                                                                                                                             atCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                   <  1:2838064/148‑1 (MQ=255)
                                                                                                                             atCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                   >  1:389909/1‑148 (MQ=255)
                                                                                                                              tCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCAc                  >  2:1565627/1‑148 (MQ=255)
                                                                                                                              tCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCAc                  >  1:928622/1‑148 (MQ=255)
                                                                                                                                 ttGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCACAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa               >  2:2498905/1‑148 (MQ=255)
                                                                                                                                 ttGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa               >  2:2206742/1‑148 (MQ=255)
                                                                                                                                 ttGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa               >  2:410870/1‑148 (MQ=255)
                                                                                                                                   ggTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAgg             <  1:2941912/148‑1 (MQ=255)
                                                                                                                                   ggTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAgg             <  1:2941720/148‑1 (MQ=255)
                                                                                                                                     ttATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGc           <  2:645044/148‑1 (MQ=255)
                                                                                                                                       aTCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCa         <  1:431086/148‑1 (MQ=255)
                                                                                                                                        tCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCAt        <  2:2790042/148‑1 (MQ=255)
                                                                                                                                         cAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATa       <  1:121591/148‑1 (MQ=255)
                                                                                                                                         cAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATa       <  1:11778/148‑1 (MQ=255)
                                                                                                                                         cAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATa       >  1:1054356/1‑148 (MQ=255)
                                                                                                                                           aTACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGc     >  1:3511019/1‑148 (MQ=255)
                                                                                                                                            tACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCCGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGcc    <  2:2080274/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAG  >  CP009273/3882232‑3882521

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: