Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,913,709 C→A 19.3% E133* (GAA→TAA)  gudX ← glucarate dehydratase‑related protein, substrate unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,913,7090CA19.3% 680.7 / 103.8 247E133* (GAA→TAA) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base C (114/85);  new base A (30/18);  total (144/103)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

GCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCT  >  CP009273/2913564‑2913855
                                                                                                                                                 |                                                                                                                                                  
gCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCa                                                                                                                                                 >  2:337483/1‑149 (MQ=255)
gCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCa                                                                                                                                                 <  1:2891285/149‑1 (MQ=255)
    tCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGAcc                                                                                                                                             >  1:2189459/1‑149 (MQ=255)
    tCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGAcc                                                                                                                                             >  1:1358563/1‑149 (MQ=255)
       gATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGa                                                                                                                                          >  2:384510/1‑149 (MQ=255)
         tGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAAc                                                                                                                                        <  1:3583139/149‑1 (MQ=255)
           ccGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACAt                                                                                                                                      <  2:1664136/149‑1 (MQ=255)
             gCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTc                                                                                                                                    <  2:2762033/149‑1 (MQ=255)
              cAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCa                                                                                                                                   >  1:863607/1‑149 (MQ=255)
              cAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCa                                                                                                                                   <  2:2910214/149‑1 (MQ=255)
               aaCTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAg                                                                                                                                  <  1:1152987/149‑1 (MQ=255)
                aCTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAgc                                                                                                                                 >  1:3702925/1‑149 (MQ=255)
                  tGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAgcgc                                                                                                                               >  1:687321/1‑149 (MQ=255)
                    aTACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCt                                                                                                                             >  1:2592454/1‑149 (MQ=255)
                        cACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTAcc                                                                                                                         <  1:2017895/149‑1 (MQ=255)
                        cACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTAcc                                                                                                                         <  1:1343336/149‑1 (MQ=255)
                        cACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTAcc                                                                                                                         <  2:558281/149‑1 (MQ=255)
                            cATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGc                                                                                                                     >  1:3988080/1‑149 (MQ=255)
                              tGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAg                                                                                                                   >  2:2974773/1‑149 (MQ=255)
                                      cGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCa                                                                                                           >  1:1740728/1‑149 (MQ=255)
                                      cGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCa                                                                                                           >  2:3920373/1‑149 (MQ=255)
                                      cGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCa                                                                                                           >  2:4015214/1‑149 (MQ=255)
                                       ggCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCaa                                                                                                          >  2:329178/1‑149 (MQ=255)
                                        gCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCaaa                                                                                                         >  1:3338723/1‑149 (MQ=255)
                                          gTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGc                                                                                                       <  2:1047639/149‑1 (MQ=255)
                                           tATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCg                                                                                                      >  1:2691769/1‑149 (MQ=255)
                                           tATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCg                                                                                                      >  1:754023/1‑149 (MQ=255)
                                              tttCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCt                                                                                                   <  2:126102/149‑1 (MQ=255)
                                                  cacaTAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCa                                                                                               >  2:3755977/1‑149 (MQ=255)
                                                         ggAAGATCGGTTTTGGTCCGATCAACGATAGAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAATCGGCTTCCAGcgccgc                                                                                        >  2:682867/1‑149 (MQ=255)
                                                            aGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCAc                                                                                     >  1:1118402/1‑149 (MQ=255)
                                                            aGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCAc                                                                                     <  1:265343/149‑1 (MQ=255)
                                                              aTCCGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACgg                                                                                   >  2:3638452/1‑149 (MQ=255)
                                                                  gTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGtt                                                                               >  2:1644380/1‑149 (MQ=255)
                                                                  gTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGtt                                                                               <  2:2490023/149‑1 (MQ=255)
                                                                   ttttGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTa                                                                              >  2:799484/1‑149 (MQ=255)
                                                                     ttGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAc                                                                            >  1:2518177/1‑149 (MQ=255)
                                                                     ttGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAc                                                                            <  1:3063374/149‑1 (MQ=255)
                                                                       ggTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAcgc                                                                          <  1:1201934/149‑1 (MQ=255)
                                                                         nccGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGt                                                                        >  1:3897834/2‑149 (MQ=255)
                                                                          ccGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTa                                                                       <  1:2242248/149‑1 (MQ=255)
                                                                           cGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTaa                                                                      >  2:1149200/1‑149 (MQ=255)
                                                                            gATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAAt                                                                     <  1:270525/149‑1 (MQ=255)
                                                                            gATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAAt                                                                     <  1:310200/149‑1 (MQ=255)
                                                                                aCCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCaa                                                                 >  2:2592882/1‑149 (MQ=255)
                                                                                   gATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAgg                                                              >  2:1244831/1‑149 (MQ=255)
                                                                                   gATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAgg                                                              >  2:2457583/1‑149 (MQ=255)
                                                                                     tataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTc                                                            >  1:1011562/1‑149 (MQ=255)
                                                                                     tataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTc                                                            <  1:1518583/149‑1 (MQ=255)
                                                                                      ataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTcc                                                           >  2:2877061/1‑149 (MQ=255)
                                                                                       taAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCa                                                          >  1:1717791/1‑149 (MQ=255)
                                                                                       taAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCa                                                          <  1:1177769/149‑1 (MQ=255)
                                                                                       taAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCa                                                          <  1:2791606/149‑1 (MQ=255)
                                                                                            cAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCAc                                                     >  1:3263612/1‑149 (MQ=255)
                                                                                              gATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCt                                                   >  1:771176/1‑149 (MQ=255)
                                                                                                tAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCttt                                                 >  1:387159/1‑149 (MQ=255)
                                                                                                tAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCttt                                                 <  1:337483/149‑1 (MQ=255)
                                                                                                  aCCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTAc                                               <  1:2339723/149‑1 (MQ=255)
                                                                                                   ccGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTAcc                                              >  1:458810/1‑149 (MQ=255)
                                                                                                     gAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACAGa                                            <  2:349877/149‑1 (MQ=255)
                                                                                                      aGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGaa                                           >  1:2286588/1‑149 (MQ=255)
                                                                                                      aGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGaa                                           >  2:429829/1‑149 (MQ=255)
                                                                                                      aGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGaa                                           <  2:2623298/149‑1 (MQ=255)
                                                                                                        gACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTCACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAgg                                         >  2:1581955/1‑149 (MQ=255)
                                                                                                          cGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTa                                       >  1:562603/1‑149 (MQ=255)
                                                                                                           ggTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGACGGtct                                      >  2:1737944/1‑147 (MQ=255)
                                                                                                            gTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATc                                     <  1:1333344/149‑1 (MQ=255)
                                                                                                             tAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCa                                    >  2:2584429/1‑149 (MQ=255)
                                                                                                                tAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCaaaa                                 <  2:596374/149‑1 (MQ=255)
                                                                                                                 aGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAAt                                >  1:3552073/1‑149 (MQ=255)
                                                                                                                 aGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAAt                                >  1:3804982/1‑149 (MQ=255)
                                                                                                                 aGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAAt                                <  2:3067917/149‑1 (MQ=255)
                                                                                                                   ccTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcg                              <  2:1812474/149‑1 (MQ=255)
                                                                                                                   ccTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcg                              >  1:1853805/1‑149 (MQ=255)
                                                                                                                    cTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcgg                             >  1:2326802/1‑149 (MQ=255)
                                                                                                                    cTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcgg                             >  1:300715/1‑149 (MQ=255)
                                                                                                                    cTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcgg                             <  2:2012788/149‑1 (MQ=255)
                                                                                                                     tCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTAACGAAGGTATCAACATcggc                            >  1:2123617/1‑149 (MQ=255)
                                                                                                                         gTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGcc                        <  1:1851774/149‑1 (MQ=255)
                                                                                                                          ttgctagcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATAAAAATCGGCGGCCt                       >  2:2050684/1‑149 (MQ=255)
                                                                                                                           tgcttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTg                      >  2:854697/1‑149 (MQ=255)
                                                                                                                            gcttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGa                     <  1:2895375/149‑1 (MQ=255)
                                                                                                                            gcttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGa                     >  2:834533/1‑149 (MQ=255)
                                                                                                                            gcttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGa                     <  2:1467029/149‑1 (MQ=255)
                                                                                                                             cttgcCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGAt                    >  2:1084560/1‑149 (MQ=255)
                                                                                                                              ttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTCTCAAAATCGGCGGCCTGAtt                   >  2:2711796/1‑149 (MQ=255)
                                                                                                                              ttgcCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGAtt                   >  2:111735/1‑149 (MQ=255)
                                                                                                                                gcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTAc                 <  1:865453/149‑1 (MQ=255)
                                                                                                                                  cTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCt               >  1:3144665/1‑149 (MQ=255)
                                                                                                                                  cTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAGGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCt               <  2:1358563/149‑1 (MQ=255)
                                                                                                                                  cTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCt               <  2:2189459/149‑1 (MQ=255)
                                                                                                                                  cTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCt               >  1:344168/1‑149 (MQ=255)
                                                                                                                                   tGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCtt              >  2:3584571/1‑149 (MQ=255)
                                                                                                                                   tGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCtt              >  1:3440918/1‑149 (MQ=255)
                                                                                                                                      cccTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTAt           >  2:4040946/1‑149 (MQ=255)
                                                                                                                                      cccTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTAt           >  2:641745/1‑149 (MQ=255)
                                                                                                                                       ccTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATg          >  2:3223863/1‑149 (MQ=255)
                                                                                                                                          aaCAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCAc       >  2:2286212/1‑149 (MQ=255)
                                                                                                                                            cAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAACGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACct     <  1:89634/149‑1 (MQ=255)
                                                                                                                                             aGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACctg    >  2:386696/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCT  >  CP009273/2913564‑2913855

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: