Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,882,377 C→A 45.5% 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,3770CA45.5% 180.8 / 271.5 232T96T (ACC→ACAtnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (66/59);  new base A (43/61);  total (109/123)
Fisher's exact test for biased strand distribution p-value = 8.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

TATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATC  >  CP009273/3882230‑3882524
                                                                                                                                                   |                                                                                                                                                   
tATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACAcc                                                                                                                                                     <  2:3451518/148‑1 (MQ=255)
 aTCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCa                                                                                                                                                    <  2:1765425/148‑1 (MQ=255)
 aTCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCa                                                                                                                                                    <  2:1241897/148‑1 (MQ=255)
  tCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAAt                                                                                                                                                   <  2:563541/148‑1 (MQ=255)
   cGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCAtt                                                                                                                                                  >  2:1263535/1‑148 (MQ=255)
   cGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAAtt                                                                                                                                                  >  1:769896/1‑148 (MQ=255)
    gATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTc                                                                                                                                                 <  1:2304681/148‑1 (MQ=255)
     aTTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTcc                                                                                                                                                >  1:790310/1‑148 (MQ=255)
     aTTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTCAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTcc                                                                                                                                                >  2:334001/1‑148 (MQ=255)
            gACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCAc                                                                                                                                         <  1:1766570/148‑1 (MQ=255)
               cGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGACAAATATCTTTGGTTATCAATACACAATTCCGACTCAcccg                                                                                                                                      >  2:2383653/1‑146 (MQ=255)
               cGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGACAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAg                                                                                                                                      >  2:2383879/1‑148 (MQ=255)
                 aCAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAggg                                                                                                                                    <  1:2903201/148‑1 (MQ=255)
                  cAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGc                                                                                                                                   >  2:3521539/1‑148 (MQ=255)
                    gCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCg                                                                                                                                 >  1:3336137/1‑148 (MQ=255)
                       gCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCTGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTgg                                                                                                                              <  2:704231/148‑1 (MQ=255)
                             gggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCaga                                                                                                                        >  1:1664482/1‑148 (MQ=255)
                             gggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCaga                                                                                                                        >  1:3302561/1‑148 (MQ=255)
                                cGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCa                                                                                                                     <  2:133989/148‑1 (MQ=255)
                                cGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCa                                                                                                                     <  2:652597/148‑1 (MQ=255)
                                  gTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCaaa                                                                                                                   <  1:2733990/148‑1 (MQ=255)
                                  gTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCaaa                                                                                                                   <  1:455690/148‑1 (MQ=255)
                                  gTGACGCAGAGCATGCAGGCAGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCaaa                                                                                                                   <  2:957230/148‑1 (MQ=255)
                                   tGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAAt                                                                                                                  >  1:1377081/1‑148 (MQ=255)
                                    gACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATc                                                                                                                 >  1:3068053/1‑148 (MQ=255)
                                      cGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCta                                                                                                               <  1:1382269/148‑1 (MQ=255)
                                         agagCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTatat                                                                                                            <  1:118681/148‑1 (MQ=255)
                                           agCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTc                                                                                                          <  1:3521539/148‑1 (MQ=255)
                                             cATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCg                                                                                                        >  1:371942/1‑148 (MQ=255)
                                              aTGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCgg                                                                                                       >  1:2661526/1‑148 (MQ=255)
                                                gCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTa                                                                                                     <  1:1650395/148‑1 (MQ=255)
                                                   ggCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTg                                                                                                  <  1:3356838/148‑1 (MQ=255)
                                                    gCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGa                                                                                                 <  2:3206248/148‑1 (MQ=255)
                                                    gCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGa                                                                                                 <  2:1719837/148‑1 (MQ=255)
                                                     cTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGAt                                                                                                <  1:1086056/148‑1 (MQ=255)
                                                          atgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTaaaa                                                                                           >  1:2961673/1‑148 (MQ=255)
                                                            gatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTaaaaaa                                                                                         <  1:2619427/148‑1 (MQ=255)
                                                                cgcgGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCAGAGTCAGTGAAAAATATCTTTGTTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAAcgcg                                                                                     <  1:1929986/148‑1 (MQ=255)
                                                                 tcggcgACGATTCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAATAATATCTGTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTTCTGATTAAAAAACGCGa                                                                                    <  2:2010393/147‑1 (MQ=255)
                                                                     cgacgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAg                                                                                >  2:2940843/1‑148 (MQ=255)
                                                                        cgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGaa                                                                             <  2:530311/148‑1 (MQ=255)
                                                                          aaGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGaaaa                                                                           >  1:3429195/1‑148 (MQ=255)
                                                                            gCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAg                                                                         >  2:1655863/1‑148 (MQ=255)
                                                                             ccTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAgg                                                                        >  1:1660519/1‑148 (MQ=255)
                                                                              cTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGc                                                                       <  1:2871344/148‑1 (MQ=255)
                                                                               taaaGCTGCAGTCGTAGCTACTAAGCGATAGCCGAGTAAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGAATCACAAGGGCAGTCTCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGAGAGAATAAAAAAAGca                                                                      >  2:2484808/4‑147 (MQ=255)
                                                                               tACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGcc                                                                      <  1:2472671/148‑1 (MQ=255)
                                                                                   gCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGa                                                                  >  1:886893/1‑148 (MQ=255)
                                                                                             tAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCaaa                                                        <  2:1039611/148‑1 (MQ=255)
                                                                                             tAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCaaa                                                        <  1:468125/148‑1 (MQ=255)
                                                                                              aGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCaaaa                                                       <  1:198274/148‑1 (MQ=255)
                                                                                              aGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCaaaa                                                       <  2:929556/148‑1 (MQ=255)
                                                                                               gCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAt                                                      >  1:328108/1‑148 (MQ=255)
                                                                                                ctactaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtg                                                     >  1:3171438/1‑148 (MQ=255)
                                                                                                ctactaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtg                                                     >  1:3478773/1‑148 (MQ=255)
                                                                                                  actaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggt                                                   >  1:2467852/1‑148 (MQ=255)
                                                                                                  actaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggt                                                   <  1:3452821/148‑1 (MQ=255)
                                                                                                   ctaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggtg                                                  <  2:74884/148‑1 (MQ=255)
                                                                                                   ctaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggtg                                                  <  2:77447/148‑1 (MQ=255)
                                                                                                    taTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggtgg                                                 >  1:532323/1‑148 (MQ=255)
                                                                                                       gCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGt                                              >  2:1515895/1‑148 (MQ=255)
                                                                                                        cGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGtt                                             <  1:1699911/148‑1 (MQ=255)
                                                                                                          ttAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtct                                           <  2:3302561/148‑1 (MQ=255)
                                                                                                             gCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTa                                        >  2:3133072/1‑148 (MQ=255)
                                                                                                             gCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTa                                        <  1:2843395/148‑1 (MQ=255)
                                                                                                                   tCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAtt                                  >  2:1761760/1‑148 (MQ=255)
                                                                                                                   tCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAtt                                  >  1:647574/1‑148 (MQ=255)
                                                                                                                    cAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAttt                                 >  2:2210648/1‑148 (MQ=255)
                                                                                                                     aGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTc                                <  2:1952699/148‑1 (MQ=255)
                                                                                                                     aGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTc                                <  2:505960/148‑1 (MQ=255)
                                                                                                                      gTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCt                               <  2:368315/148‑1 (MQ=255)
                                                                                                                       tGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCtt                              >  1:3057819/1‑148 (MQ=255)
                                                                                                                       tGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCtt                              >  2:1182542/1‑148 (MQ=255)
                                                                                                                         aaaaaTATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTg                            <  1:3363776/148‑1 (MQ=255)
                                                                                                                          aaaaTATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                           <  1:1379981/148‑1 (MQ=255)
                                                                                                                          aaaaTATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                           >  1:865705/1‑148 (MQ=255)
                                                                                                                           aaaTATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGAt                          <  1:1951084/148‑1 (MQ=255)
                                                                                                                           aaaTATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGAt                          <  1:3112616/148‑1 (MQ=255)
                                                                                                                           aaaTATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGAt                          >  2:2587509/1‑148 (MQ=255)
                                                                                                                           aaaTATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGAt                          <  1:549899/148‑1 (MQ=255)
                                                                                                                            aaTATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATa                         <  2:1060861/148‑1 (MQ=255)
                                                                                                                              tatCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAcc                       >  1:1020857/1‑148 (MQ=255)
                                                                                                                               atCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                      >  2:360705/1‑148 (MQ=255)
                                                                                                                               atCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                      <  1:905954/148‑1 (MQ=255)
                                                                                                                                tCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCAc                     >  1:2174675/1‑148 (MQ=255)
                                                                                                                                 cTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACg                    >  2:2335760/1‑148 (MQ=255)
                                                                                                                                 cTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACg                    >  2:2536377/1‑148 (MQ=255)
                                                                                                                                   ttGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa                  >  1:2700018/1‑148 (MQ=255)
                                                                                                                                      gTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAggg               <  1:2553762/148‑1 (MQ=255)
                                                                                                                                        tATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGcc             >  2:3512590/1‑148 (MQ=255)
                                                                                                                                          tCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCAt           <  2:872959/148‑1 (MQ=255)
                                                                                                                                          tCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCAt           <  2:3478773/148‑1 (MQ=255)
                                                                                                                                           cAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATa          >  2:2133729/1‑148 (MQ=255)
                                                                                                                                             aTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGc        >  1:666748/1‑148 (MQ=255)
                                                                                                                                              tACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGcc       >  2:189278/1‑148 (MQ=255)
                                                                                                                                               acacAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCa      >  1:3202773/1‑148 (MQ=255)
                                                                                                                                                cacaATTCCGACTCACCCGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAg     >  1:726015/1‑148 (MQ=255)
                                                                                                                                                  ccATTACGACTCACCAGGGCCGTGGCTGAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAACAGGCCTTGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGAt   <  1:154804/148‑1 (MQ=255)
                                                                                                                                                   cATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATc  >  2:2038870/1‑148 (MQ=255)
                                                                                                                                                   cATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATc  >  2:1531137/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATC  >  CP009273/3882230‑3882524

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: