Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,851,456 G→T 9.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 9.2% 565.4 / 24.1 162L103L (CTC→CTAyidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (88/59);  new base T (7/8);  total (95/67)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 166 total aligned reads displayed.

CTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGC  >  CP009273/3851309‑3851602
                                                                                                                                                   |                                                                                                                                                  
cTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCg                                                                                                                                                    >  1:2135434/1‑148 (MQ=255)
  gTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCgag                                                                                                                                                  >  1:391345/1‑148 (MQ=255)
   tAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCtaga                                                                                                                                                 <  2:1335647/148‑1 (MQ=255)
       tCTTCGACGGTGTTTTGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCTGTTCTCTAAGATAGTTCGCCCCATCTAAGCAGTTATCAGCGTCCCACTCCGGCGGACTCTCGCCAGTGCCGATATAGTCATGACCGTGGAGATGcc                                                                                                                                             <  2:53557/148‑1 (MQ=255)
         ttCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCAt                                                                                                                                           <  1:3141403/148‑1 (MQ=255)
         ttCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCAt                                                                                                                                           <  2:722762/148‑1 (MQ=255)
          tCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCAtt                                                                                                                                          <  2:195129/148‑1 (MQ=255)
          tCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCAtt                                                                                                                                          <  2:195053/148‑1 (MQ=255)
           cGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCAttt                                                                                                                                         <  1:1191668/148‑1 (MQ=255)
           cGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCAttt                                                                                                                                         >  2:1707661/1‑148 (MQ=255)
              cGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGcc                                                                                                                                      >  1:3306177/1‑148 (MQ=255)
                     tAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTaa                                                                                                                               <  1:1784755/148‑1 (MQ=255)
                       ggCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACa                                                                                                                             <  2:2141807/148‑1 (MQ=255)
                             tGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGtgtg                                                                                                                       >  2:2136343/1‑148 (MQ=255)
                                gcCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATa                                                                                                                    <  1:708650/148‑1 (MQ=255)
                                  cacaGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGc                                                                                                                  >  1:3050258/1‑148 (MQ=255)
                                     aGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCgg                                                                                                               >  2:2639626/1‑148 (MQ=255)
                                       gCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCa                                                                                                             >  1:538340/1‑148 (MQ=255)
                                       gCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCa                                                                                                             >  1:150334/1‑148 (MQ=255)
                                        cTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCAt                                                                                                            >  2:3222998/1‑148 (MQ=255)
                                         tAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATCGCCGGCATc                                                                                                           >  2:3427288/1‑148 (MQ=255)
                                          aaTCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATcc                                                                                                          >  1:1559287/1‑148 (MQ=255)
                                          aaTCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATcc                                                                                                          >  2:313206/1‑148 (MQ=255)
                                          aaTCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATcc                                                                                                          >  2:313269/1‑148 (MQ=255)
                                          aaTCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATcc                                                                                                          >  2:3294402/1‑148 (MQ=255)
                                              tctTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTaa                                                                                                      >  1:1538687/1‑148 (MQ=255)
                                              tctTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGGAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTaa                                                                                                      <  2:2019380/148‑1 (MQ=255)
                                               ctTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTaaa                                                                                                     <  2:3380450/148‑1 (MQ=255)
                                                   ttCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAg                                                                                                 >  1:2660586/1‑148 (MQ=255)
                                                   ttCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAg                                                                                                 >  1:1120535/1‑148 (MQ=255)
                                                    tCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAgc                                                                                                <  1:1091912/148‑1 (MQ=255)
                                                     ccGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAgcg                                                                                               >  2:849843/1‑148 (MQ=255)
                                                      cGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAgcgc                                                                                              <  1:580137/148‑1 (MQ=255)
                                                       gTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGcc                                                                                             >  2:2899319/1‑148 (MQ=255)
                                                       gTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGcc                                                                                             <  1:797591/148‑1 (MQ=255)
                                                       gTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGcc                                                                                             >  2:840402/1‑148 (MQ=255)
                                                        tCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGccc                                                                                            >  1:2978951/1‑148 (MQ=255)
                                                         cAGTTCGCTAAGATAGTTAGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGcccc                                                                                           >  2:2300603/1‑148 (MQ=255)
                                                          aGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCa                                                                                          <  1:2807037/148‑1 (MQ=255)
                                                          aGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCa                                                                                          <  1:2868438/148‑1 (MQ=255)
                                                          aGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCa                                                                                          >  1:2249738/1‑148 (MQ=255)
                                                          aGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCa                                                                                          <  1:1414846/148‑1 (MQ=255)
                                                           gTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCAt                                                                                         <  2:2308913/148‑1 (MQ=255)
                                                           gTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCAt                                                                                         >  2:1859118/1‑148 (MQ=255)
                                                             tCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAg                                                                                       <  1:3486875/148‑1 (MQ=255)
                                                             tCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAg                                                                                       >  2:97748/1‑148 (MQ=255)
                                                                cTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTgg                                                                                    >  2:1412462/1‑148 (MQ=255)
                                                                cTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTgg                                                                                    >  2:1413329/1‑148 (MQ=255)
                                                                 tAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGga                                                                                   >  1:717705/1‑148 (MQ=255)
                                                                   aGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGgaga                                                                                 >  1:405982/1‑148 (MQ=255)
                                                                     aTAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATg                                                                               >  1:2747265/1‑148 (MQ=255)
                                                                     aTAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCAGAGTGGAGATg                                                                               >  1:2747031/1‑148 (MQ=255)
                                                                        gTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGttt                                                                            >  1:2075065/1‑148 (MQ=255)
                                                                          tCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGttttt                                                                          >  2:2606943/1‑148 (MQ=255)
                                                                          tCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGttttt                                                                          >  1:837680/1‑148 (MQ=255)
                                                                          tCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGttttt                                                                          >  2:1827902/1‑148 (MQ=255)
                                                                            gCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGc                                                                        <  1:19114/148‑1 (MQ=255)
                                                                             ccccATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGcc                                                                       >  1:944453/1‑148 (MQ=255)
                                                                              cccATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCt                                                                      >  2:1510003/1‑148 (MQ=255)
                                                                               ccATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTg                                                                     >  1:83973/1‑148 (MQ=255)
                                                                                  tCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGgcg                                                                  >  1:1013748/1‑148 (MQ=255)
                                                                                  tCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGgcg                                                                  >  2:3004562/1‑148 (MQ=255)
                                                                                   cGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGgcgc                                                                 >  2:449/1‑148 (MQ=255)
                                                                                   cGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGgcgc                                                                 >  2:432/1‑148 (MQ=255)
                                                                                   cGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGgcgc                                                                 >  2:1193886/1‑148 (MQ=255)
                                                                                       ccAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACg                                                             >  1:1381409/1‑148 (MQ=255)
                                                                                       ccAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACg                                                             <  2:3267492/148‑1 (MQ=255)
                                                                                       ccAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACg                                                             <  1:1022875/148‑1 (MQ=255)
                                                                                          gTAATCAGCGTCCCACTCCGGCGGACACTCGCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACgttgt                                                         <  1:3468284/148‑1 (MQ=255)
                                                                                           tAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACgttgt                                                         >  2:530633/1‑148 (MQ=255)
                                                                                              tCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTgg                                                      >  2:1099093/1‑148 (MQ=255)
                                                                                               cAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGt                                                     <  1:3224590/148‑1 (MQ=255)
                                                                                                     cccACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATCACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCAGAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACgg                                               <  2:146080/148‑1 (MQ=255)
                                                                                                      ccACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACCGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGACCACggg                                              >  2:3543753/1‑148 (MQ=255)
                                                                                                      ccACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACggg                                              >  2:287568/1‑148 (MQ=255)
                                                                                                       cACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGc                                             >  2:2474109/1‑148 (MQ=255)
                                                                                                        aCTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGcc                                            <  1:495359/148‑1 (MQ=255)
                                                                                                        aCTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGcc                                            >  1:2051196/1‑148 (MQ=255)
                                                                                                         cTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCg                                           >  2:2662737/1‑148 (MQ=255)
                                                                                                          tCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCgg                                          >  1:2474261/1‑148 (MQ=255)
                                                                                                              gcggACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTg                                      <  1:2973223/148‑1 (MQ=255)
                                                                                                               cggACACTCGCCAGTGCCGAAATAGTCATGACCGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTgg                                     <  2:2978951/148‑1 (MQ=255)
                                                                                                                ggACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGt                                    <  1:494703/148‑1 (MQ=255)
                                                                                                                 gACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGtt                                   >  1:1933090/1‑148 (MQ=255)
                                                                                                                  acacTCGCCAGTGCCGAAATCGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTa                                  <  2:546970/148‑1 (MQ=255)
                                                                                                                     cTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCg                               >  1:3418915/1‑148 (MQ=255)
                                                                                                                      tCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTGGGGGGGGGGGGCTGGTTAGCGt                              >  2:1625450/1‑148 (MQ=255)
                                                                                                                          cAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGAt                          >  2:3369110/1‑148 (MQ=255)
                                                                                                                             tGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATAcc                       >  1:2595918/1‑148 (MQ=255)
                                                                                                                             tGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATAcc                       >  1:2596440/1‑148 (MQ=255)
                                                                                                                             tGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATAcc                       >  2:496648/1‑148 (MQ=255)
                                                                                                                                 gAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTa                   >  2:2603416/1‑148 (MQ=255)
                                                                                                                                  aaaTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTaa                  >  2:85486/1‑148 (MQ=255)
                                                                                                                                    aTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAAt                <  2:3450455/148‑1 (MQ=255)
                                                                                                                                    aTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAAt                <  2:3450029/148‑1 (MQ=255)
                                                                                                                                        tCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTc            >  2:3319551/1‑148 (MQ=255)
                                                                                                                                               cGTCTAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCg     <  1:2658159/148‑1 (MQ=255)
                                                                                                                                               cGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCg     <  2:717705/148‑1 (MQ=255)
                                                                                                                                               cGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCg     >  2:2240468/1‑148 (MQ=255)
                                                                                                                                               cGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGGAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCg     <  2:305018/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGC  >  CP009273/3851309‑3851602

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: