Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 51,352 A→T 8.4% G85G (GGT→GGA apaG ← protein associated with Co2+ and Mg2+ efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927351,3520AT 8.4% 615.3 / 38.8 203G85G (GGT→GGAapaGprotein associated with Co2+ and Mg2+ efflux
Reads supporting (aligned to +/- strand):  ref base A (116/70);  new base T (9/8);  total (125/78)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 204 total aligned reads displayed.

TAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAGATTG  >  CP009273/51207‑51499
                                                                                                                                                 |                                                                                                                                                   
tAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTcc                                                                                                                                                   >  1:1144916/1‑148 (MQ=255)
tAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTcc                                                                                                                                                   >  2:2384823/1‑148 (MQ=255)
tAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTcc                                                                                                                                                   >  1:1606672/1‑148 (MQ=255)
   ggTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCt                                                                                                                                                >  1:1415077/1‑148 (MQ=255)
     tATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTgg                                                                                                                                              >  1:2466980/1‑148 (MQ=255)
        gTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGgtgt                                                                                                                                           >  1:499064/1‑148 (MQ=255)
        gTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGgtgt                                                                                                                                           <  1:426021/148‑1 (MQ=255)
        gTCGACATATTCTTTTAATGAATGAGTGTGTGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGgtgt                                                                                                                                           >  2:1677690/1‑148 (MQ=255)
             cATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTgg                                                                                                                                      >  2:2446878/1‑148 (MQ=255)
               tatTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTa                                                                                                                                    <  1:1314836/148‑1 (MQ=255)
               tatTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTa                                                                                                                                    <  2:2603704/148‑1 (MQ=255)
               tatTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTa                                                                                                                                    <  1:1081498/148‑1 (MQ=255)
               tatTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTa                                                                                                                                    >  1:893459/1‑148 (MQ=255)
                atTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTAc                                                                                                                                   >  2:2778756/1‑148 (MQ=255)
                atTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTAc                                                                                                                                   >  1:2551579/1‑148 (MQ=255)
                 ttCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTAct                                                                                                                                  <  2:2623651/148‑1 (MQ=255)
                 ttCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTAct                                                                                                                                  >  2:413381/1‑148 (MQ=255)
                 ttCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTAct                                                                                                                                  >  1:2884991/1‑148 (MQ=255)
                   cTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTActct                                                                                                                                >  2:2416989/1‑148 (MQ=255)
                    ttttAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCtt                                                                                                                               >  1:3176130/1‑148 (MQ=255)
                         atgaatgaGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGACATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCa                                                                                                                          >  2:74472/1‑148 (MQ=255)
                          tgaatgaGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAg                                                                                                                         >  1:1357981/1‑148 (MQ=255)
                           gaatgaGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAgg                                                                                                                        >  2:1166092/1‑148 (MQ=255)
                           gaatgaGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAgg                                                                                                                        >  2:1463654/1‑148 (MQ=255)
                              tgaGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGgcgc                                                                                                                     >  2:84889/1‑148 (MQ=255)
                                  tgtgGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATa                                                                                                                 >  2:954070/1‑148 (MQ=255)
                                    tgGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAg                                                                                                               <  2:2009398/148‑1 (MQ=255)
                                       gAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTgg                                                                                                            >  2:1822902/1‑148 (MQ=255)
                                            gCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGa                                                                                                       <  2:307110/148‑1 (MQ=255)
                                            gCGAGTCGGAATACGGGAATGTCGATGATGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGa                                                                                                       <  1:3224484/148‑1 (MQ=255)
                                                 tCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCa                                                                                                  >  1:2965417/1‑148 (MQ=255)
                                                     aaTACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCa                                                                                              >  2:354868/1‑148 (MQ=255)
                                                     aaTACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCa                                                                                              <  1:168536/148‑1 (MQ=255)
                                                     aaTACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCa                                                                                              <  1:1976381/148‑1 (MQ=255)
                                                         cGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTcc                                                                                          <  1:231879/148‑1 (MQ=255)
                                                             aaTGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCg                                                                                      >  1:1919806/1‑148 (MQ=255)
                                                               tGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGcc                                                                                    <  2:368914/148‑1 (MQ=255)
                                                                gTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGccc                                                                                   >  2:2276702/1‑148 (MQ=255)
                                                                gTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGccc                                                                                   <  2:2617782/148‑1 (MQ=255)
                                                                 tCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCt                                                                                  >  2:2351321/1‑148 (MQ=255)
                                                                     tGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGAc                                                                              >  2:3177287/1‑148 (MQ=255)
                                                                      gCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACt                                                                             >  2:2850563/1‑148 (MQ=255)
                                                                      gCTGAAAGGGACGCCATTTTCAACGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGAAGCCAACCACTCCTTCGCCCTGGACt                                                                             >  1:71583/1‑148 (MQ=255)
                                                                       cTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACtt                                                                            <  1:390878/148‑1 (MQ=255)
                                                                         gAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCg                                                                          >  1:3080623/1‑148 (MQ=255)
                                                                          aaaGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCgg                                                                         <  1:946816/148‑1 (MQ=255)
                                                                          aaaGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCgg                                                                         <  1:358386/148‑1 (MQ=255)
                                                                          aaaGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCgg                                                                         >  1:1523573/1‑148 (MQ=255)
                                                                          aaaGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCGGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCgg                                                                         <  1:2778756/148‑1 (MQ=255)
                                                                             gggACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGttt                                                                      >  2:1616270/1‑148 (MQ=255)
                                                                             gggACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGttt                                                                      >  2:3025804/1‑148 (MQ=255)
                                                                              ggACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTc                                                                     >  2:1689673/1‑148 (MQ=255)
                                                                                   ccATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgc                                                                <  1:1674442/148‑1 (MQ=255)
                                                                                   ccATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgc                                                                <  1:84889/148‑1 (MQ=255)
                                                                                   ccATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgc                                                                >  1:55329/1‑148 (MQ=255)
                                                                                        ttCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattg                                                           <  2:1760965/148‑1 (MQ=255)
                                                                                         tCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgc                                                          <  1:477116/148‑1 (MQ=255)
                                                                                            tCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgccat                                                       >  2:800503/1‑148 (MQ=255)
                                                                                            tCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgccat                                                       >  1:2044647/1‑148 (MQ=255)
                                                                                             cGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgccatt                                                      >  2:1882726/1‑148 (MQ=255)
                                                                                             cGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGgccattgccatt                                                      >  1:3046801/1‑148 (MQ=255)
                                                                                              gATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGccattgccattg                                                     <  1:727633/148‑1 (MQ=255)
                                                                                                tCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGt                                                   <  2:2117219/148‑1 (MQ=255)
                                                                                                tCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGt                                                   >  1:2478650/1‑148 (MQ=255)
                                                                                                tCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGt                                                   >  1:2617372/1‑148 (MQ=255)
                                                                                                tCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGt                                                   >  1:2664300/1‑148 (MQ=255)
                                                                                                tCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGt                                                   >  1:1011434/1‑148 (MQ=255)
                                                                                                 cATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTg                                                  >  2:1378112/1‑148 (MQ=255)
                                                                                                 cATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTg                                                  >  1:1864183/1‑148 (MQ=255)
                                                                                                 cATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTg                                                  >  2:1989457/1‑148 (MQ=255)
                                                                                                   tttCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGAt                                                >  1:1027822/1‑148 (MQ=255)
                                                                                                   tttCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGAt                                                <  2:471545/148‑1 (MQ=255)
                                                                                                    ttCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATc                                               <  2:2290164/148‑1 (MQ=255)
                                                                                                      cGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAg                                             <  2:1042533/148‑1 (MQ=255)
                                                                                                            gACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTa                                       <  1:1230247/148‑1 (MQ=255)
                                                                                                            gACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTa                                       <  2:529895/148‑1 (MQ=255)
                                                                                                                cTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGc                                   >  1:604259/1‑148 (MQ=255)
                                                                                                                  gCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGccc                                 >  2:1898808/1‑148 (MQ=255)
                                                                                                                  gCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGccc                                 <  1:2424930/148‑1 (MQ=255)
                                                                                                                  gCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGccc                                 <  2:2485831/148‑1 (MQ=255)
                                                                                                                     tGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaa                              <  1:1886543/148‑1 (MQ=255)
                                                                                                                     tGGTGCCCAGCGGGGTTTCAATGATTGCTCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaa                              <  1:217657/148‑1 (MQ=255)
                                                                                                                     tGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaa                              >  2:2817422/1‑148 (MQ=255)
                                                                                                                       gTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACAACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaaca                            <  1:2282460/148‑1 (MQ=255)
                                                                                                                        tGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaacaa                           >  2:619355/1‑148 (MQ=255)
                                                                                                                        tGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCcaacaa                           >  2:619219/1‑148 (MQ=255)
                                                                                                                           ccAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTg                        >  2:252818/1‑148 (MQ=255)
                                                                                                                            cAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGc                       >  2:1074451/1‑148 (MQ=255)
                                                                                                                              gCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCAc                     >  2:1149509/1‑148 (MQ=255)
                                                                                                                                ggggTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTg                   <  1:2401015/148‑1 (MQ=255)
                                                                                                                                    tttCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGgcgc               >  1:2332860/1‑148 (MQ=255)
                                                                                                                                     ttCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCt              >  1:255505/1‑148 (MQ=255)
                                                                                                                                       cAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCg            >  2:1627442/1‑148 (MQ=255)
                                                                                                                                          tGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGccc         <  2:124893/148‑1 (MQ=255)
                                                                                                                                          tGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGccc         <  2:1235790/148‑1 (MQ=255)
                                                                                                                                          tGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGccc         <  1:1247441/148‑1 (MQ=255)
                                                                                                                                          tGAGTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGccc         <  2:2884981/148‑1 (MQ=255)
                                                                                                                                              tGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAGa     <  2:2904105/148‑1 (MQ=255)
                                                                                                                                               gCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAGAt    <  1:2209379/148‑1 (MQ=255)
                                                                                                                                                 aCCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAGATTg  >  2:2706385/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TAAGGTATGTCGCCATATTCTTTTAATGAATGAGTGTGGGAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGCCCTGGACTTCGGTTTCACGGCCATTGCCATTGGTGATCAGCCAGTAACGCCCCAACAACTGCACTGGCGCTCGCCCCAGATTG  >  CP009273/51207‑51499

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: