Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,851,456 G→T 5.2% L103L (CTC→CTA)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,851,4560GT 5.2% 810.9 / 12.9 211L103L (CTC→CTA) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (105/95);  new base T (4/7);  total (109/102)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 220 total aligned reads displayed.

CTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAG  >  CP009273/3851309‑3851604
                                                                                                                                                   |                                                                                                                                                    
cTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCga                                                                                                                                                     <  2:3969027/149‑1 (MQ=255)
   tAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGAt                                                                                                                                                  >  2:3064363/1‑149 (MQ=255)
   tAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGAt                                                                                                                                                  <  1:2276275/149‑1 (MQ=255)
    aaaTCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATg                                                                                                                                                 >  1:3448823/1‑149 (MQ=255)
      aTCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGcc                                                                                                                                               >  2:612520/1‑149 (MQ=255)
       tCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCa                                                                                                                                              >  1:62182/1‑149 (MQ=255)
       tCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCa                                                                                                                                              >  1:1566786/1‑149 (MQ=255)
               gCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGa                                                                                                                                      >  2:25221/1‑149 (MQ=255)
               gCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGa                                                                                                                                      >  1:1430263/1‑149 (MQ=255)
                cTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGAt                                                                                                                                     <  2:1405793/149‑1 (MQ=255)
                cTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGAt                                                                                                                                     >  2:1771415/1‑149 (MQ=255)
                cTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGAt                                                                                                                                     >  2:507872/1‑149 (MQ=255)
                  gTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGt                                                                                                                                   <  1:1228475/149‑1 (MQ=255)
                   tttAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTa                                                                                                                                  >  1:181335/1‑149 (MQ=255)
                   tttAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTa                                                                                                                                  <  2:2448146/149‑1 (MQ=255)
                      aGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCGGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACa                                                                                                                               >  1:2420898/1‑149 (MQ=255)
                      aGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACa                                                                                                                               <  2:888060/149‑1 (MQ=255)
                      aGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACa                                                                                                                               >  1:349395/1‑149 (MQ=255)
                        gCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAgg                                                                                                                             >  2:1595471/1‑149 (MQ=255)
                          cATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGgtg                                                                                                                           >  1:2123988/1‑149 (MQ=255)
                          cATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGgtg                                                                                                                           <  2:1312590/149‑1 (MQ=255)
                            ttGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGtgtg                                                                                                                         >  2:3205510/1‑149 (MQ=255)
                                 ccACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGc                                                                                                                    >  2:640016/1‑149 (MQ=255)
                                  cacaGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGcc                                                                                                                   <  1:871301/149‑1 (MQ=255)
                                   acaGGCTAATCTCTTTGTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCg                                                                                                                  <  2:3365706/149‑1 (MQ=255)
                                        cTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATc                                                                                                             >  1:3548855/1‑149 (MQ=255)
                                        cTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATc                                                                                                             <  1:3108462/149‑1 (MQ=255)
                                         tAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATcc                                                                                                            >  1:2944300/1‑149 (MQ=255)
                                         tAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATcc                                                                                                            >  1:3042397/1‑149 (MQ=255)
                                         tAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATcc                                                                                                            >  1:3042392/1‑149 (MQ=255)
                                            tctcTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTa                                                                                                         <  2:2951072/149‑1 (MQ=255)
                                             ctctTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTaa                                                                                                        >  2:845994/1‑149 (MQ=255)
                                               ctTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAg                                                                                                      >  2:3965946/1‑149 (MQ=255)
                                               ctTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAg                                                                                                      >  1:1725253/1‑149 (MQ=255)
                                                     ccGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAgcgc                                                                                                >  1:286118/1‑149 (MQ=255)
                                                      cGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGcc                                                                                               <  2:3391099/149‑1 (MQ=255)
                                                      cGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGcc                                                                                               >  2:557950/1‑149 (MQ=255)
                                                       gTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGccc                                                                                              >  1:3132327/1‑149 (MQ=255)
                                                           gTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATa                                                                                          >  1:73828/1‑149 (MQ=255)
                                                              cGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTg                                                                                       >  2:1166567/1‑149 (MQ=255)
                                                               gCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTgg                                                                                      >  2:3164607/1‑149 (MQ=255)
                                                                cTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGga                                                                                     >  1:1627755/1‑149 (MQ=255)
                                                                cTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCAGAGTGga                                                                                     >  2:945918/1‑149 (MQ=255)
                                                                    gATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATg                                                                                 >  1:2328923/1‑149 (MQ=255)
                                                                      ttgTTCGCACCATCTAACCATTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTTCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGtt                                                                               <  1:1893753/147‑1 (MQ=255)
                                                                        gTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGtttt                                                                             <  2:288626/149‑1 (MQ=255)
                                                                        gTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCACACCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGtttt                                                                             >  1:1462537/1‑149 (MQ=255)
                                                                         ttCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGttttt                                                                            <  2:2376404/149‑1 (MQ=255)
                                                                          tCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTg                                                                           >  2:697781/1‑149 (MQ=255)
                                                                           cGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGc                                                                          <  1:2553544/149‑1 (MQ=255)
                                                                            gCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGcc                                                                         >  2:3694909/1‑149 (MQ=255)
                                                                             ccccATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCCTAGTGGAGATGTTTTTGCCt                                                                        >  2:3040612/1‑149 (MQ=255)
                                                                              cccATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTg                                                                       <  1:3770510/149‑1 (MQ=255)
                                                                              cccATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTg                                                                       >  2:3310871/1‑149 (MQ=255)
                                                                               ccATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTgg                                                                      <  2:4075371/149‑1 (MQ=255)
                                                                                cATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGgc                                                                     >  1:537094/1‑149 (MQ=255)
                                                                                cATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGgc                                                                     >  1:3254505/1‑149 (MQ=255)
                                                                                     aaCCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGAc                                                                >  2:2667909/1‑149 (MQ=255)
                                                                                       ccAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACgt                                                              >  2:2262924/1‑149 (MQ=255)
                                                                                           tAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACgttgtt                                                          >  1:3372350/1‑149 (MQ=255)
                                                                                             aTCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTgg                                                        >  1:2246671/1‑149 (MQ=255)
                                                                                               cAGCGTCCCACTCCGGCGGACACTCTCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTc                                                      <  1:1073308/149‑1 (MQ=255)
                                                                                                 gCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCa                                                    >  1:1067019/1‑149 (MQ=255)
                                                                                                 gCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCa                                                    >  2:1420443/1‑149 (MQ=255)
                                                                                                  cGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCAc                                                   <  2:1044155/149‑1 (MQ=255)
                                                                                                      ccACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGc                                               >  1:2993644/1‑149 (MQ=255)
                                                                                                      ccACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGc                                               <  1:3966020/149‑1 (MQ=255)
                                                                                                        aCTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATCACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGACCCATAGTGGAGATGTTTTTGCCTGGCGAGACGTTGTTGGTCCTCGGGCCg                                             >  2:214842/1‑149 (MQ=255)
                                                                                                          tCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGa                                           >  2:1667340/1‑149 (MQ=255)
                                                                                                            cggcggACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACt                                         <  1:1484586/149‑1 (MQ=255)
                                                                                                            cggcggACACTCGCCAGTGCCGAAATAGTCAGGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGTTAGCCGGCATCCTTAAAGGAGCGCCCCATAGGGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACt                                         <  1:3520233/149‑1 (MQ=255)
                                                                                                             ggcggACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTCAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTg                                        <  2:3775670/149‑1 (MQ=255)
                                                                                                                     cTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGt                                <  1:58206/149‑1 (MQ=255)
                                                                                                                      tCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTa                               <  2:3832036/149‑1 (MQ=255)
                                                                                                                         ccAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGAt                            <  2:62182/149‑1 (MQ=255)
                                                                                                                         ccAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGAt                            <  2:1566786/149‑1 (MQ=255)
                                                                                                                          cAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATa                           <  1:639984/149‑1 (MQ=255)
                                                                                                                          cAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATa                           <  1:640370/149‑1 (MQ=255)
                                                                                                                           aGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATAc                          >  2:1007831/1‑149 (MQ=255)
                                                                                                                           aGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATAc                          <  1:845994/149‑1 (MQ=255)
                                                                                                                             tGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCg                        >  1:2143307/1‑149 (MQ=255)
                                                                                                                             tGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCg                        >  2:820391/1‑149 (MQ=255)
                                                                                                                               ccGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGt                      <  2:3547634/149‑1 (MQ=255)
                                                                                                                                  aaaTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTaaa                   >  1:1266238/1‑149 (MQ=255)
                                                                                                                                  aaaTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTaaa                   >  1:709144/1‑149 (MQ=255)
                                                                                                                                    aTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATa                 >  1:1119798/1‑149 (MQ=255)
                                                                                                                                    aTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATa                 >  1:3176435/1‑149 (MQ=255)
                                                                                                                                     tAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAg                <  2:3748902/149‑1 (MQ=255)
                                                                                                                                     tAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAg                <  2:2045029/149‑1 (MQ=255)
                                                                                                                                     gaGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAg                <  2:1703767/148‑1 (MQ=255)
                                                                                                                                       gTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTc              <  1:482058/149‑1 (MQ=255)
                                                                                                                                       gTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTc              >  2:1871863/1‑149 (MQ=255)
                                                                                                                                            gACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGcc         >  2:2072783/1‑149 (MQ=255)
                                                                                                                                             aCCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGccc        <  2:181335/149‑1 (MQ=255)
                                                                                                                                              ccGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCg       <  2:2230242/149‑1 (MQ=255)
                                                                                                                                              ccGTCGAGATGCCATGTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCg       <  2:3658991/149‑1 (MQ=255)
                                                                                                                                               cGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGt      >  2:4073260/1‑149 (MQ=255)
                                                                                                                                                gTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTg     >  2:3071521/1‑149 (MQ=255)
                                                                                                                                                   gagaTGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAg  >  1:411122/1‑149 (MQ=255)
                                                                                                                                                   gagaTGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAg  >  2:1602691/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAG  >  CP009273/3851309‑3851604

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: ‑