Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,825,702 C→A 91.2% K732N (AAG→AAT yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,825,7020CA91.2% 960.0 / 76.0 324K732N (AAG→AATyicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (11/18);  new base A (151/144);  total (162/162)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

CGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTC  >  CP009273/3825555‑3825848
                                                                                                                                                   |                                                                                                                                                  
  tACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcg                                                                                                                                                >  1:617973/1‑150 (MQ=255)
     ccGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcct                                                                                                                                             >  1:4918476/1‑150 (MQ=255)
      cGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTtcgcctc                                                                                                                                            <  1:2432274/150‑1 (MQ=255)
      cGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTtcgcctc                                                                                                                                            <  1:1857761/150‑1 (MQ=255)
      cGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctc                                                                                                                                            <  1:1417421/150‑1 (MQ=255)
      cGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctc                                                                                                                                            >  1:4899426/1‑150 (MQ=255)
      cGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctc                                                                                                                                            >  1:333536/1‑150 (MQ=255)
      cGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctc                                                                                                                                            <  2:2942475/150‑1 (MQ=255)
         cATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctcgcc                                                                                                                                         <  2:3614329/150‑1 (MQ=255)
          aTTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGccc                                                                                                                                        <  2:4169587/150‑1 (MQ=255)
              cATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCAc                                                                                                                                    >  1:2097721/1‑150 (MQ=255)
                   aTAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTc                                                                                                                               >  2:3841459/1‑150 (MQ=255)
                    tAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTca                                                                                                                              <  1:2817798/150‑1 (MQ=255)
                     aGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTcac                                                                                                                             <  2:2612251/150‑1 (MQ=255)
                      gATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTcaca                                                                                                                            <  1:2707264/150‑1 (MQ=255)
                      gATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTcaca                                                                                                                            >  1:199083/1‑150 (MQ=255)
                       aTTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAg                                                                                                                           >  1:4240731/1‑150 (MQ=255)
                       aTTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAg                                                                                                                           >  2:334145/1‑150 (MQ=255)
                       aTTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAg                                                                                                                           >  2:3506078/1‑150 (MQ=255)
                        ttACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGt                                                                                                                          <  1:2474472/150‑1 (MQ=255)
                           cAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAAt                                                                                                                       >  2:180295/1‑150 (MQ=255)
                            aaCGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATc                                                                                                                      <  2:3118723/150‑1 (MQ=255)
                               gTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTCATCgtg                                                                                                                   >  2:3411651/1‑150 (MQ=255)
                               gTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCgtg                                                                                                                   >  1:1737844/1‑150 (MQ=255)
                                tAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCgtgt                                                                                                                  <  1:1815220/150‑1 (MQ=255)
                                 aaTTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGtt                                                                                                                 <  1:3165897/150‑1 (MQ=255)
                                    tGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGcc                                                                                                              <  1:1170515/150‑1 (MQ=255)
                                     gTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCa                                                                                                             <  1:4743695/150‑1 (MQ=255)
                                        aGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTa                                                                                                          <  2:187196/150‑1 (MQ=255)
                                         gcgcATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTAc                                                                                                         <  2:953392/150‑1 (MQ=255)
                                         gcgcATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTAc                                                                                                         <  1:4024420/150‑1 (MQ=255)
                                            cATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtg                                                                                                      <  2:2244406/150‑1 (MQ=255)
                                              ttCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtgct                                                                                                    <  2:3718963/150‑1 (MQ=255)
                                               tCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGAGCAGGCACAGTGTCCAGTTCTTCTCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtgctg                                                                                                   >  2:5053266/1‑150 (MQ=255)
                                                  cTTGAGTCTTCACCAAAAGCCGCTGCTCAATTTCAGACTGCGAACCGTATTGCAGACCATTTACTTTCACAACATTGCGCAGGCACTGTGTCCAGTTATTCGCCTCGCCCTCACCAGTCACATTAATAGTGTTGCCTGTACGTGcttctt                                                                                                <  2:2419661/150‑1 (MQ=255)
                                                   ttGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCttt                                                                                               <  1:153276/150‑1 (MQ=255)
                                                     gAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTta                                                                                             >  2:879952/1‑150 (MQ=255)
                                                      aGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaa                                                                                            <  2:579093/150‑1 (MQ=255)
                                                      aGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaa                                                                                            >  1:26098/1‑150 (MQ=255)
                                                        gtttCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaag                                                                                          <  2:1646025/148‑1 (MQ=255)
                                                        gCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaag                                                                                          <  1:3987572/150‑1 (MQ=255)
                                                         cTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaagt                                                                                         >  2:1469853/1‑150 (MQ=255)
                                                         cTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaagt                                                                                         >  2:4211283/1‑150 (MQ=255)
                                                         cTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaagt                                                                                         >  2:4141528/1‑150 (MQ=255)
                                                          ttCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaagta                                                                                        >  2:2880170/1‑150 (MQ=255)
                                                           tCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCAAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaagtaa                                                                                       <  1:883775/150‑1 (MQ=255)
                                                                 ccaGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCa                                                                                 <  2:1277897/150‑1 (MQ=255)
                                                                   aGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCAcc                                                                               >  1:4607402/1‑150 (MQ=255)
                                                                   aGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCAcc                                                                               >  1:2194596/1‑150 (MQ=255)
                                                                   aGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCAcc                                                                               <  1:3751711/150‑1 (MQ=255)
                                                                      cccTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGAt                                                                            <  1:4141484/150‑1 (MQ=255)
                                                                        cTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATcc                                                                          <  2:2045172/150‑1 (MQ=255)
                                                                        cTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATcc                                                                          >  2:3197828/1‑150 (MQ=255)
                                                                        cTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATcc                                                                          >  2:2138610/1‑150 (MQ=255)
                                                                         tGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCg                                                                         >  1:2969288/1‑150 (MQ=255)
                                                                             cACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTcag                                                                     <  1:3522138/150‑1 (MQ=255)
                                                                              aCTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTGTTAAAGTAAAGATCACCGATCCGTcagc                                                                    >  2:3065113/1‑150 (MQ=255)
                                                                               cTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTcagca                                                                   >  2:1115084/1‑150 (MQ=255)
                                                                                  tCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCg                                                                >  2:4014453/1‑150 (MQ=255)
                                                                                   cAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCgg                                                               <  1:4421609/150‑1 (MQ=255)
                                                                                   cAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCgg                                                               >  2:1734287/1‑150 (MQ=255)
                                                                                   cAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCgg                                                               <  2:2896297/150‑1 (MQ=255)
                                                                                   cAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCgg                                                               >  1:2485487/1‑150 (MQ=255)
                                                                                     gCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGc                                                             <  1:976418/150‑1 (MQ=255)
                                                                                       cTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCAc                                                           >  1:3395216/1‑150 (MQ=255)
                                                                                        tGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACt                                                          >  2:4952160/1‑150 (MQ=255)
                                                                                         gCGAACCGTCTTGCAGACCATTTCCTGTCACAACCGTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACnt                                                         <  1:1026218/150‑1 (MQ=255)
                                                                                           gAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTca                                                       >  2:4571288/1‑150 (MQ=255)
                                                                                             aCCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTcaca                                                     >  2:4934880/1‑150 (MQ=255)
                                                                                              ccGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAg                                                    >  1:4359575/1‑150 (MQ=255)
                                                                                                gTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGAc                                                  >  1:90575/1‑150 (MQ=255)
                                                                                                  cTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACgg                                                >  2:2690775/1‑150 (MQ=255)
                                                                                                  cTTGCAGACCATTTACCTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACgg                                                >  1:3541648/1‑150 (MQ=255)
                                                                                                   ttGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGc                                               <  1:2287935/150‑1 (MQ=255)
                                                                                                       aGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCa                                           <  2:1181496/150‑1 (MQ=255)
                                                                                                         aCCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGCGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATg                                         >  1:1657661/1‑150 (MQ=255)
                                                                                                           cATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGcc                                       <  2:3994607/150‑1 (MQ=255)
                                                                                                           cATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGcc                                       >  2:1943820/1‑150 (MQ=255)
                                                                                                           cATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGcc                                       >  2:4497605/1‑150 (MQ=255)
                                                                                                           cATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGcc                                       >  1:1371812/1‑150 (MQ=255)
                                                                                                             tttACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCg                                     >  1:187630/1‑150 (MQ=255)
                                                                                                                  tttCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGACCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTg                                <  2:720554/150‑1 (MQ=255)
                                                                                                                   ttCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGc                               >  1:2808125/1‑150 (MQ=255)
                                                                                                                     cacaACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAg                             >  1:2057169/1‑150 (MQ=255)
                                                                                                                     cacaACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAg                             >  1:4611446/1‑150 (MQ=255)
                                                                                                                        aacaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGAtt                          <  1:487844/150‑1 (MQ=255)
                                                                                                                        aacaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGAtt                          >  1:3587636/1‑150 (MQ=255)
                                                                                                                             tGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAaga                     >  1:2755300/1‑150 (MQ=255)
                                                                                                                                  aGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTgg                >  2:2375286/1‑150 (MQ=255)
                                                                                                                                   ggCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGa               <  2:4074962/150‑1 (MQ=255)
                                                                                                                                    gCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGaa              >  2:1588294/1‑150 (MQ=255)
                                                                                                                                        aGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGcc          >  1:568064/1‑150 (MQ=255)
                                                                                                                                        aGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGcc          >  2:2973096/1‑150 (MQ=255)
                                                                                                                                         gtgtCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCg         >  1:4486615/1‑150 (MQ=255)
                                                                                                                                          tgtCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGt        >  1:510051/1‑150 (MQ=255)
                                                                                                                                           gtCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTg       >  1:1598654/1‑150 (MQ=255)
                                                                                                                                              cAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCt    <  1:723192/150‑1 (MQ=255)
                                                                                                                                              cAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCt    <  1:2952370/150‑1 (MQ=255)
                                                                                                                                                gTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTc  <  2:4754975/150‑1 (MQ=255)
                                                                                                                                                gTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTc  <  1:3948103/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTC  >  CP009273/3825555‑3825848

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: